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Abstract 



Chiral effective field theory complements numerical simulations of quantum 
chromodynamics on a spacetime lattice. It provides a model-independent formal- 
ism for connecting lattice simulation results at finite volume, and at a variety of 
quark masses, to the physical region. Knowledge of the power-counting regime of 
chiral effective field theory, where higher-order terms of the expansion may be re- 
garded as negligible, is as important as knowledge of the expansion. Through the 
consideration of a variety of renormalization schemes, techniques are established to 
identify the power-counting regime. Within the power-counting regime, the results 
of extrapolation are independent of the renormalization scheme. 

The nucleon mass is considered as a benchmark for illustrating this approach. 
Because the power-counting regime is small, the numerical simulation results are 
also examined to search for the possible presence of an optimal regularization scale, 
which may be used to describe lattice simulation results outside of the power-counting 
regime. Such an optimal regularization scale is found for the nucleon mass. The 
identification of an optimal scale, with its associated systematic uncertainty, mea- 
sures the degree to which the lattice QCD simulation results extend beyond the 
power-counting regime, thus quantifying the scheme-dependence of an extrapola- 
tion. 

The techniques developed for the nucleon mass renormalization are applied to 
the quenched p meson mass, which offers a unique test case for extrapolation schemes. 
In the absence of a known experimental value, it serves to demonstrate the ability of 
the extrapolation scheme to make predictions without prior phenomenological bias. 
The robustness of the procedure for obtaining an optimal regularization scale and 
performing a reliable chiral extrapolation is confirmed. 

The procedure developed is then applied to the magnetic moment and the elec- 
tric charge radius of the isovector nucleon, to obtain a consistent optimal regulariza- 
tion scale. The consistency of the results for the value of the optimal regularization 
scale provides strong evidence for the existence of an intrinsic energy scale for the 
nucleon-pion interaction. 
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Chapter 1 



Introduction 

"One measure of the depth of a physical theory is the extent to which it poses seri- 
ous challenges to aspects of our worldview that had previously seemed immutable." 
(Greene, B. 1999. The Elegant Universe p.386) [IGre99ll 

1.1 Prologue 

The theoretical physicist challenges previous theory, using original research that en- 
ables alternative coherence to emerge, as outlined by Bohm IIBoh92ll (p. 223). The 
underlying theory behind the strong force of particle interactions, which is the force 
responsible for the binding of protons and neutrons together in atomic nuclei, had 
been a persistent mystery throughout the first half of the Twentieth Century. This 
hitherto unknown force acts in opposition to the electric Coulomb force that repels 
positively charge protons from each other, but is at least two orders of magnitude 
stronger at the distance scale of an atomic nucleus. The strong interaction between 
protons and neutrons, or nucleons, is currently most successfully described by the 
theory of quantum chromodynamics (QCD). The advent of the quark model, and 
the theory of the colour force by which the quarks interact, opened a new field of 
research into the internal structure of matter. 

In 1964, Gell-Mann and Zweig independently proposed the existence of a new 
constituent particle, the quark, in order to classify the bewildering array of subatomic 
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particles called hadrons IIGM64L It was discovered that the hadrons can be arranged 
into families that correspond to representations of the group SU(3), and that three 
quark types, or flavours, were required to form the fundamental representation of 
this group. It was not until 1968 that the results of deep inelastic scattering experi- 
ments provided the first evidence of the existence of these new elementary particles. 
As more hadrons were discovered, additional quark flavours were proposed. It is 
currently accepted that six flavours of quark are required to produce the full range 
of hadrons observed in particle accelerator experiments. Their names, in ascending 
order of mass, are: up, down, strange, charm, bottom and top. Of the six flavours, the 
most recent to be discovered experimentally was the top quark, in 1995 at Fermilab, 
with a mass of 172 GeV [|Pro96ll . 

Each quark has a unit of charge equal to +2/3 or —1/3 times the charge of 
a proton (units of +e). An an example, the proton consists of two up quarks and 
a down quark for a total charge of +le, whereas a neutron consists of two down 
quarks and an up quark for a total charge of zero. However, because quarks have 
a certain spatial distribution inside the nucleon, or indeed any hadron, the internal, 
high energy dynamics as described by the behaviour of quarks gives rise to proper- 
ties such as non-zero magnetic moments for the neutron and anisotropic momentum 
distributions. It is clear that in order to describe the internal behaviour of a hadron, 
one cannot assume that a quark behaves as a static source. Instead, the dynamics of 
quarks must be described by a theory, the most successful of which is QCD. 

QCD connects the quark model of nuclear physics to quantum gauge field the- 
ories by introducing the quarks as the relevant degrees of freedom inside a hadron. 
The hadrons are formed by confined colour singlets of three quarks called baryons, 
or quark-antiquark pairs, known as mesons. Quarks are spin- 1/2 fermions, which 
also have the properties of colour and approximate flavour symmetry. Since fermions, 
by definition, must satisfy Fermi-Dirac statistics, the fact that each baryon contains 
three bound quarks in the same state violates the Pauli Exclusion Principle. There- 
fore, it was necessary to suppose an the existence of an additional quantum number, 
known as colour charge, so that each quark may be assigned one of three, orthogonal 
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basis states, labelled red, green and blue. Colour is mediated by the related gauge 
particles of the strong force; the gluons, and the quarks also form a representation of 
the colour gauge group SU(3)c, with eight group generators. 

Mathematically, QCD is a non-Abelian theory. That is, the gauge connection 
of the gluons is non-commutative. The fact that isolated, unbound quarks are never 
found in experiment is one of the striking consequences of a non-Abelian theory. 
Confinement of the quarks within a hadron is a result of the gluon fields exert- 
ing a linear potential that increases with distance between quarks [|Wil74l . That 
is, for a large distance scale, the strong coupling parameter also becomes large. 
This behaviour contrasts with electromagnetism, where the electric Coulomb force 
diminishes as two charged particles are separated. However, quarks experience 
only a small force from the gluon fields as becomes small at short distances 
IIGW73b[ iPoTTSl IGW73a[ IGW74II . This asymptotic freedom is observed when prob- 
ing the internal structure of hadron at high energies, where the small de Broglie 
wavelength of the probe is able to resolve the short distances within the composite 
particle. Near this asymptotically free regime, the methods of perturbative quantum 
field theory are suitable for constructing amplitudes, cross-sections and scattering 
matrices. However, it leads to a difficulty in finding an appropriate method for per- 
forming a calculation with QCD in the low-energy region. Two of the most suc- 
cessful methods that will be discussed in this thesis are chiral effective field theory 
(%EFT) and lattice QCD. 

Using %EFT, one is able to encapsulate the dynamics of a quantum system by 
writing down an 'effective' action of low-energy degrees of freedom. By imposing 
symmetries satisfied by QCD, one can expand out the formula for an observable 
property into a series of quantum amplitudes that can be arranged in order of the 
importance of their contribution by a choice of power-counting scheme: usually in 
increasing powers of mass/energy. These amplitudes can alter, or renormalize the 
calculation of an observable from its naive value, and landmark success has been 
made in confirming these real and measurable effects by experiment. For example, 
the value of the anomalous magnetic dipole moment of the electron agrees with ex- 
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periment to better than twelve significant figures. The Casimir effect (1948), which 
describes the forces arising from the quantum vacuum fluctuations, were success- 
fully predicted by the gauge field theory of quantum electrodynamics (QED). In the 
low-energy, non-perturbative region of QCD, many phenomena can be explained 
by the emergent properties of quark confinement and the behaviour of their bound 
states as hadrons. For example, the proton and neutron also have a large anomalous 
component of their magnetic moment. This is due to the cloud of interacting fields, 
which renormalize the core of the observable. This 'hadron cloud' is one of the 
unique properties of a quantum field theory. Of the available low-energy effective 
theories of QCD, chiral perturbation theory (%PT) is the most notable, in its careful 
incorporation of the fundamental symmetries of QCD. However, the robustness of 
%PT is confined only to a restrictive region called the power-counting regime. Within 
the power-counting regime, the perturbative expansions that occur in %PT are con- 
vergent; the terms of the expansion series are ordered such that higher-order terms 
are sufficiently smaller than lower-order terms. The details of the power-counting 
regime are discussed in more detail in Chapter [3l 

Lattice QCD is a discretized version of QCD, where the dynamics are evaluated 
on a finite-sized box with only certain allowed values of position (or momentum) 
separated by a fixed spacing. Thus, lattice QCD is equivalent to QCD in the limit 
of infinite box size and vanishing lattice spacing. Using lattice QCD, one is able to 
access the non-perturbative, low-energy regime of QCD and provide reliable predic- 
tions of hadronic behaviour. In addition, lattice QCD simulations do not suffer from 
the common problems of quantum field theory associated with renormalization. The 
discrete lattice spacing and the finite box size of the lattice act as an ultraviolet and 
infrared regulator, respectively. Thus, observable quantities evaluated on the lattice 
are finite and calculable. Nevertheless, it can be computationally expensive to evalu- 
ate observables at large box sizes, small lattice spacings and physical quark masses. 
To be able to obtain a result using quark masses as small as their physical values, an 
extrapolation is a practical alternative to a brute-force approach. In addition, the cor- 
rections to finite- volume effects ought also to be calculated for a realistic comparison 
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with experiments. 

1.2 Overview and Aims 

The framework of QCD provides a rich selection of possibilities for inquiry. Among 
these, the low-energy, chiral dynamics of hadrons provides us with a uniquely suc- 
cessful understanding of many of their imporant properties. 

This thesis explores the properties of the aforementioned power-counting regime 
by considering how low-energy constants, which occur in a calculation using the 
methods of xPT, are renormalized, or altered, at different energy scales. This knowl- 
edge of the power-counting regime, in turn, yields insight into the repercussions of 
chiral symmetry breaking in QCD. 

The results of lattice QCD simulations provide an important application for the 
investigation into xPT and the power-counting regime. Lattice QCD results are typ- 
ically produced at a variety of quark masses larger than the physical quark mass. As 
such, a chiral extrapolation to the physical point is required before the result can be 
compared to experiment. In addition, experimental results are not constrained by the 
boundaries of a small box only a few fermi in length. It is important to be able to 
quantify how the finite-volume nature of lattice QCD affects calculations. Analysis 
shows that the finite-volume behaviour of QCD on the lattice can affect the result 
of a calculation in non-trivial ways. Being able to perform an extrapolation that 
takes into account finite-volume effects is also an important step in understanding 
the effects of a finite- volume box on the dynamics of QCD. 

The investigation of the power-counting regime has additional importance. Few 
lattice QCD results in the literature are evaluated at quark masses that lie within 
the power-counting regime. As such, the powerful tools associated with xPT may 
not be used legitimately, since the chiral power-counting expansion of an observable 
would not be convergent. If higher-order terms in the series expansion are not small 
with respect to some power-counting scheme, the result of an extrapolation will be 
scheme-dependent. This thesis describes the construction of an extended effective 
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field theory that can be applied outside the power-counting regime by extracting a 
quantative estimate of the extent of the energy scale-dependence, associated with the 
process of regularization in %PT calculations. It is discovered that lattice simulation 
results themselves can provide guidance on an optimal choice of regularization scale. 
This optimal scale indicates a possible connection with the finite-size of the hadron 
cloud in the form of an intrinsic scale. 

Thus, by analyzing the results from the supercomputer simulations of lattice 
QCD, an intrinsic scale will be discovered that characterizes the finite size of the 
interaction between the hadron cloud and the core of the hadron. 



Chapter 2 
Lattice QCD 



"While the classical vision of the world is intrinsically limited, nothing restricts the 
scientific representation. During the conception stage, the method is free to consider 
all hypotheses, even the most far-fetched, in order to mimic Reality." (Omnes, R. 
2002. Quantum Philosophy : Understanding and Interpreting Contemporary Science 
p.268) IIOmn02ll 

The inception of a discrete, lattice approach to quantum chromodynamics (QCD) 
in 1974 by Wilson marked the beginning of a robust, investigative technique into the 
previously inaccessible low-energy region of strong force interactions IIWil74ll . By 
simulating the behaviour of quarks on a lattice, bound states of hadrons are formed, 
exhibiting confinement, and the behaviour of particle interactions is correctly pre- 
dicted: a testament to the achievement of QCD as a theory of the strong force. 

Lattice QCD provides non-perturbative techniques for obtaining results from the 
low-energy, chiral dynamics of hadrons. It involves the construction of a finite- 
volume box of discrete momenta, with calculations performed from first principles. 
The finite box size of the lattice removes any infrared divergences that would occur 
in infinite- volume QCD, and the lattice spacing acts to regulate the ultraviolet be- 
haviour of observable quantities by limiting the lattice momenta to discrete values. 

In lattice QCD, a Euclidean hypercube is constructed with finite length and dis- 
crete lattice spacing. The quantum field theory can then be represented by the func- 
tional integrals defined on such a box. The momenta can only take the discrete values 
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in the four-box: 

271 

k^t = —n^: (2.1) 
aN 

where a is the lattice spacing, n^j is an integer array representing the lattice sites, 
and A'^ is the number of lattice sites in each direction, such that —N/2 < <N /2. 
Thus, the maximum value can take is 7r/a. This means that the ultraviolet physics 
included in our lattice is entirely determined by the lattice spacing, which thus acts to 
regulate arbitrarily hard momentum contributions to quantum field theoretical quan- 
tities. The real- world dynamics of QCD are recovered in the limit of vanishing lattice 
spacing (the continuum limit) and the infinite- volume limit. 

The dynamics of QCD are encoded in the QCD Lagrangian: a quantity in quan- 
tum field theory that extends the classical notion of the difference between the ki- 
netic and potential energy terms to a density in spacetime. The generalized kinetic 
and potential terms are constructed from the relevant degrees of freedom: quantum 
fields IIWei95ll . The QCD Lagrangian includes a sum of Fermi-Dirac Lagrangians 
for all quark flavours, an interaction term and a Yang-Mills term. In tensor form (and 
summing over repeated indices), the Lagrangian reads: 

^QCD = ^Dirac + ^int + ^YM (2.2) 

= E {%{f -m,)x|/^ - a.xi/^/7y^>^} - \Gl,Cr: . (2.3) 

The fields and are Dirac spinors representing different quark flavours and 
colours, with mass niq. (Dirac spinor algebra was introduced in References [|Dir28b[ 
IDir28aL and some of the basic properties of a Dirac spinor can be found in Appendix 
IA.3I ) The fields G^^ are the non-Abelian field strength tensors corresponding to the 
gluon field J^^, via the equation: 

G^v = - - iasfabc^iK , (2-4) 



where the structure constants fabc are defined in Appendix lA.2[ The Yang-Mills term 
describes the self-interaction of the gluon fields, such that the result is invariant with 
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respect to a special type of symmetry known as the gauge symmetry. In QCD, the 
gauge symmetry is realized in the Lagrangian by forming representions of a mathe- 
matical group, in this case, SU(3)c (where c stands for 'colour'). Each term in the 
Lagrangian must be invariant under transformations involving this group. The quark 
spinors form a basis for the fundamental representation of the group. The gluon 
fields, however, are defined in the eight-dimensional representation of SU(3), and 
the index a runs from 1 through 8. The matrices J a are the generators of the gauge 
group SU(3). A detailed review of the symmetries of QCD is included in Chapter 
|3j Suffice to say, the Lagrangian in Equation (12.31) will be assumed in defining the 
QCD Action in the following Section. 



Lattice QCD relies on a variety of techniques to obtain information about the dy- 
namics of QCD. In particular, the path integral method of quantization serves as 
a starting point, where complex valued Grassmann fields are used to represent the 
quark spinors \|/ and their adjoints (For a short summary on the properties of 
Grassmann algebra and Berezin integrals, refer to Appendix I A.3I ) Before introduc- 
ing the procedure for calculating the expectation values of observables using lattice 
QCD, it is helpful to review the functional methods required for defining the gener- 
ating functionals and the n-point Green's Functions. In the following Section, use is 
made of the functional derivative g^^, the properties of which follow analogously 
from the standard derivative of a function IIRS80L 

Consider the generating functional technique, choosing a set of fields 4> = {^^, \|/, iff/}, 
defined by a set of gauge fields lA'^ and Dirac spinors \|/ & lij/, and integrating over 
all possible paths. In general, for a Lagrangian L (4>,3'"4>), the corresponding action 
can be written as follows: 



2.1 Functional Methods 




(2.5) 
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The generating functional with source terms J7 [xi) takes the form: 
^ i^i)] = * exp |/5[4>] - J d^xj {xi)Mxi) } , l'D^ = flJ d^'' (2-6) 



with normalization: 



(2.7) 



The calculation of the n-point Green's Functions is performed by taking functional 
derivatives of the generating functional with respect to sources y (xi) , and then setting 
each source to zero: 

xW(xi, . . . ,x„) = ^1©$ 4>i ■ ■ ■4>,exp{/S[4)]} . (2.8) 

In order to obtain only the connected diagrams for the generating functional, one can 
define the connected generating functional 'W : 



(2.9) 



The connected (or irreducible) n-point Green's Functions can then be calculated as 
the time-ordered vacuum expectation values of the fields, with respect to the inter- 
acting vacuum | ^1) : 



gW (xi, ■ ■ ■ ,xn) = {a\T [4>(xi) ■ ■ ■4>(x„)] I a) 



1 



(2.10) 



y=0 



The generating functional of Equation (12.61) is useful for constructing an expan- 
sion of amplitudes. This expansion is obtained from the Schwinger-Dyson equa- 
tions, the set of differential equations satisfied by the generating functional: 



1_5_ 



Z[j{xi)]+J{xi)z[3]^0. 



(2.11) 



The Schwinger-Dyson Equations are simply the Euler-Lagrange equations of mo- 
tion for the n-point Green's Functions of the gauge field theory. They provide a 
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continuum persective to the challenging problems of non-perturbative QCD, as sum- 
marized by Roberts and Williams IIRW94II . The investigation of the analytic prop- 
erties of these equations form a crucial component of the study of quark confine- 
ment: where the strong coupling parameter becomes large. The Schwinger-Dyson 
Equations also shed light onto the process of dynamical chiral symmetry breaking, 
discussed in detail in Chapter|3j 

Physical observables of a system can be obtained conveniently using Equation 
(12.101) . To evaluate expectation values (o) numerically, it is common practice to 
remove the difficulties of Minkowski spacetime by an analytic continuation to imag- 
inary Euclidean time, or a Wick rotation, t — )• —it, and S = iSe- Thus the expectation 
values become numerically soluble, since the highly oscillatory behaviour of the 
n-point Green's Functions have been exponentially damped. Thus: 

^ /^4>0exp{-5^[4>]} 
^""^ /2)4>exp{-5^[4>]} ' ^^-'^^ 

which is of the same form as the correlation function in statistical mechanics. Using 
the Euclidean Action, the fermionic part of the partition function can be calculated 
explicitly, leaving an expression in terms of a fermion correlation matrix !M : 

Z = I^Jl^dot{^[R^])oxp{-SE[Ji;,]} . (2.13) 
2.1.1 Wilson Fermions 

In constructing an action on the lattice, such as that of Equation (12.51) . there is a 
difficulty in implementing the fermion field. This difficulty is known as the fermion 
doubling problem. The problem occurs when solving the kinetic part of the Dirac 
Equation of motion, {0 — m)\\r = 0, on the lattice. The derivative d is taken as an av- 
erage (or a forward-backward average so that the result is Hermitian), and the prop- 
agator derived is of the form: sin(^ + m)^^. The correct behaviour of the Green's 
Function is exhibited as p —^0, but as p —^n the propagator also vanishes at the edge 
of the Brillouin Zone: the fundamental cell of a lattice theory with a periodic bound- 
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ary. Thus for sin(|>) = there are 2^^™ degenerate quarks for each flavour, which 
corresponds to sixteen degenerate quarks in four-space. In order to amend this, Wil- 
son introduced a five-dimensional operator, which increases the mass of the doubler 
species proportional to lattice spacing a [|Wil74|| . Note that as a —t- in the con- 
tinuum limit, the Wilson term disappears and does not alter the dynamics of QCD. 
However, the Wilson Action violates chiral symmetry. This important symmetry en- 
sures the consistent renormalization of the low-energy constants of the Lagrangian 
via the chiral Ward Identities, which describe the conservation of a symmetry as ap- 
plied to quantum amplitudes. Chiral symmetry is described in more depth in Chapter 
[3l Additionally, a so-called Clover term is often added to the Lagrangian, which is 
proportional to yifJ'^G'l^y^. This term is also a five-dimensional object, and, like the 
Wilson Action, is suppressed in the continuum limit. In addition, errors of 0{a) 
can be removed, and higher-order errors of 0{a^) can be suppressed by using non- 
perturbatively improved actions IINN951 lLSSW96[ IZ+02l . Lattice QCD simulation 
results relying on a variety of actions are presented in Chapters |4] through |71 and the 
benefits and shortcomings of each one will be addressed as they arise. 

2.1.2 Correlation Functions and the Effective Mass 

Consider the following example regarding the construction of a correlation matrix 
element, and the extraction of the effective mass. In applying lattice QCD to the 
extraction of the mass of the nucleon, one defines interpolating fields % and %, which 
incorporate the structure of a nucleon in terms of its constituent quarks. For ex- 
ample, in the case of a proton, % = £f^^'^{u^Cy^d},)uc is a suitable choice, since the 
maximally anti-symmetric Levi-Civita symbol £ ensures a colour- singlet state, and 
the Dirac spin matrix 75 (defined in Appendix lA.ll) preserves the spinor properties 
of the interpolating field. The fields u, d are Dirac spinors representing the up and 
down quarks, respectively, and the charge conjugation matrix C = zyoY2 ensures that 
the product of a spinor and its transpose satisfies Lorentz invariance. 

The two-point Green's Function for a proton, or more generally, the nucleon, can 
be expanded by inserting both a complete set of momentum- and spin-dependent 
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eigenstates \A{q,s)), and a translation operator on the % field: 

G(2)(x,0 = (n|x(x,05C(0)|^2), (2.14) 

G^^\p,t)^ £ e-'P-'{a\e-'^-^X{0)e''i-^\A{q,s)){A{q,s)\xiOm (2.15) 

= e-'MV^-^®^(n|x(0)|A(^,5))(A(^,5)|x(0)|n) (2.16) 

= L Hp-q)e-^'^^'mxm\Mq,s))\' (2.17) 

A,q,s 

= Ll^A,p,.| V^^(^Mp,*)¥(p,^), (2.18) 

for complex-valued scalar coefficients XA,p,i. ^1 p * ^^'^ spinor fields vfr defined by 
the matrix elements: 



Xa,p,.X|/= (r2|x(0)|A(p,.)); XXp^,^={A{p,s)\x{0)\^). (2.19) 

The mass of the nucleon can then be extracted from the two-point Green's Function 
at zero 3 -momentum, that is, Ea{p = 0) = Ma. To obtain a measure of this quantity 
from the exponential, one defines the effective mass Mgff by comparing the behaviour 
of the Green's Function at times t and t+l: 

Note that Mgff is a dimensionless quantity, and the calculation of the mass of the 
nucleon must involve the conversion to physical units from lattice units by dividing 
by the lattice spacing a. Since the Green's Function incorporates the full quantum 
mechanical spectrum of modes, the behaviour of Mgfr is strongly influenced by the 
excited states of the nucleon at small t. In the limit of large t, however, the ground- 
state nucleon mass can be recovered: 

Mn = lim (2.21) 
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2.1.3 Quenching and Computational Alternatives 

Computing the quantity det(5Vf [^^]) is the most time-consuming operation in the 
calculation of the partition function in Equation (12.131) . For this reason, calculations 
are performed at fermion masses larger than their physical value, thus decreasing 
the Compton wavelength of a fermion and significantly reducing the computational 
resource requirement of the summation over all paths and the time required to ex- 
ecute all the necessary fermion matrix inversion algorithms. Usually, results from 
lattice QCD are obtained at multiple fermion masses, so an extrapolation can be 
used to obtain the result at the physical value, or at zero mass (the chiral limit). A 
complementary computational simplification known as quenching exists, whereby 
det{M [^^]) is set equal to a constant. This has the effect of removing from the 
theory all vacuum polarization diagrams, changing the dynamics of the quantum 
field theory in a non-trivial way. For this reason, quenched QCD (QQCD) should 
be considered, in essence, a different theory from QCD. The results from QQCD 
calculations can nonetheless be interesting points of investigation, as they offer a 
unique testing ground for extrapolation schemes. This is because results from the 
unphysical QQCD calculation cannot be known in advance from experiment. 

Several other alternatives to quenching have been used in the literature to date. 
Sometimes, the vacuum polarizations, normally omitted in QQCD, are calculated 
for a different (usually larger) quark mass than the valence quarks, which couple to 
external sources. The quarks that appear in the disconnected loops are known as sea 
quarks. This distinction between sea quark mass and valence quark mass provides 
some of the dynamics of QCD, albeit altered, whilst still ameliorating the computa- 
tional intensity of the calculation of det(fW [.^^]). An alternative, particularly used 
in electromagnetic contributions to QCD, is to omit diagrams that include indirect 
couplings, that is, external fields coupling to sea quark-antiquark pairs, as shown in 
Figure 12. 1[ The computation of the indirect couplings to disconnected quark loops 
is by far the most time-consuming portion of the calculation of a diagram. Valence 
QCD (VQCD) therefore only includes diagrams where any external particles, such 
as incoming photons, couple directly to valence quarks in the relevant hadron. Al- 
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Figure 2.1: An external photon coupling to a sea quark-antiquark pair. Diagrams including this 
kind of coupling are omitted in Valence QCD. 

though the resulting theory differs from full QCD, often properties of particles are 
calculated using an isovector combination. In the case of an isovector, a linear com- 
bination of isospin partners is formed so that the resultant combination transforms 
as a vector in isospin space. For example, in the case of the nucleon, the combi- 
nation of the fermion fields: p — n (proton minus neutron) is isovectorial with total 
isospin of 1, and all diagrams containing indirect couplings cancel. This is because 
diagrams that contain indirect couplings to disconnected loops are exactly the same 
for the proton and neutron, and thus disappear in the combination: p — n. It is only 
the valence quark composition that differs between the proton (uud) and the neu- 
tron (ddu). Thus the distinction between full QCD and VQCD disappears for this 
observable, and the calculations of its properties are less computationally intensive. 



2.2 Lattice QCD Applicability and Issues 

It is important to identify clearly the constraints of lattice gauge theory. Lattice 
QCD is well defined over all box sizes, lattice spacings and quark masses, and it 
is also infinitely scalable. However, the computational cost of the calculation of an 
observable is generally proportional to the square of the lattice volume and inversely 
proportional to the sixth power of the lattice spacing. To avoid major finite-volume 
effects, the literature suggests that the lattice box length should be about 2.5 to 3.0 fm 
[ISW85[ ILTTWOOi IFKOU951 IDLL96[ lLS96]l . This is the typical size of most current 
lattice QCD calculations. Nevertheless, finite-volume effects can still be significant 
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at these box sizes, and ought not to be neglected. In fact, for many observables, a box 
length of 3.0 fm is insufficient to avoid large finite-volume corrections at physical 
quark masses. This will be demonstrated in Chapters H] through |7] for a variety of 
observables. 

While continued supercomputing advances in numerical simulations of lattice 
QCD are important, one ought to recognize its limitations in providing a thorough 
understanding of the internal structure of hadrons, which can be aided, in part, 
by complementary techniques such as chiral perturbation theory (xPT). For exam- 
ple, consider the effects of the mesons known as kaons, vital for understanding 
strangeness in the nucleon, which appear in the meson octet (see in Appendix I A. 41) . 
One must either use xPT calculated to significantly high order in the relevant per- 
turbative expansion, or develop new non-perturbative approaches which utilize the 
non-perturbative information expressed in the lattice simulation results. Since the 
former is likely to be compromised by the asymptotic nature of the expansion, atten- 
tion is given to the latter approach in Chapter HI 

The computation of observables in lattice QCD provides great insight into the 
non-perturbative region of QCD. As long as one can account for finite-volume and 
momentum discretization effects, lattice QCD provides excellent predictions of the 
behaviour of quarks at low-energy. In simulating the interactions of hadrons, and 
demonstrating confinement, lattice QCD is a landmark achievement in the realm of 
chiral dynamics. 

The complementary methods obtained from effective field theory offer guidance 
in the calculation of observables on the lattice. They provide estimates of finite- 
volume effects and extrapolations to physical quark masses, and providing a deeper 
understanding of the applicable regions of lattice QCD. This can serve to ameliorate 
the otherwise unseen difficulties encountered in a brute-force approach to calculation 
by considering symmetries, renormalization, power-counting, and other techniques 
built into the formalism of chiral effective field theory. Each method presents its 
own challenges, but also brings enlightenment through the significantly different 
approaches to a given problem. 



Chapter 3 

Chiral Effective Field Theory 



"Everything can be tried, a bold abstraction of something that has succeeded else- 
where, the exploration of the faintest clue, or a leap through empty spaces. . . 

Thus, the method exists, boundless, its ultimate foundation being the freedom of 
the mind." (Omnes, R. 2002. Quantum Philosophy: Understanding and Interpreting 
Contemporary Science p. 268) [|Omn02ll 

In an effective field theory, one identifies the relevant degrees of freedom at a 
particular energy, and encodes the behaviour of these degrees of freedom in a suit- 
able Lagrangian. For a low-energy effective field theory corresponding to quantum 
chromodynamics (QCD), such as chiral effective field theory (xEFT), effective La- 
grangians may take on different forms to the QCD Lagrangian, but the physics of 
the strong interaction must remain the same in each case. That is, results for cal- 
culating elements of the S-matrix must agree among effective field theories, up to 
some order. In order to construct such a theory, terms in the effective Lagrangian are 
chosen so that they satisfy the fundamental symmetries of QCD. The coefficients of 
the terms in the effective Lagrangian are new coupling constants, the values of which 
are determined from experiments. 

The method of effective Lagrangians provides alternative machinery to lattice 
QCD for understanding the low-energy behaviour of QCD, and physical theories in 
general at a specific energy level. The dynamics of the low-energy degrees of free- 
dom, such as mesons and baryons in the case of xPT, are incorporated directly into 
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the Lagrangian, whereas very massive particles are treated as static sources [|DGH96[ 
IBor07L Examples of important effective field theories include the Sigma Model 
and its various instructive representations, the MIT Bag Model [|CJJ^74ilJoh78ll and 
Cloudy Bag Model HMTTBOll . as well as quantum electrodynamics (QED) and QCD 
themselves IIWei95i 

Recall that the QCD Lagrangian comprises a Yang-Mills term involving vector 
potentials R^^, their field strength tensors G^y = d^Jl"-^ — iasfabc-^j^-^v ^^'^ ^ Dirac 
term of quark spinors \\f corresponding to a mass matrix 5Vf . The spinors and the 
mass matrix are extended to contain the six flavour and three colour components 
of QCD. Using the slash-notation y"/)^ = l/>, the QCD Lagrangian may written out 
conveniently in matrix fomJl]: 

^ 1 _ _ 

LQCD = Wi^-^)^-2^AGMV&''], (3.1) 

where the trace acts over colour indices for the matrix-valued versions of the gluon 
field strength tensor G^v, defined by summing over the generators J" of SU(3): 

G^v = ^''G^v = ^\p^v] + '«.[-S^„ -^v] ■ (3.2) 
The generators J" in the eight-dimensional representation of SU(3) are related to the 



Gell-Mann matrices X", defined in Appendix lA. 21 by a factor of a half: 

The quark-gluon interaction vertex is incorporated into the covariant derivative, de- 
fined as: 

D^ = a^ + /a,y.^;, (3.4) 

which acts as a parallel transport in gauge- space, so that the QCD Lagrangian of 
Equation (13.11) is gauge-invariant. By substituting into Equation (13.21) . it can be seen 



The double-headed arrow indicates the difference between the derivative acting to the right and 
to the left. i.e. D = D D. 
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that the anti-symmetrization of the covariant derivative is the field strength tensor. 
This is a consequence of the gauge connection lacking torsion: 



The fundamental symmetries of QCD are built into the QCD Lagrangian of 
Equation (13.11) . In particular, chiral symmetry will be important in the subsequent 
analyses of observables using %EFT. The consequences of chiral symmetry break- 
ing ultimately have a profound effect on the behaviour of subatomic particles, their 
masses, magnetic moments and other properties. Therefore, it will be beneficial to 
describe some of the subtleties of chiral symmetry with care in the discussion below. 



In general, a symmetry, or an invariance of a dynamical quantity under a transfor- 
mation of one of its parameters, leads to important physical insights into a system. 
Noether's Theorem demonstrates that a conserved current can always be constructed 
from a (non-anomalous) symmetry of a field theory. 

Chirality is defined as the handedness of the representations of the Poincare 
group (which encodes the isometries of Minkowski spacetime) under which the 
quark spinors transform. It is related to the helicity of a particle: the projection 
of its spin on its direction of linear momentum, which is equivalent to chirality if the 
quarks are massless. Helicity is not in general a Lorentz-invariant quantity. Its value 
in one frame may be flipped with respect to its value in a boosted frame. 

The QCD Lagrangian in Equation (13.11) can be split into separate left- and right- 
handed chiral states under the projections Tij^ = j{l ±75). The left- and right- 
handed spinors are written as: 



(3.5) 



3.1 Chiral 




¥L,i? = rL,/?V|/. 



(3.6) 



Note that the resultant chirality of the quark fields is decoupled only for zero mass 
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IIDGH96II : 

mass 

= /% + m - 2 Tr [G^v G^"] - ( Vl + ^ ¥l) • (3.7) 

The quark fields transform under the chiral rotations L and R, which are elements of 
the left- and right-handed Lie Algebra, respectively, defined for the group generators 
Qli^ and arbitrary, continuous, real parameters oc^^: 

L = exp(/<e£)e5f/(3)L, (3.8) 
R^cxp{ia],Q%)eSUi3)R. (3.9) 

The transformation laws for each of the spinor fields can thus be written: 

Vl Lm^l = Vl + 5\|/L , (3.10) 

^R^R^R = ^R + hM^R. (3.11) 

Noether's Theorem allows one to construct the left and right symmetry currents, with 
the corresponding time-independent charges forming the eight unique invariants of 
the group. These invariants are the generators, and are found by integrating over 
a spacelike surface a. Note that in the case of S\]{3)l,r the generators are related 
to the previously mentioned Gell-Mann matrices X", after chiral projection by the 
group elements (up to a minus sign and a factor of a half, by convention): 

QLr = jA<5,Jl% = = -Tl.r^. (3.13) 

An equivalent convention to that of left/right chirality is the construction of vec- 
tor and axial vector transformations. The group action can be written out explicitly 
for either convention, using the definition of a Lie group with continuous group pa- 
rameters tty^. The charges Qy and 2^ simply count the sum and the difference of 
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left- and right-handed fermions, respectively: 



y = exp(/a^(2f.) eSU(3)v, 
A = exp(/aiel)eSU(3)A, 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



These sets of rotations are the most convenient for asserting the invocation of an 
important theorem known as Goldstone's Theorem. Goldstone's Theorem, described 
below, is crucial in understanding the connection between axial charges and the 
origin of mesons in QCD. 

3.1.1 Spontaneous Symmetry Breaking 

In QCD, particles are believed to utilize the Nambu-Goldstone mode of spontaneous 
breaking of a continuous global gauge symmetry. This symmetry breaking occurs 
in flavour space, and only the lightest three quark flavours will be considered: up, 
down, strange. Since the up and down quarks have relatively low mass (m„ ~ 2-6 
MeV, nis ~ 100 MeV) compared to the other quarks (m^ ~ 1.3 GeV), they contribute 
the most strongly to symmetry breaking effects. 

Goldstone's Theorem states that the symmetry group SU(3)v ® SU(3)a is not 
respected by the (no-particle) vacuum state |0), even though this group is a symmetry 
of the massless QCD Lagrangian. One might naively expect that the vacuum state is 
invariant under the group transformations: 



y|0)=A|0)=|0). 



(3.18) 



Chapter 3. Chiral Effective Field Theory 



Hall 22 



Noether's Theorem entails that the symmetry charges are time-independent: 

IQv,a=0=Wqcd.QU- (3-19) 

This means that the charges should annihilate the vacuum state |0), since the QCD 
Hamiltonian i^QCD annihilates the vacuum state. In 1984, Witten and Vafa proved 
this result for vector charges even without assuming chiral symmetry [IVW84II . How- 
ever, were this the case for axial charges, a spectrum of mass-degenerate partners 
with opposite parity would be expected to exist for all hadrons. This is because 
the axial charges are odd under parity transformations, and any state acted on by 
the axial charges will also retain the same energy eigenvalue (but with a flipped 
parity eigenvalue), because of the commutation relation in Equation (13.191) . There 
is a stark lack of experimental evidence for such particles llY"'"06ll . Thus, physical 
hadrons merely observe the symmetry group SU(3)v/. 

Instead of annihilating the vacuum state, the axial charges transform it to an 
element of a new Hilbert space: 

e^|0)=0, (3.20) 
Q"A\0)=qA\n%p))y^O. (3.21) 

This new state (with axial eigenvalue qA) has the same energy as the vacuum state as 
long as the symmetry is not also explicitly broken by terms in the Lagrangian. Gold- 
stone's Theorem states that new particles are created, the number of which corre- 
sponds to the number of generators for the relevant representation of SU(3)a. These 
new particles must be massless and spinless pseudoscalar mesons, called Goldstone 
bosons. 

If the physical manifestation of a symmetry of a Lagrangian involves the spon- 
taneous breaking of one or several local continuous transformations, the theory pre- 
dicts a massive spin zero boson called a Higgs field, and the Higgs mode is said to 
be realized. Although the Higgs mode is not expected to occur in the strong nuclear 
force sector of the Standard Model, its actualization in the electroweak sector would 
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result in a mass term for the {W ,Z} weak gauge bosons. Such a mass is observed 
in experiments, and also explains how the charged fermion fields gain mass, through 
the following argument. By considering the Lagrangian for an SU(2) complex dou- 
blet of bosons, which can be expanded about its minimum potential energy in the 
same manner as the Goldstone bosons, one must arbitrarily choose a direction in 
isospin space in which to expand. Three of the Higgs degrees of freedom combine 
to become the longitudinal spin modes of the three weak gauge bosons, and the mass 
of the fermions is produced by the vacuum expectation value of the remaining Higgs 
boson, which remains in the theory IIGHK64L It should be noted that the Higgs 
mechanism contributes only a small amount to the mass of hadrons in QCD, and 
that the dominant process for their mass acquisition is dynamical chiral symmetry 
breaking [|GN74llRCR10ll . A more detailed analysis of the consequences of dynam- 
ical chiral symmetry breaking for the mass of the nucleon is discussed in ChapterHl 
in the context of varying quark masses in lattice QCD results. 

3.1.2 Partial Conservation of the Axial Current 

Before discussing the powerful techniques associated with effective Lagrangians, a 
brief overview is now presented for the current algebra method for obtaining the 
low-energy matrix elements of pion decay. It is known that the SU(3) axial currents 
7^" are non-zero. But in order to know exactly how these matrix elements vary and 
how they depend on the octet meson masses, one requires a current algebra tech- 
nique known as Partial Conservation of the Axial Current (PC AC). The statement of 
Goldstone's Theorem in Equation (13.211) can be re-expressed as a matrix element: 

(0|7^|7^^(p))=/^;.^6^^ (3.22) 

from which follows the divergence: 



(3.23) 
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Equation (|3.22l) serves as a suitable definition of the pion decay constant /71. Taking 
the value from experiment: ^ 92.4 MeV. 

Equation (13.231) . together with the Haag Theorem, forms the principal statement 
of PCAC, that either tt^ or 3^/^" can be used equivalently, and that if the pion mass 
becomes zero then the axial current is totally conserved. Thus the following relation 
may be written: 

:i'' = v^a^/r- (3-24) 

This situation is a special case of the Soft-Pion Theorem for a matrix element in- 
volving a general local operator : 

lim {n%p)^\ o \a) = [Jl 0] \a) . (3.25) 

While they are useful in obtaining specific information about the low-energy matrix 
elements of pion decay, the methods of PCAC can be subtle in determining possible 
momentum dependence in an amplitude of a low-energy process. One must also 
make the assumption that matrix elements vary continuously in taking the soft pion 
limit, — !• 0. The method of effective Lagrangians is less awkward in obtaining the 
appropriate momentum dependence and any quantum corrections to a low-energy 
amplitude. This is because the effective Lagrangians are ordered by a systematic ex- 
pansion in momentum or mass, which encodes the relative importance of corrections 
to an amplitude in question. 



3.1.3 The Sigma Model 

The Linear Sigma Model IIGML60II is a useful pedagogical tool, because with it im- 
portant theoretical techniques such as the construction of symmetry currents, spon- 
taneous symmetry breaking and changes in parameterization can be demonstrated 
easily [IDGH96II . First, consider an SU(2) Sigma Model Lagrangian consisting of a 
massless spinor field \|/, a so-called pion field S spanning the triplet representation of 
SU(2) and a massive scalar field o. The consideration of isospin symmetry in SU(2) 
provides a simple and instructive example for investigating symmetries [IDGH96L 
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The Lagrangian takes the following form: 

11 ^ X 

L„ = ^fi'i\\f+-d^n-d^Ti+-d^od^o-g^{o-n-ny5)\\f+^{o^+n^^ 

(3.26) 

(for constant coupling parameters g, /j and X, and SU(2) Pauli Spin matrices x de- 
fined in Appendix lA.il) . 

Spontaneous symmetry breaking occurs in the Lagrangian of Equation (|3.26l) for 
> O.ln minimizing the potential: 

V(0,7r) = -y (O^ +7t2) + ^(o2 + ^ (3.27) 

a ground state is found that is non-trivial (unlike the case < 0, for which the only 
ground state solution is: o = S = 0). This ground state is defined by: 

o^ + Te = ^. (3.28) 

By redefining the a field and expanding the Lagrangian about the new ground state 
(a)o = V, the Linear Sigma Model exhibits spontaneous symmetry breaking, evident 
in the acquisition of mass for the d field: 

6 = a-v, (3.29) 
La = - gv)\\f+ ^d^n ■ a^'S -f ^ (a^da^o - 2p^d^) - g\\f{d - IX ■ Sy5)\^ 
- Xvd(d2 + S^) - ^ [(a^ + n^f- v^] . (3.30) 

Nevertheless, SU(2) isospin symmetry is preserved in this Lagrangian. 

The active degrees of freedom in an effective field theory do not necessarily cor- 
respond to elementary particles of nature, and so it is expected that changes in the 
representation do not alter the outcome of physical processes. This notion is formal- 
ized in the Haag Theorem flHaaSSI . which states that for two field variables derived 
from (unitarily) equivalent representations, if one is a free field, then the other is 
also free, regardless of how they are related and whether the associated diagrams 
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and Lagrangian vertices change IIGue66ll . As a corollary, an interacting quantum 
field theory 'does not exist', in the sense that its fields do not transform covariantly 
under the interacting Poincare group. Weinberg suggested that only free fields are 
required to construct the S-matrix from the relativistic Hamiltonians in QED, but in 
QCD one must simply resort to writing down the most general Lagrangian [|Wei95i 
An alternative approach demonstrates that an interaction picture can be constructed 
consistently if time evolution is taken to be only locally unitarily implementable 
IIGue66L 

By redefining the scalar field in either linear or non-linear combinations of the 
other involved fields, different sets of interaction vertices can be assembled. For 
example, using the Linear Sigma Model, two particularly instructive representations 
are considered for later adaptation to low-energy QCD. By rewriting the heavy a 
field and pion triplet as a matrix quantity E = a + /x ■ S, the resultant new field E 
transforms as an object in the adjoint representation, which forms left cosets of the 
group SU(2)l ® S\J{2)r, as described by Scherer [ISS05H : 

L^LLR\ (3.31) 

In this representation, the Lagrangian becomes: 

= w'i^L + w'i^R + ^Tr [a^Ea^E^] + i^^Tr pL] - ^Tr pL]^ 

- g (^fLY\\fR + V/^eVl) . (3.32) 

This form is useful because it allows one to identify easily the terms involved in 
spontaneous chiral symmetry breaking (~ Tr[E^E]). Terms responsible for explicit 
chiral symmetry breaking (eg. o ~ Tr [E -f E^]) do not occur in this case, but will be 
considered in the context of xPT, in Section [X2l 

The exponential representation is most commonly employed for its application to 
low-energy QCD. By defining a matrix- valued field U = exp(zx ■ S/v) that transforms 
the same way as the previous E field, and a massive scalar field S, the Lagrangian 
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becomes: 




Tr[a^f/a^[/'^] 



(3.33) 



for an arbitrary coupling constant v. This representation combines the matrix form 
with a heavy scalar degree of freedom, which can be integrated out of the theory eas- 
ily using the prescription provided by Donoghue, Golowich and Holstein [IDGH96L 
This is exactly the form needed to construct a low-energy effective field theory. 



The formalism of chiral perturbation theory (%PT) will take advantage of Gold- 
stone's Theorem and the study of symmetries discussed in the previous section. In 
this case, however, the global gauge group considered is flavour SU(3). In order 
for the effective field theory to emulate physical results, one must write down the 
mechanics of a Lagrangian field theory incorporating the necessary symmetries and 
degrees of freedom at the observed scale. To represent particles such as pions and 
kaons obeying Bose-Einstein statistics, one can write the standard massless scalar 
Lagrangian: 



and interpret Tl" as the octet of Goldstone bosons (whose explicit form can be found 
in Appendix IA.4I) . By defining a matrix- valued function U , and its transformation 
law, one can collect together the interaction terms in the exponential representation 
in a similar way to the Sigma Model: 



3.2 Chiral Perturbation Theory 




(3.34) 




(3.35) 



U ^LUR\ 



(3.36) 
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with constant /. Now the effective Lagrangian can be written down as an expan- 
sion of successive orders of momenta. The two derivatives in the scalar Lagrangian 
mean that only even chiral powers are admitted for particles such as mesons. For 
the lowest-order free mesonic Lagrangian, there is only one low-energy coupling 
constant, /: 

^iff=^Tr[a^t/a^[/t]. (3.37) 

(Higher-order mesonic Lagrangians can be found in References i SS05llBor07llBer08ll .) 
The coefficient / from the definition of the field U appears here as a low-energy 
constant (LEC), since it is expected that the expanded effective Lagrangian for the 
pseudo-Goldstone fields will have the standard normalization for bosons, Zeff = 
^dfjTt^d^n'^ + O (tt"^) . This LEC can further be identified with the pion decay con- 
stant /ti by first considering the Fermi weak interaction Lagrangian as a left-handed 
source field and computing the decay rate from the resultant invariant S-matrix &\&- 
ment IIDGH96[[SS05l . 

The second-order Lagrangian of Equation (13.371) will be the starting place for the 
consideration of the low-energy meson sector of QCD. 

3.2.1 Meson Sector 

In the theory of mesons, one considers a set of Goldstone boson fields and interprets 
them as the meson sector of QCD. One can use the knowledge of explicit symmetry 
breaking from Section [J. 1.1 1 to provide the fields with a (small) mass. Using the ex- 
ponential representation, U {x) can be systematically expanded in powers of its small 
momentum and mass with respect to some energy scale A^^. In 1984, Manohar et al. 
identified this scale of chiral symmetry breaking as Aj^ ~ Anf ~ 1 GeV [IMG84II . In 
renormalization, this is the scale at which the next-order loop contribution retains 
the same effective coupling strength (see Section [33l) . 

The total mesonic Lagrangian can be written out in the expanded form of even 
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chiral powers HBorOVII : 

oo 

L^{U{x),M ) = £ i:f ) {U{x),M). (3.38) 

i=l 

In order to quantify the extent of the chiral symmetry breaking caused by the mass 
terms in the expansion, initially let M transform as a field {M — )• LM K^), so that 
the Lagrangian will remain invariant under global SU(3)l ® SU(3)^. At lowest non- 
trivial order: 

, (2) , kin I I- mass 

- ^Tr [d^Ud^U^] + ^Tr [MU^ + Um"^] , (3.39) 

where Bq is a constant (with dimensions of mass) included for generality. Chiral 
symmetry breaking then results from imposing the Hermitian condition for the quark 
mass matrix M = Thus the constant Bq directly corresponds to the extent of 
chiral symmetry breaking [|SS05llBor07L 

Some terms in the Lagrangians of either QCD or %PT explicitly break chiral 
symmetry. For example, i^mass involving the quark mass in Equation (13.71) is invari- 
ant under an axial group action A = exp(— ^a^X^ys). The associated axial Noether 
currents J'^" encountered in PCAC will not be conserved, but diverge according to 
the equation: 

a^yf = 2iWY5^v. (3.40) 

To relate the meson masses to the quark masses, consider chiral SU(3). It is 
expected that the vacuum expectation values of the scalar quark densities are the 
same in each theory: QCD and %PT. That is, the quark condensate (qq), where q 
stands for u, d or s quarks, should be an observable independent of representation. 
Consider the explicit chiral symmetry breaking terms of each Lagrangian, namely, 
-^mass defined in Equation (13.71) for QCD, and the mass term i of Equation (13.391) 
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for xPT. By expanding out the exponential field U in one obtains: 

J^r' = BoflTv - ifioTr [M n^] + (tt^) . (3.41) 

For approximate isospin symmetry Mu ^m^^ j {niu -\-md) = m^ms, expanding the 
first of these terms yields the relations: 



r) r r mass 

(W) = -(0|%^|0) = -(0|%-|0) = -B,fl (3.42) 
am am 



Thus there is a profound connection between quark condensation and the process of 
dynamical chiral symmetry breaking. The second term yields the Gell-Mann— Oakes— Renner 
Relations [IGMOR68II relating meson masses to quark masses: 

ml = 2Bom, (3.43) 
m| = 5o(mrf + m,), (3.44) 

m\ = -Bo{m + 2ms) ■ (3.45) 

Just like PCAC, Equation (13.431) shows that if the light quark masses are zero, then 
the pion mass must also be zero, and thus chiral symmetry holds. This leads to the 
Gell-Mann— Okubo Mass Relation: 



3m^ = Am\ - ml . (3.46) 



By additionally enforcing local chiral symmetry, the set of all chiral Ward Identi- 
ties become an invariant of the generating functional encountered in Section (|2.1I) . as 
long as no anomalies are present nSSOSII . The chiral Ward Identities simply encode 
the statement of symmetry preservation and the existence of conserved quantities as 
a consequence, much like Noether's Theorem; but applied to quantum amplitudes. 
Consider QED, a U(l) gauge theory, as an example. The Ward Identity amounts 
to a statement of charge conservation, and the existence of a conserved electric cur- 
rent. In QCD, to be able to generate all the Green's Functions for the theory, the 



Chapter 3. Chiral Effective Field Theory 



Hall 31 



Lagrangian must include pseudo-scalar (p) and vector (/^, r^) source fields, which 
vanish to recover the standard QCD Lagrangian in Equation (13.11) . and a scalar field 
(s) that assumes the role of the quark masses 5Vf . This is known as the method of 
external sources. This generalization of the QCD Lagrangian is vital for calculating 
the divergence of Green's Functions. These fields obey the following transformation 
laws for the local chiral rotations L{x),R{x) E SU{3)i^^r: 

L{x) l^L^ix) + i{d^L{x))L\x) , (3.47) 
r, ^ R{x) r,R\x) + i{d ^R{x))R\x) , (3.48) 
{s + ip) -> L{x) {s + ip)R^ {x) . (3.49) 

The QCD Lagrangian, invariant under local SU(3)l ® SU(3)^, becomes: 

^QCD = - ^Tr [G^vG^"] - + ip)^R -^r{s + ip)^L 

- WA-rL/^¥ - W^r^^¥ ■ (3.50) 

In the case of the low-energy effective Lagrangian, one must define a covariant 
derivative with transformation law: 

V^U = dfJJ + il^U -iUr^, (3.51) 
V^U -> L{x)V^UR\x) . (3.52) 

Therefore, the lowest-order non-trivial Lagrangian for mesons obeying local chiral 
symmetry can now be written with mass source defined using the convention % = 
2Bo{s + ip), functioning as a field, as before: 

^(2) ^ J^^^ [V^UV^U^ + XU^ + f/%+] . (3.53) 
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3.2.2 Baryon Sector 

Since the Lagrangian of a low-energy theory can be expanded out in a convergent 
series of small momenta p/Aj^, the mass of the baryons themselves cannot be treated 
as an expansion parameter, since their mass and momenta are of the same order of 
magnitude as the scale A^^; thus the perturbation theory diverges. That is, the mass 
of a baryon Mb ~ A^^, so the expansion parameter Mb/A-^ cannot be small. To 
overcome this difficulty in ordering the chiral series in the baryon sector of %PT, 
consider the heavy-baryon approximation. 

Define some alternative fields By{x) to the SU(3) octet baryons B{x) = B"{x)'Xf^, 
with velocity v^ largely unchanged by pion interactions [|Geo90l |JM91a[ |JM91bL 
These fields By{x) are only just off-shell by a small amount k ■ v: 

p'^=MBV^ + kf. (3.54) 

A perturbation theory about this small momentum can now be constructed. In 
addition, the difficult spin structure of the new fields By can be handled by using 
the particle projection operator = j{l +f), thus absorbing the effects of virtual 
baryon loops into higher chiral orders of the theory: 

5vW =i've'^*^'"'5(x). (3.55) 

This procedure can be repeated in exact analogy for the totally symmetric Rarita- 
Schwinger tensor t^"''^ (^x), which represents the spin-3/2 decuplet fields, as long as 
all spin-1 /2 components are removed (y- Ty"^*^' = 0). It is defined by: 

Ty{x) = /'ve'(^^+^^)'*"-^r(;c) . (3.56) 

The sum of the octet and decuplet masses is used, by convention, in the exponential 
in order to avoid extra factors of mixed octet-decuplet fields in the final Lagrangian. 
This results in a positive term proportional to the mass splitting A = \Mt — Mg| 
njM91all . (The explicit representation of these fields in SU(3) can be found in 
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Appendix IA.4I) . Treating the mass splitting A ^ Aj^ as a small perturbation, the 
new velocity-dependent fields and (indices suppressed) obey a massless Dirac 
Equation, and a Dirac Equation with small mass splitting, respectively: 

0By{x)=O, (3.51) 
{0 - A)Ty{x) = . (3.58) 

To write out a completely velocity-dependent Lagrangian for baryons and their 
interactions with mesons, it now remains to rewrite all Dirac bilinears in terms of 
a covariant spin operator 5v = — |Y5[Y", Y^]vv, which has the useful property that 
its commutation and anti -commutation rules depend only on the four-velocity v^. 
The meson interactions are incorporated into the theory by coupling baryon fields to 
the axial current encountered in PCAC (Section (13.1.11) ). which is equivalent to the 
Goldstone bosons as per the Haag Theorem. The convention is to define exponential 
fields = i7, which follow the transformation rule [|JM91bll : 

^^L^H\x)=H{x)^R''. (3.59) 

The transformation matrix H = H{x) is a spacetime dependent combination of the 
chiral transformation matrices and the Goldstone bosons themselves. This means 
that the octet and decuplet fields' transformation rules also involve H, and in fact, 
the axial current and the octet baryon field By are exactly analogous to the ^ field 
in their transformations. The additional subtlety with the decuplet field is that each 
of its three indices transforms separately: 

B^HBH\ (3.60) 

jabc _^ fjaa'fjbb'fjcc'ja'b'c' (g ^j) 

Because the transformation matrix // is a spacetime-dependent object, a vectorial 
connection needs to be included to preserve the gauge invariance of the Lagrangian. 
Similarly, an axial vector combination of exponential fields can be defined. Under 
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the Haag Theorem, these axial vectors are equivalent to the pseudo-Goldstone boson 
fields: 

y^ = l(^a^r+r5,^), (3.62) 

V^^HV„H'^-{d^H)H\ (3.63) 
A, = ^i^d,^^^-^,%Q, (3.64) 
A^^HA^H^. (3.65) 

Thus the covariant derivative can now be included for both octet and decuplet fields. 
As before, the decuplet requires a separate connection to act on each index: 

©^5v=a^5v + [y^,5v], (3.66) 
^^j^abc ^ ^^jaabc ^ ydjadbc ^ yd^j^adc ^ ydjaabd _ ^^^^-^ 

The most general lowest-order Lagrangian for the baryon octet and decuplet 
fields, including transition vertices, can be now written by identifying the relevant 
SU(3) invariants [IJM9 1 ai lJM9Tbl lJen92l lLS96l IWL051 IWLTY07H : 

^oct&dec = /Tr[5v(v- ©)5v] +2DTr[5v5^{A^,5v}] +2FTr[5v5^;[A^,5v]] 

-iT,^{Y ■ 'D)T,^ + df^^A^B.+B^A^T^^) 

+2^T^^S^aA''T^^ + AT,^T,^. (3.68) 

The so-called D-style and F-style couplings for the octet occur simply as linear com- 
binations of the most general first-order invariants of flavour SU(3) symmetry. The 
reversed sign of the kinetic term of the decuplet simply encodes the spacelike nature 
of its positive energy spinors (1/^ < 0), and the Rarita-Schwinger field propagators 
contain a polarization projector that sums over these spinors njM91all : 

n = £ ilj = {v^v^ - gn - ^S^,Sl (3.69) 
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When considering the mass renormalization of the nucleon in Chapters |4] and 

(2) 

[51 contributions from the second-order octet Lagrangian L^^J are required, which 
correspond to an NNrni vertex. This gives rise to a tadpole contribution. In full 
SU(3) form, the vertices required from L^;^J are HWLOSII : 

4d = 2aMTr[^U+]Tr[BB]+2DMTr[B{^U+,B}]+2FMTr[B[^U+,B]] , (3.70) 

where fM+ = ji^^^M + ^iW ^) is the Hermitian mass source constructed from the 
quark mass ?W [IGSS88[|WL05]I . 

Consider now the lowest-order Lagrangian for the nucleon-pion interaction by 
simplifying Equation (13.681) to involve only the nucleon doublet field ^ = (/>, n)^ 
and the SU(2) pion triplet (see Appendix IA.4I) . This is a useful approach when 
kaon loop contributions are neglected. The axial coupling constant below is simply 
defined 2iSgA=D + F llSSOSll : 

4i^=^f0-M^+^rY5X-a^S 1^. (3.71) 

The tadpole Lagrangian now takes the form: 

lI^^''"' = C2Tt[^+]^'¥, (3.72) 

where the coeffcient is a combination of the LLCs <3m. Dm and Fm, labelled C2 in 
anticipation of the analysis presented in Chapter IH 

A local, chirally symmetric form of Equation (|3.71l) can be recovered simply 
with the replacement: 

^ = (3^ + - i(Z^ -f r^)) , (3.73) 
= \{^Hd, - ir,)^ + - il,)^^) , (3.74) 

and also by replacing the product x ■ 3^,7r with a more general object: the Hermitian 

, , o o 

axial combination = — ir^)t, — ^(3^ — ?7^)^' }. The values of Mn^ Sa and 
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/;j are taken to be the nucleon mass, the axial coupling strength and the pion decay 
constant, respectively, in the chiral limit. The Goldberger-Treiman Relation relates 
the nucleon-pion interaction strength to the axial coupling gA [|GT581, and can be 
obtained by comparing the matrix elements {p\n(x)\n) and {p\dijA^{x)\n) using the 
relation between the pion field and the axial current in Equation (13.231) as per PCAC 
[ICol85]l : 

gnNN^gA-^- (3.75) 

This equation becomes exact in the chiral limit gAijn^ — 0) = gA- 



3.2.3 Electromagnetic Contributions 

The baryon form factors comprise a parameterization for the matrix element ob- 
tained from the isovector quark current = \j/Q^y^\|/, where Q, is the SU(3) quark 
charge matrix Ci = diag(2/3, — 1/3, — 1/3). To evaluate this matrix element, one 
must calculate the fully-amputated vertex for a baryon-photon interaction, wedged 
between the usual in- and out-going fermion spinors u'\p) and u^' [p'): 

(5(/)|7,|5(p)) = /(/)|Y,Fi(e2) + !^i72(e2)|^^(p), (3.76) 

for the tensor quantity o^v = H'iV'Yv}- is a positive momentum transfer = 
— (p' — p)^, and Fi and F2 are called the Dirac and Pauli form factors, respectively. 
The Sachs electromagnetic form factors Ge,m are the linear combinations: 

GEiQ^) = Fi{Q^) - ^HQ") , (3-V7) 

B 

GM{Q^)=Fi{Q^)+F2iQ^). (3.78) 



Thus, in the non-relativistic, heavy-baryon formulation: 

{B{p')\J,\B{p)) = u^\p') L,Ge{Q^) + '^^^P"^^^-^^' Gm{Q^) \u^{p). (3.79) 
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By considering the behaviour of the Sachs form factors at zero momentum transfer, 
one can construct moments and charge radii. Two such important examples that 
will be considered are the magnetic moment, and the electric charge radius of the 
isovector nucleon. Recall from Section 12.1.31 that the isovector nucleon is simply 
the combination p — n, which transforms as a vector in isospin space, chosen so that 
diagrams containing indirect couplings will cancel, and the computation will be less 
intensive. The magnetic moment //^ is simply the value of G^^ at = 0: 

p:. = GIi(Q^ = 0) (3.80) 
= !+<. (3.81) 

The first term is simply the value of the Dirac form factor of the proton at = 0, 
and the second term kJ, is the anomalous magnetic moment originating from the 
finite-size behaviour of the hadron interactions of the effective quantum field theory: 
the hadron cloud, which surrounds the nucleon. 

The electric charge radius is obtained by taking a derivative with respect to 
in the limit that equals zero: 



r2)^= lim -e —^Xi ' ■ (3-82) 



,,,,, ^^^{0 " 



For octet baryons, the magnetic moments obey the Coleman- Glashow SU(3) 
relations, related to the following Lagrangian of two independent terms IIJLMS93[ 
IWLTY071 IWLTY09al IWLTY09bll : 

i^oct" = ^(/^oTr5vO^^{F+,5v} +A/FTr5vO^^[F+,5v]) . (3.83) 

For an electromagnetic gauge field ^2^^ with field strength tensor F^v = the 
quantity has been chosen such that it is invariant under local chiral symmetry 
transformations: 

F+ = i(^V^va^ + ^i>va^^). (3.84) 
In the case of decuplet baryons, there is only a single invariant term that can be 



Chapter 3. Chiral Effective Field Theory 



Hall 38 



obtained from the group product 10 (g) 10 ® 8 that is proportional to their electric 
charge tensor qij^ [IJLMS93II : 

''de"" = ''^-A/c q.jk Tl^kl . (3.85) 
The transition Lagrangian can be written out likewise: 

i^tranT = ^^^tF,, (e,,, 5^ 5^ T.^'^- + e'^' oj fi;,,^ 5- ) . (3.86) 

These electromagnetic Lagrangians are obtained simply by collecting the photon- 
baryon terms from the electromagnetic covariant derivative. This new covariant 
derivative can be expressed by updating Equation (13.661) so that the electromagnetic 
field is included in both the vector connection Vu from Equation (13.621) and the axial 
combination from Equation (|3.64l) : 

V,-^V, + ^ie^.i^'^d^ + ^d^^) , (3.87) 
A, ^A,-^- es^.i^d^^ - d^) . (3.88) 

The covariant derivative for the pseudo-Goldstone Lagrangian can be updated in a 
similar fashion: 

Vf,U ^Vf,U + ies^^[d,U]. (3.89) 

3.3 Regularization and Renormalization 
3.3.1 Historical Overview 

The calculation and interpretation of amplitudes from a quantum field theory proved 
more subtle than other theories due to their divergent behaviour. Despite success 
in predicting hitherto unexplained phenomena, many quantities calculated using the 
relevant quantum field theory become infinite, though the known experimental value 
is finite. Consideration of the Lamb Shift in the electron energy levels in hydrogen 
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atoms (1947) prompted the first real insight into this problem. It was conceived that 
if a quantity were altered infinitely by quantum corrections so that the final result was 
finite, the initial 'bare' quantity should never have been expected to be finite. That is, 
the bare quantity becomes renormalized. For example, the bare core of an electron 
has certain properties, such as electric charge, which become altered by an infinite 
amount due to vacuum polarizations. This polarization cloud surrounding the un- 
physical, bare electron core contains all possible diagrams of electron-positron pair- 
production from virtual (off-shell) photons, which serve to screen the electron core's 
infinite charge, so that the observed, long-range charge is —1.6 x 10^^^ Coulomb, 
or —e (in units of the charge of the proton). This 'running' of the electron's charge 
to large values under deep probing from hard momenta in Bhabha scattering was 



confirmed in 1997 by the TOPAZ Collaboration at TRISTAN [|A+05l . The virtual 
particles of a quantum field theory are simply consequences of the Green's Functions 
of the equations of motion. The Fourier transform of a particle propagator integrates 
the whole momentum spectrum, with a pole on the mass shell ]^ = nP' (up to factors 
of c and h). Heisenberg's Uncertainty Principle for energy and time, (A£')(A?) > h, 
allows the extra energy of pair-production, and other processes, for sufficiently small 
time. As a corollary, the virtual interactions take place over a spacelike time interval. 



3.3.2 The Power- Counting Regime 

The Lagrangians of %PT are constructed with the intention that they can be expanded 
in a series of some expansion scale, such as small momenta or masses. Although, 
ideally, the series is convergent for a sufficiently small expansion scale, it need not 
necessarily be convergent, and instead will often take the form of an asymptotic (or 
Poincare) series. Nevertheless, in a realistic calculation, which involves calculating 
the expansion series only up to some finite order, it is desirable to be able to ensure 
that the uncertainty in the truncation is small. Thus, a knowledge of the applicable 
region of the expansion is as crucial as knowledge of the terms of the expansion 
series themselves. The range of values of the expansion scale for which a chiral 
expansion is convergent is known as the power-counting regime (PCR), and the ex- 
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pansion series is generally known as the chiral expansion. 

The PCR is the region where the quark masses are small, and higher-order terms 
in the chiral expansion are negligible beyond the order calculated. Within the PCR, 
the truncation of the chiral expansion is reliable to some prescribed precision. A 
chief focus of this thesis is to establish a formal approach to determining the PCR of 
a truncated chiral expansion quantitatively. The chiral expansion will be examined, 
and the individual low-energy coefficients of the chiral expansion will be analyzed. 
The approach involves the examination of these low-energy coefficients as they un- 
dergo the process of renormalization. This approach provides a determination of the 
PCR for a truncated expansion in %EFT. 

First, it is essential to discuss methods of regularization in the chiral loop inte- 
grals, so that the renormalization can take place. In order to renormalize a quantity, 
one must find a way to make the divergent amplitudes tractable, using a process 
called regularization. This involves solving an integral over propagators in such 
a way as to isolate the divergent piece, ready for handling with a suitable renor- 
malization scheme. There is a wide variety of regularization schemes available. 
Pauli-Villars regularization (1949) involves the introduction of fictitious, 'auxiliary' 
particles, associated with some mass scale, into a Lagrangian with a quadratic in- 
teraction. The extra formal terms in the Lagragian vanish as the mass scale is taken 
to infinity, and then a subtraction can take place. However, because Pauli-Villars is 
not a gauge-covariant scheme, it is not applicable directly to Yang-Mills theory. In 
Slavnov's regularization scheme (1971) of higher covariant derivatives, once again, 
additional terms are added to the Lagrangian, but these do not render all amplitudes 
finite, thereby requiring a Pauli-Villars or other scheme to be used for divergent 
fermion-loop Feynman diagrams. In this thesis, a finite-range regularization scheme 
is used, which has powerful benefits in establishing the PCR, as will become appar- 
ent in Chapter m 
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3.3.3 Dimensional Regularization 

Dimensional regularization (DR) (1972) is an important procedure whereby loop in- 
tegrals are analytically continued to generalized fractional dimensions and shown to 
converge [|tHV72L The infinitesimal four- volume box d^k is replaced with d^^^k, 
and the limit as E — )■ 0+ is then taken. For example, the integral over a single (Eu- 
clidean) pion propagator is easy solved in spherical polar coordinates, evaluating the 
angular part explicitl>a: 



^ lim, / (3.90) 



{2n)H^ + ml E^Q+Jo (27r)4-e/t2 + m2r(2-£/2) 

Thus the minimal subtraction scheme result is recovered correctly. 

Since there is no explicit scale-dependence in the interaction, this minimal sub- 
traction scheme makes DR suitable for use with elementary fields, where the absence 
of new degrees of freedom at higher energies is assumed. This is a powerful tech- 
nique by which the divergent term(s) of a loop integral can be obtained, and then 
handled using a renormalization scheme. 

Nevertheless, in the case of effective field theories, there exists an energy scale 
beyond which the effective fields are no longer the relevant degrees of freedom, and 
so DR is not ideally suited. Selecting a hard energy scale in the renormalization 
group equation, changes the relevant degrees of freedom in the Lagrangian. At high 
energy scales, the high de Broglie frequency would resolve the internal structure 
of the hadrons, which would be the quarks (and beyond, if such higher degrees of 
freedom exist). However, quarks are integrated out of the low-energy xEFT La- 
grangian by construction. When one calculates quantum amplitudes over this high 
energy domain, there is no guarantee that one can efficiently subtract the model- 
dependent, ultraviolet physics with a finite number of counter- terms, as is required 
for successful renormalization, unless the perturbative expansion is convergent. In- 
deed this problem of beginning with rapidly varying loop contributions, which must 
then be removed with a finite number of counter-terms, can easily be overcome. The 



^The Gamma function on C is defined as r(z) — J^dse "s^ 
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hard momentum contributions to the chiral loop integrals can be suppressed via the 
introduction of a finite-range regulator. 

3.3.4 Finite-Range Regularization 

One alternative to DR is finite-range regularization (FRR), in which one introduces 
a functional form u{k\A), known as a finite-range regulator, which controls the di- 
vergent integral at high momentum values. In this case, the integral over a single 
pion propagator would be modified as follows: 



FRR involves the choice of a finite-valued momentum cutoff A. Allowing hard, 
internal momenta to flow through a loop integral yields unphysical results, in the 
form of a divergence. The high de Broglie frequency would resolve the internal 
structure of the hadrons, which would be the quarks (and beyond, if such higher 
degrees of freedom exist). Therefore, a finite value of A is suitable for an effective 
field theory, where quarks are integrated out of the Lagrangian by construction. The 
choice of parameter A determines how fast the integral will now converge, and the 
regulator function should satisfy u\k=Q = 1 and u\k^oo = 0. The exact functional 
form chosen for the regulator should be independent of the result of calculation, as 
long as the perturbative expansion is convergent, that is, one works within the PGR. 

FRR has already been shown to be a powerful technique in solving the chiral 
extrapolation problem and identifying the PGR. The infinite series is resummed so 
that leading-order terms are large and the series converges. A variety of choices of 
functional forms for the regulator have been demonstrated to agree with each other, 
and with DR, in extrapolating lattice QGD results for the mass of the nucleon to 
physical quark masses ULTYOSI . Thus, the results of calculations using FRR are 
consistent with DR within the PGR. 

Gonsider the example of a one-pion loop contribution for a nucleon, denoted E/^, 
with constant coefficient Xn- (This type of calculation is considered in more detail 




(3.91) 
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in Chapter HI) The chiral expansion for the mass of the nucleon in this simple case, 
with one pion loop only, takes the form: 

Mn = {ao + aiml + a4m^ + O (m^) } + , (3.92) 

working to chiral order 0{m^). The chiral expansion comprises a polynomial ex- 
pansion in and the contribution from the one-pion loop. Each of the coefficients 
ao, a2 and a4 is renormalized by the contributions from the loop integral, at each or- 
der. The result of the integral using DR is equivalent to a massless renormalization 
scheme, with no explicit momentum cutoff: 



A 



dk-^ ^ (3.93) 



n Jo k^ + i ^ 

= fd^'^'+^-'f -"'')+'"■ (3.94) 

n Jo k^+ 

In a massless renormalization scheme, there is no explicit momentum cutoff, so each 
of the coefficients a, undergoes an infinite renormalization or none at all: 



2Xn 



CQ = ao + -^ I dkk^, (3.96) 
7t Jo 

C2 = a2-^ f dk, (3.97) 



Tl Jo 

C4 = a4 + 0, etc. (3.98) 

By contrast, in a ERR scheme, a momentum cutoff A is introduced, and the chiral 
expansion is resummed. Using a sharp momentum cutoff A: 

E;v(A) = % t^kj^ (3.99) 

r A i\ 

(3.100) 



TT Jo k'^ + i 

2%yv ^a3 .2^3 , r A 
Amjj + arctan 



m 



- ^Am^ + ^,rnl - ^^\ml + (m°) . (3.101) 
TT 3 TT TT A 



Chapter 3. Chiral Effective Field Theory 



Hall 44 



The result obtained from DR can be recovered in an FRR scheme by taking the 
regularization scale parameter A to infinity: 

co = «o + ^A^ (3.102) 

C2 = a2-^A, (3.103) 
n 

C4 = ^4 — , etc. (3.104) 

n A 

Thus, DR applied outside the PGR could be considered equivalent to a model with an 
arguably injudicious choice of cutoff scheme. The polynomial expansion of hadron 
mass is not expected to converge, and indeed it does not, using DR xPT, as men- 
tioned by Young, et al. [IYLT03II . Outside PGR, the expansion breaks down since 
the chiral expansion is truncated without an attempt to estimate the higher-order 
contributions [|LTY05llLTY06ll . 

In addition, because FRR involves the resummation of the higher-order terms 
of the chiral expansion, it affords an opportunity to perform a calculation beyond 
the PGR. Using FRR, one must select a value for the ultraviolet regularization scale 
A. The choice in the value of A is irrelevant within the PGR, where the results of 
extrapolations are scheme-independent (so long as A is not chosen to be too small, as 
explained in Section l4.2.1l) . Nevertheless, the principal exercise of this thesis will be 
to handle any scheme-dependence occuring in a %EFT calculation outside the PGR. 
By quantifying the scheme-dependence one arrives at a rigorous procedure for using 
FRR beyond the PGR. 

In Ghapter HI a variety of finite-range regulators are used and compared. For 
example, the Heaviside Step Function u^{k;A) = 0(A — fc) is an acceptable choice; 
however, it is unfavorable for finite-volume considerations because discrete lattice 
momenta are either fully included in the integral or not included at all. This results 
in inconvenient finite-volume artefacts. In the investigation of the nucleon mass, the 
family of smoothly attenuating dipole regulators will be considered. The general 
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Figure 3.1: Behaviour of three dipole-like regulators as a function of momentum k, for a regulator 
parameter A = 1 .0 GeV. 

multiple-dipole function of order n takes the following form, for a cutoff scale of A: 

u„{k;A)=\l + ^j . (3.105) 

The standard dipole is recovered for n = I. The cases n = 2,3 are the 'double-' 
and 'triple-dipole' regulators, respectively. The behaviour of each of these three 
attenuators is shown in Figure 13. 1[ These functional forms allow one to interpolate 
between the dipole regulator and the step function (which corresponds to n — )■ oo). 
These n— tuple-dipole regulators generate extra non-analytic terms. 

It has been suggested in the literature that the only FRR scheme consistent 
with chiral symmetry uses the step function regulator [BH M04II . Djukanovic et 
al. [IDSGS05II have demonstrated that more general functional forms can be gen- 
erated by proposing a scheme in which the regulator function is interpreted as a 
modification to the propagators of the theory, obtained from a new chiral symmetry- 
preserving Lagrangian. Higher derivative coupling terms are built into the Lagrangian 
to produce a regulator from the Feynman Rules in a symmetry-preserving manner. 
Alternatively, one can choose the regulator judiciously such that any extra scheme- 
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dependent non-analytic terms are removed to any chosen order. The regulators used 
in the present investigation follow the latter approach. For example, the n— tuple 
dipole regulators generate extra non-analytic terms in the chiral expansion of Equa- 
tion (|4.2I) in Chapter HI at higher chiral orders. An explicit example of this is shown 
for the quenched p meson mass in Section |6.1.2[ once the renormalization scheme 
has been introduced. 



Chapter 4 

The Intrinsic Scale of the Nucleon 



"There lies the originality of our approach: to deduce common sense from the quan- 
tum premises, including its limits— that is, to demonstrate also under which condi- 
tions common sense is valid, and what is its margin for error. . . 

[W]e no longer explain reality from our mental representation of it, taken for 
granted without question: but it is this representation, . . . that we want to explain/J " 
(Omnes, R. 2002. Quantum Philosophy: Understanding and Interpreting Contem- 
porary Science p. 165) [|Omn02ll 

4.1 Renormalization Issues for the Nucleon Mass 

In chiral effective field theory (xEFT), the nucleon mass may be written as an or- 
dered, chiral expansion in the quark mass. The Gell-Mann— Oakes— Renner Rela- 
tion from Equation (13.431) entails the proportionality m^. By considering the 

o 

renormalization of the nucleon mass Mn^ Mn from the Lagrangian in Equation 
(13.711) . the chiral expansion will generally include a polynomial in m| and non- 
analytic terms obtained from the chiral loop integrals. In addition, to establish a 
model-independent framework in %PT, the expansion must display the properties of 
a convergent series for the terms considered. Recall that within the power-counting 
regime (PGR) the higher-order terms of the expansion may be regarded as suffi- 
ciently small for the truncation of the chiral expansion to be reliable to a prescribed 
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precision. However, truncated expansions are typically applied to a wide range of 
quark (or pion) masses, with little regard to a rigorous determination of the PGR. In 
the case of the nucleon mass, evidence suggests that the PGR is small: limited to 

< 200 MeV at 1% accuracy at the chiral order (^(m^logmjt) nBea04all IILTY05L 
This estimate of the PGR of %PT was identified by comparing the results of infrared 
regularization, dimensional regularization (DR) and a variety of finite-range regula- 
tors in analyzing lattice quantum chromodynamics (lattice QGD) simulation results. 
The different regularization schemes constitute different ways of summing higher- 
order terms in the chiral expansion. Thus, the PGR is manifest when the pion mass 
dependence of the nucleon mass is independent of the renormalization scheme. In 
addition, the asymptotic nature of the chiral expansion places the focus on the first 
few terms of the expansion. 

A survey of the literature for the baryon sector of %EFT illustrates the rarity 
of calculations beyond one-loop IIMB991 |MB06[ ISDGSOTl . and there are currently 
no two-loop calculations that incorporate the effects of placing a baryon in a finite 
volume. With only a few terms of the expansion known for certain, knowledge of the 
PGR of %EFT is as important as knowledge of the expansion itself. Though scheme- 
dependent, it is worthwhile to note that, using a dipole regulator with A = 0.8 GeV, 
the coefficient of the induced m| term compares favorably with the infinite-volume 
two-loop calculation IIMB991 ILTY041 ILTYOSi IMB061 ISDGS07H . 

4.1.1 Chiral Expansion of the Nucleon Mass 

The nucleon mass expansion formula can be expressed in a form that collects the 
non-analytic behaviour into the loop integral contributions: 



= co + C2ml + XNtnl + CAtnl + ( - 7— tXa + Xt ) m^log — + {ml). (4.2) 



(4.1) 
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Figure 4.1: The pion loop contribution to the self-energy of the nucleon, providing the leading 
non-analytic contribution to the nucleon mass. All charge conserving transitions are implicit. 
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Figure 4.2: The pion loop contribution to the self-energy of the nucleon, allowing a transition to a 
nearby and strongly-coupled decuplet baryon. 

f \ 
\ / 

^ 

N N 
Figure 4.3: The O [niq) tadpole contribution to the nucleon self energy. 

The superscript A denotes the scale-dependence of the af coefficients. The analytic 
terms in can be written as a polynomial with renormalized coefficients c,. The 
non-analytic contributions arise from the self energy integrals (E), which correspond 
to the diagrams in Figures HTTI through 14.31 

It is essential to note that the degrees of freedom present in the residual se- 
ries coefficients af are sufficient to eliminate any dependence on the regulariza- 
tion scale parameter A, to the order of the chiral expansion calculated: in this case 
0(m^logm;t). Any differences observed in results obtained at the same chiral order, 
but with different regularization schemes, are a direct result of considering data that 
lie outside the PGR (provided that the scale A is not chosen so small that it introduces 
an unphysical low-energy scale). 
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4.1.2 Chiral Loop Integrals 

Each of the loop integral contributions to the nucleon mass can be simplified to a 
convenient form by taking the non-relativistic heavy-baryon limit, and performing 
the pole integration for ko. The integrals may be expanded out to a particular chiral 
order, in this case order o{m^\ogmTi), to obtain an analytic polynomial with coef- 
ficients i»f , and the leading-order non-analytic term. Using a finite-range regulator 
u{k;A): 

= b^'"^ + b^/ml + + Z^J'^m^ + o [ml) , (4.4) 



EA(m^;A) = ^ 



kV{k;A) 
(0{k) {A + (o{k)) 



+ b2'^ml + b 



A,A 4 



4^%Am^ log — + 0(m^) 



(4.5) 
(4.6) 



= C2»4('i>J'+(.fmi + x,"4logy + 0(m5)V (4.8) 

where p is an implicit mass scale from the logarithm, (o{k) = ^/k'^ + and A is 
the nucleon-delta baryon mass-splitting, treated as a perturbation in the approximate 
flavour symmetry. The mass of the A baryon is chosen to be the centre of its Breit- 
Wigner resonance. 

The bf coefficients renormalize the residual coefficients of the chiral expansion 
of Equation (14.11) . to obtain the scale-independent coefficients c,. Though both the 
af^ coefficients and the bf^ coefficients are scale-dependent, adding them together 
at each order results in a scale-independent coefficient. These are the renormalized 
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coefficients c,-. Explicitly: 

Co 

C2 
C4 



A I 7 A,A' , , A,A 

A I 7 A.A' I , A.A I , A.t' 

= + b^'^ + b^'^ + b^'' , etc. 



(4.9) 
(4.10) 
(4.11) 



This is simply a slight generalization from the worked example in Section 13.3.41 
Dimensional analysis reveals that the coefficients b^ are proportional to A^-^^'^ . Thus 
it can be realized that as the cutoff scale A tends to infinity, the result from DR, as 
described in Section 13.3.31 is recovered. At any finite A, a partial resummation of 
higher-order terms is introduced. Previous studies indicate that extrapolation results 
show very little sensitivity to the precise functional form of the regulator IILTY04L 

A modification is now made to the integrals of Equations (14.31) through (14.81) . 
by subtracting terms from their Taylor expansion, thus absorbing them into the 
corresponding low-energy coefficients c,. This achieves the renormalization to a 
chosen chiral order. In this case, only the low-energy coefficients cq and C2 will be 
analyzed. The amplitudes for each process are thus altered: 

E^(m^;A) = — j d k^^^—^ -b^^ -b^ (4.12) 
= XNml + bf^^ml + 0{ml), (4.13) 



2%^] &{k){A + &{k)) 



J A.A J A,A 2 



XaK \og — + 0{ml), 



AnA 



(4.14) 
(4.15) 
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E,arfK;A) = ^2m^[j^j^^—^^ ^2 J (4-16) 

= C2ml(^^^'ml + 'itml\og'^ + 0{ml)^ (4.17) 
= C2mlotad{ml\A) . (4.18) 



Note that the coefficient of the tadpole amplitude contains the renormalized low- 
energy coefficient C2- This is because the same coefficient C2 from the chiral expan- 
sion occurs in the tadpole Lagrangian in Equation (13.721) . The tilde (~) denotes that 
the integrals are written out in renormalized form to chiral order o(m^). As the bf 
coefficients are regulator and scale-dependent, the subtraction reshuffles this depen- 
dence into higher-order terms. The coefficients Uq and of the analytic terms in 
the chiral expansion in Equation (14.21) automatically become the scale-independent 
renormalized coefficients cq and C2- 

With the renormalized integrals specified, the finite-range regularization (ERR) 
modified version of the chiral expansion in Equation (|4.2I) takes the form: 

= CQ + C2ml{l + dtad{ml,A))+a^ml + tN{ml,A) +EA(m2,A) . (4.19) 

The ^4 term is left in unrenormalized form for simplicity. Indeed, the coefficient b^ 
can be evaluated by expanding out corresponding loop integrals, such as in Refer- 
ence I1YLT03L However, the focus here is on the behaviour of cq and C2. 



4.1.3 The Sigma Term 

In addition to the mass of the nucleon in the chiral limit cq, the low-energy con- 
stant (LEC) C2 corresponding to the tadpole vertex is of interest phenomenologically 
because, by inspection of Equation (13.721) . it is a measure of the explicit chiral sym- 
metry breaking of the relevant flavour symmetry group. That is, a sigma term can 
be defined for the light quarks up and down, and the explicit breaking of the group 
SU(2)v (g) SU(2)^ may be investigated. In order to obtain a value for the sigma term 
relating to the heavier strange quark, %EET has been used to study the explicit break- 



Chapter 4. The Intrinsic Scale of the Nucleon 



Hall 53 



ing of the baryon octet representation of SU(3) IIGasSli INK871 IBM97[ lYTlObll . 

It can be seen that the term C2m^, and higher-order terms in the nucleon mass 
expansion formula of Equation (14.21) . will disappear if the chiral symmetry breaking 
quark mass is zero. To investigate this, one can consider how the QCD Hamilto- 
nian behaves under commutation with the three-component axial charge operator of 
flavour SU(2). If the symmetry were unbroken, all quantities 2^|0) would vanish, 
so the commutator [Q\tO-CqcT)\ would also vanish. In the more general case, con- 
sider two applications of the commutator, which yield the symmetry breaking mass 
term from the total Hamiltonian. This defines the pion-nucleon sigma term 
[ICDTTI lETWOOl IWLTOOi IHMRW061 lYTlObll : 

= \m[Q\^ [ei,^QCD]]|A^) (4.20) 

= {N\{muUu + mJd)\N) = {N\^,^,\N) . (4.21) 

Under the simplification of mass degeneracy between quark fields (which is approx- 
imately true under flavour SU(2)), one can apply the Feynman— Hellmann Theorem 
[ |Fey39[ and recover the important result for small m\: 

O-siN = mq^^ = C2ml + O {rnj^) . (4.22) 
aniq 

That is, the value of the sigma term is dominated by the leading-order term with 
coefficient C2- The violation of this axial symmetry is therefore important for un- 
derstanding the behaviour of hadrons, because a non-zero sigma term affects the 
structure of the interaction between hadrons and the meson cloud which surrounds 
them, and provides a small, but not statistically insignificant contribution to the total 
mass of the hadron. 

The standard result for the sigma term using SU(2) xPT, incorporating meson 
loop corrections, is: Etia? = 35 ± 5 MeV UGasSlll . By analyzing data from np and 
rni scattering experiments [|H6h83L an early analysis by Gasser suggests a value of 
^ttA' = 45 ± 8 MeV [IGLS91II . The currently accepted value of the sigma term, due 
to the work of Koch, is larger than the theoretical value: E^yv = 64 ± 8 MeV [|KP80[ 
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IKoc82ll . A more recent analysis of the experimental data by Pavan, incorporating 
a partial wave and dispersion relation analysis, suggests an even higher value of 
^nN = 79 ± 7 MeV [IPSWA02L In contrast, calculations from two-flavour dynamical 
quark lattice QCD comparatively underestimate the value of the sigma term. In a 
study by Giisken, it was found that E^ia? = 18 ± 5 MeV, by direct calculation of the 



scalar matrix element in Equation (14.211) [|G"'"99l . This apparently low value for the 



sigma term was found to be a consequence of its sensitivity to chiral extrapolation, 
and large pion masses (above 500 MeV) were used in the extrapolation ULTWOOi 
lYTlObL 

4.1.4 Scheme-Independent Coefficients 

The chiral coefficients Xn, Xa and Xt for each integral are defined in terms of the 
pion decay constant, which is taken to be fn = 92.4 MeV, and the axial coupling 
parameters D, F and C which couple the baryons to the pion field, as shown in 
the Lagrangian -^oct&dec Equation (13.681) . The coeffcient C2, which occurs in the 
tadpole loop integral of Equations (14.161) through (14.181) . is a combination of the 
LECs Om, Dm and Fm, which occur in the tadpole Lagrangian of Equation (13.701) . 
Though C2 is treated as a fit parameter, the phenomenological values for the D, F and 
C couplings are used, applying the SU(6) flavour-symmetry relations IIJen92[|Leb95ll 
to yield C = -2D, F ^ ^D and the value D = 0.76: 

These coefficients are constant and remain unaffected by renormalization scale or 
finite- volume effects. Ultimately, one may try to determine these directly from lattice 
simulation results. Nevertheless, because of the limited number of lattice simulation 
results currently available, this analysis will focus on the determination of cq, C2 and 
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the nucleon mass M^y. 



4.1.5 Finite- Volume Effects 

In lattice QCD, the introduction of finite- volume effects become significant for small 
box sizes. The expansion parameter 1 /L contributing to finite- volume effects should 
be of the same order of magnitude as the momenta for the perturbation scheme 
to remain valid. If L is small, the exponential factor e^"^^^ no longer suppresses 
the finite- volume corrections HBNSIOL As a general rule, the characteristic di- 
mensionless quantity mjiL specifies the £-regime through the condition m-nL < 1 
IIHan90[ IHL91[ IHL90L This is a breakdown region in xPT, since divergences in 
the leading-order pion contributions cannot be approximated by standard perturba- 
tive techniques HENS 1 OH . 

Since the results of lattice simulations reflect the presence of discrete momentum 
values associated with the finite volume of the lattices, the formalism of XEFT must 
also take into account these finite-volume effects. xEFT is ideally suited for exam- 
ining finite-volume effects, because of its accurate characterization of the dominant 
infrared physics. In order to accommodate the effect of the finite volume, the con- 
tinuous loop integrals occurring in the meson loop calculations in an infinite volume 
are transformed into a sum over discrete momentum values. The difference between 
a loop sum and its corresponding loop integral is defined to be the finite- volume cor- 
rection, which should vanish for all integrals as m^L becomes large [|GL88[|Bea04bll . 
While Equation (14.191) is useful in describing the pion mass evolution of the nucleon 
mass, for the consideration of lattice QCD results, one also needs to incorporate 
corrections to allow for the finite- volume nature of the numerical simulations. As 
the pion is the lightest degree of freedom in the system, it is the leading-order pion 
loop effects that are most sensitive to the periodic boundary conditions. The correc- 
tions can be determined by considering the transformation of each loop integral in 
Equations (14.121) . (14.141) and (14.161) into a discrete sum for a given lattice size. The 
three-dimensional integrals can be replaced by summations over all possible mo- 
mentum values llAAL"'"06l . It is useful to define the finite- volume correction to the 
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loop integral, by convention, by subtracting the integral from the sum quantity. This 
technique will be used to correct for finite- volume effects encountered in Chapters H] 
through |7J 

The finite- volume correction 6^^*-' can be written as the difference between the 
finite sum and the integral: 



(4.26) 



where / = A^, A, and the integrands are denoted Ii(k,m^,A). The finite- volume cor- 
rections to the tadpole contribution are not considered in this investigation because 
of subtleties in their behaviour at large m-n. Details regarding the finite- volume be- 



haviour of the tadpole amplitude are discussed in Appendix IB. 3 1 and a more general 



discussion of its convergence properties occurs in Section 18.21 By adding the rel- 
evant finite-volume correction to each loop contribution, the finite- volume nucleon 
mass can be parameterized: 

M^^co + C2m\{\ + Otad) + 4m\ + (E^ + 6^^^) + (Ea + 5^^) . (4.27) 

It is also shown that the finite- volume corrections are independent of the regular- 
ization scale A in this domain. In Figures |44] and 1431 the scale-dependence of the 
finite-volume corrections is shown for a dipole regulator (from Equation (|3.105|) in 
Chapter[3]) and a 2.9 fm box (the same box size used for the PACS-CS data [|A+09|| ). 
It is of note that choosing A too small suppresses the very infrared physics that one 
is trying to describe. Thus, caution should be exercised in choosing a suitable value 



of A. Figures 14.51 and 14.71 show the behaviour of the finite- volume correction for a 
4.0 fm box, and the corrections are much smaller, as expected. 

For large A, the finite- volume corrections, displayed in Figures 1441 through I4.7[ 
saturate to a fixed value. Provided that A > 0.8 GeV, the estimated finite- volume cor- 
rections are stable for light pion masses. In order to preserve the scale-independence 
of the finite-volume corrections, their asymptotic result will be used. This approach 



has been demonstrated to be successful in previous studies [lAK+041 . Numeri 
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Figure 4.4: Behaviour of the finite- 
volume corrections 5^^*^ vs. A on a 
2.9 fm box using a dipole regulator. 
Results for two different values of 
are shown. 



Figure 4.6: Behaviour of finite- 
volume corrections 5^^*-" vs. A on a 
2.9 fm box using a dipole regulator. 
Results for two different values of 
are shown. 
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Figure 4.5: Behaviour of finite- 
volume corrections 5^^^ vs. A on a 
4.0 fm box using a dipole regulator. 
Results for two different values of 
are shown. 



Figure 4.7: Behaviour of finite- 
volume corrections 5^^*^ vs. A on a 
4.0 fm box using a dipole regulator. 
Results for two different values of 
are shown. 



cally, this is achieved by evaluating the finite-volume corrections with a parameter, 
A' = 2.0 GeV, 5f^^ = 6f^'^(A') . It should be noted that this is equivalent to the more 
algebraic approach outlined by Beane IIBea04bL 



4.2 The Intrinsic Scale: An Example by Construction 



This xEFT extrapolation scheme to order o(m^ log m^t) will be used in conjunc- 
tion with lattice QCD data from JLQCD [IO+08l . PACS-CS [|A+09l and CP-PACS 
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IIAK+021 to predict the nucleon mass for any value of m^. The full set of data from 
each of these collaborations is listed in Appendix O Tables IC.ll through IC.3I The 
JLQCD data were generated using overlap fermions in two-flavor QCD, but the lat- 
tice box size for each data point is ~ 1.9 fm, smaller than the other two data sets. 
The PACS-CS data were generated using non-perturbatively (a) -improved Wilson 
quark action at a lattice box size of ~ 2.9 fm, but the data set only contains five data 
points and a large statistical error in the smallest point. The CP-PACS data were 
generated using a mean field improved clover quark action on lattice box sizes for 
each data point varying from ~ 2.2 fm to ~ 2.8 fm. 

The lattice data used in this analysis will be used to extrapolate to the phys- 
ical point by taking into account the relevant curvature from the loop integrals in 
Equations (14.131) . (14.151) and (14.171) . As an example, a regularization scale of A = 1.0 
GeV was chosen for Figures 14.81 through I4.10[ where the finite- volume corrected 
effective field theory appears concordant with previous QCDSF-UKQCD results 



||AK"'"041 . An extrapolation or interpolation is achieved by subtracting the finite- 
volume loop integral contributions from each data point and then fitting the result to 
obtain the coefficients cq, C2 and using Equation (14.191) . The finite- or infinite- 
volume loop integrals are then added back at any desired value of m^. 

If the regularization scale is altered from the choice A = 1.0 GeV, the extrapo- 
lation curve also changes. This signifies a scheme-dependence in the result due to 
using lattice QCD data beyond the PCR. To demonstrate this, consider the infinite- 
volume extrapolation of the CP-PACS data, as shown in Figure 14.1 1[ Figure 14.1 II 
also shows that the curves overlap exactly when is large, where the lattice data 
reside, and they diverge as the chiral regime is approached. 

Consider an insightful scenario, whereby a set of ideal 'pseudodata' with known 
low-energy coefficients is produced, using the formula from Equation (14.271) . A 
particular regularization scale is selected and a dense and precise pseudodata set is 
generated, which smoothly connects with the lattice simulation results. In this case, 
the pseudodata are converted to infinite-volume results in order to ensure that the 
following analysis is not simply a consequence of finite- volume effects. If all the 
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Figure 4.8: Example dipole extrapolation based on JLQCD data Il0+08l . box size: 1 .9 fm. 
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Figure 4.9: Example dipole extrapolation based on PACS-CS data l|A+09l . box size: 2.9 fm. 
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Figure 4. 10: Example dipole extrapolation based on CP-PACS data ||AK+02| . lattice sizes: 2.3 - 
2.8 fm. 
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data considered lie within the PGR then the choice of regularization scale is irrele- 
vant, and the finite-range regularized chiral expansion is mathematically equivalent 
to scale-invariant renormalization schemes, including DR. This scenario will form 
the basis of the investigation of the PGR, and ultimately, will lead to determining the 
existence of an intrinsic scale hidden within the lattice QGD simulation results. 

The pseudodata are produced by performing an extrapolation such as shown in 
Figures 148) through 14. 101 The difference is that 100 infinite- volume extrapolation 
points are produced close to the chiral regime. The exercise is to treat these pseu- 
dodata as if they were lattice QGD data. Glearly, a regularization scheme must be 
chosen in generating the pseudodata. In this case, a dipole regulator was chosen and 
pseudodata were created at Ac = 1 .0 GeV. 

The regularization-dependence of the extrapolation is characterized by the scale- 
dependence of the coefficients c,-. These coefficients are obtained from fitting the 
pseudodata. Gonsider how cq and C2 behave when analyzed with a variety of regu- 
larization scales in Figures l4.12l and l4.13l By using infinite- volume pseudodata, one 
eliminates the concern that the variation in c/ with respect to A is merely a finite- 
volume artefact. 

Three pseudodata sets are compared, each with different upper bounds on the 
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Figure 4.12: Behaviour of cq vs. regularization scale A, based on infinite-volume pseudodata 
created with a dipole regulator at Ac = 1.0 GeV (based on lightest four data points from PACS-CS). 
Each curve uses pseudodata with a different upper value of pion mass 
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Figure 4.13: Behaviour of C2 vs. A, based on infinite-volume pseudodata created with a dipole 
regulator at Ac = 1 .0 GeV (based on lightest four data points from PACS-CS). 
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range of considered in the fit. An increasing regulator-dependence in cq and 
C2 is seen as the data extend outside the PGR. In Figures 14.121 and 14.131 the be- 
haviour of the fit parameters cq and C2, respectively, are shown as functions of the 
regularization scale A for different values of j^^^. A steep line indicates a strong 
scheme-dependence in the result, and this occurs for data samples extending far out- 
side the PGR. Scheme-independence will appear as a horizontal line, as is apparent 
for jjjax < 0-04 GeV^, in Figures |4. 1 2l and l4~T3l This indicates that the pseudodata 
lie within the PGR. 

Note that in both figures all three curves (corresponding to different values of 
'^n,max) arrive at stable values for cq and C2 on the right-hand side of the plot, corre- 
sponding to large A. To read off the values of cq and C2 for large A is tempting, but 
this does not yield the correct values of cq and C2, which are known by construction. 
The correct values of cq and C2 are recovered at A = 1.0 GeV. 

The analysis of the pseudodata in Figures 14.121 and 14.131 shows that even as the 
value of jjjj^x is changed, the correct value of cq is recovered at exactly A = Ac, 
where the curves intersect. The same value of A for the intersection point is obtained 
by analyzing C2- This suggests that when considering lattice QGD results extending 
outside the PGR, there may be an optimal finite-range cutoff. Physically, such a 
cutoff would be associated with an intrinsic scale reflecting the finite size of the 
source of the pion dressings. Mathematically, this optimal cutoff is reflected by an 
independence of the fit parameters on ^^j^. 

By analyzing the pseudodata with a different regulator, for example, a triple- 
dipole regulator. Figures |4. 141 and |4. 151 show that the scale of the intersection is no 
longer a clear point, but a cluster centred about 0.5 to 0.6 GeV. The triple-dipole will 
of course predict a different optimal scale, since the shape of the regulator is different 
from that of the dipole used to create the pseudodata. The essential point of this 
exercise is that clustering of curve intersections identifies a preferred renormalization 
scale that allows one to recover the correct low-energy coefficients. In this case, the 
crossing of the dash and dot-dash curves (from fitting) clearly identifies Af^^^ = 0.6 
GeV as a preferred regularization scale, which reflects the intrinsic scale used to 
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param. 


input 


A scale _ 
''Hrip 
2 

fYl 

jr,max 


= 0.6 
= 0.25 


A scale _ 
2 

jr.max 


= 0.6 
= 0.5 


Atrip = 2.4 
'"Imax = 0-25 


Atrip = 2.4 


Co 


0.902 


0.901 




0.902 




0.899 


0.896 


C2 


3.00 


3.07 




3.07 




3.17 


3.23 



Table 4.1: A comparison of the parameters cq (GeV) and C2 (GeV^') at their input value (pseu- 
dodata created with a dipole at Ac = 1.0 GeV) with the values when analysed with a triple-dipole 
regulator. Different values of Atrip (GeV) and (GeV^) are chosen to demonstrate the scheme- 

dependence of Co and C2 for data extending outside the PGR. Note: the values of cq and C2 are 
calculated from an ideal model and thus they are exact; there are no statistical uncertainties. 

create the data. Table l4. 1 I compares the values for cq and C2 recovered in this analysis 
for two different regularization scales: the preferred value A^^^^ = 0.6 GeV, and a 
large value At^p = 2.4 GeV reflecting the asymptotic result recovered from DR. The 
input values of cq and C2 used to create the pseudodata are also indicated. 

4.2.1 Lower Bounds for the Regularization Scale 

Figures 14.131 and 14.151 clearly indicate that the finite-range renormalization scheme 
breaks down if the FRR scale is too small. This is because A must be large enough to 
include the chiral physics being studied. The exact value of a sensible lower bound 
in the FRR scale will depend on the functional form chosen as the regulator. 

Figure |4T3] shows that the renormalization for C2 breaks down for small values 
of A. FRR breaks down for a value of Ajip much below 0.6 GeV, simply because the 
coefficients of the loop integral expansion in Equations (14.131) . (14.151) and (|4.17l) 
are proportional to A^^^') . For higher-order terms with large /, the coefficients will 
become large when A is small. In theory, these very large terms add up to zero, 
and so the limit A — )• amounts to neglecting the infrared physics of the hadron. 
In practice, the finite curvature and higher-order terms of the residual series are not 
large enough to cancel the small-A behaviour of the bf" coefficients, which dominate. 
This adversely affects the convergence properties of the chiral expansion. On the 
other hand, one obtains a residual series expansion with good convergence properties 
when A reflects the intrinsic scale of the source of the pion dressings of the hadron 
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Figure 4.14: Behaviour of cq vs. A, based on infinite-volume pseudodata created with a dipole 
regulator at Ac = 1 .0 GeV but subsequently analyzed using a triple-dipole regulator. 
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Figure 4.15: Behaviour of C2 vs. A, based on infinite-volume pseudodata created with a dipole 
regulator at Ac = 1 .0 GeV but subsequently analyzed using a triple-dipole regulator. 
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in question. 

The pseudodata analysis provides a good indication of a lower bound for A using 
a dipole regulator: A^ip > 0.6 GeV. Similarly, Figure |4T5] suggests a lower bound 
for the triple-dipole regulator: Atrip ^0-3 GeV. The same analysis can be repeated 
for the double-dipole regulator to obtain Adoub ^ 0.4 GeV. 

One can also estimate the lowest reasonable value of A by considering arguments 
from phenomenology. Based on the physical values of the sigma commutator and 
the nucleon mass, a pion mass of m^t ^ 0.5 GeV is a suitable upper bound for the 
radius of convergence [|BHLO02[ IYHL091 lYTlOaL This follows from the estimate 
of the two-flavour pion-nucleon sigma term due to Gasser IIGLS91II . Using the Gell- 
Mann— Oakes— Renner Relation niq oc m\: 



For good convergence, it is expected that the sigma term is dominated by the leading- 
order C2 term. The second and third terms in the expansion are as large as the leading- 
order C2 term for m^^ ~ 0.5 GeV. Therefore, in order to maintain good convergence 
of the chiral expansion whilst ensuring the inclusion of important contributions to 
the chiral physics, one should choose a scale Ajharp ~ 0.5 GeV for a sharp cutoff 
(step function) regulator. To compare this estimate for the sharp cutoff to that of 
dipole-like regulators, one can calculate the regularization scale required such that 
uf^{k\A) = 1/2 when the momentum takes the energy scale of Asharp- This results in 
a rough estimate for a sensible value for the dipole, double-dipole and triple-dipole 
regulators. These values are Ajip ~ 1.1 GeV, Adoub ~ 0.76 GeV and At^p ~ 0.66 
GeV, respectively. 

In the forthcoming chapter, a range of regularization scales will be considered, 
and the intersections of the curves for the low-energy coefficients will be used to 
construct fits that include data sets that extend outside the PGR. This is done in order 
to identify the presence of an intrinsic scale for the pion source and an associated 
preferred regularization scale. 



ciml + iNtnl + CAm\ -f O (m^) ^ 45 MeV . 



(4.28) 
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Chapter 5 

Results for the Mass of the Nucleon 



"The datum is a classical property concerning only the instrument; it is the expres- 
sion of a fact. The result concerns a property of the quantum world. The datum is an 
essential intermediary for reaching a result." (Omnes, R. 2002. Quantum Philoso- 
phy: Understanding and Interpreting Contemporary Science p. 209) [ )Omn02| 

This quotation, and those introduced in Chapters [6] to [8l contain an argument that 
links data, results, theory and experience. 

The aim of this chapter is to apply an analysis that allows a reliable extrapolation 
of the nucleon mass to the physical point by obtaining an optimal regularization 
scale, using lattice quantum chromodynamics (lattice QCD) simulation results. The 
identification of an optimal regularization scale, along with its associated systematic 
uncertainty, indicates the degree to which the lattice QCD simulation results extend 
beyond the power-counting regime (PCR). This quantifies and effectively handles 
the scheme-dependence of chiral extrapolations. Ultimately, the agreement among 
optimal regularization scales obtained from different simulation results indicates the 
existence of an intrinsic scale that characterizes the interaction between the pion 
cloud and the core of the nucleon. Such an agreement will be demonstrated through 
the results in this chapter, and Chapter |7J In Chapter [6l the procedure developed 
in this thesis for analyzing the renormalization flow of the low-energy coefficients, 
obtaining a possible intrinsic scale (or a range of acceptable regularization scales), 
and performing a robust chiral extrapolation will be tested. 
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In the previous Chapter, extrapolation of the lattice results was discussed in the 
context of finite-range regularized chiral effective field theory (xEFT). The scheme- 
dependence of the various extrapolations was analyzed. A method was developed for 
extracting an optimal finite-range regularization scale from ideal pseudodata. Since 
the pseudodata were generated at a known scale Ac, they contain an intrinsic scale by 
construction, and so it was demonstrated that an optimal finite-range regularization 
scale could be extracted from the pseudodata by analyzing the scale-dependence of 
the low-energy coefficients. This optimal scale was the same value as the intrinsic 
scale built into the pseudodata. 

The pseudodata example leads the researcher to consider whether actual lattice 
QCD simulation results have an intrinsic cut-off scale embedded within them. That 
is, by analyzing lattice QCD data in the same way as the pseudodata, can a similar 
intersection point be obtained from the renormalization- scale flow of the low-energy 
coefficients? If so, it would indicate that the lattice data contain information regard- 
ing an optimal finite-range regularization scale, and thus provide evidence for the 
existence of an underlying intrinsic scale in the nucleon-pion interaction. 

5.1 Evidence for an Intrinsic Scale in the Nucleon Mass 
5.1.1 Renormalization Flow Analysis 

Consider the mass of the nucleon as extrapolated from the results of lattice QCD 
simulations. The results for cq and C2 as a function of the regularization scale A are 
now presented for lattice QCD data from the collaborations: JLQCD, PACS-CS and 
CP-PACS. Initially, the chiral expansion, calculated to chiral order o(m^), should 
be used for fitting: 

MN = co + C2ml{\+ d,ad (m2 , A) ) + (m2 , A) + Ea (m^ , A) . (5.1) 

Thus, the relevant fit parameters used in the extrapolation are cq and C2 only. Results 
for the higher chiral order of (m^logm^i) will be discussed in Section [5. 1.3[ The 
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resultant renormalization flows, using a dipole regulator, are shown in Figures 15.11 
through [5]6j the results for the double-dipole case are shown in Figures [577] through 
15.121 and the results for the triple-dipole are shown in Figures [5T3] through l5. 181 On 
each plot of the renormalization flow in Figures [STI through 15 . 1 81 there are multiple 
curves, each corresponding to different values of the upper bound of the fit window, 
'^n.max- ^ example points are selected in Figures [STI through 15. 181 to indicate 
the general size of the statistical error bars. 

It should be noted that none of the curves in Figures 15.11 through 15.181 is flat 
to within 1% accuracy. All of the fits have lattice data included beyond the PGR. 
Clearly, there is a well-defined intersection point in each plot. Also, the value of A 
at which the intersection point occurs is the same, even for different data sets, and for 
different c,. The tight groupings of the curve crossings lend credence to the notion 
of an intrinsic scale that can be interpreted as a finite size of the source of the pion 
dressings of the nucleon. This is a central result of the analysis. 

Using the method described in Chapter |4l the intersection point of the renor- 
malization flow curves for different values of j^^^ is estimated from Figures 15.11 
through 15.181 As an initial estimate, by inspection, a mean value for the optimal 
regularization scale of A^^p'^ ^1.3 GeV was obtained for the dipole, a value of 
-^doub ~ 1-0 GeV was obtained for the double dipole, and a value of A^^^^^ ^ 0.9 
GeV was obtained for the triple-dipole. These values differ because the regulators 
have different shapes, as evident in Figure ISTTl and thus different values of A^'^^^'^ are 
required to create a similar suppression of large loop momenta. In order to determine 
an estimate of the systematic uncertainty in an extrapolation due to the choice of reg- 
ularization scale A^'^^^'^, one should use a robust method for estimating the systematic 
uncertainty of A^'^^''^ itself. In the following section, a chi-square-style analysis will 
be introduced to fulfill this requirement. 

5.1.2 Analysis of Systematic Uncertainties 

The optimal regularization scale A^'^^^'^ can be obtained from the renormalization 
flow curves using a chi-square-style analysis. In addition, the analysis will allow the 
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Figure 5.1: Behaviour of co vs. A, based on 
JLQCD data. The chiral expansion is taken to or- 
der o(ot^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 
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Figure 5.4: Behaviour of C2 vs. A, based on 
JLQCD data. The chiral expansion is taken to or- 
der (m^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 
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Figure 5.2: Behaviour of cq vs. A, based on 
PACS-CS data. The chiral expansion is taken to 
order (m^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 
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Figure 5.5: Behaviour of C2 vs. A, based on 
PACS-CS data. The chiral expansion is taken to 
order (m^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 
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Figure 5.3: Behaviour of cq vs. A, based on 
CP-PACS data. The chiral expansion is taken to 
order (m^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 



Figure 5.6: Behaviour of C2 vs. A, based on 
CP-PACS data. The chiral expansion is taken to 
order (m^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 
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Figure 5.7: Behaviour of cq vs. A, based on 
JLQCD data. The chiral expansion is taken to or- 
der (m^) and a double-dipole regulator is used. 
A few points are selected to indicate the general 
size of the statistical error bars. 



Figure 5.10: Behaviour of C2 vs. A, based on 
JLQCD data. The chiral expansion is taken to or- 
der (m^) and a double-dipole regulator is used. 
A few points are selected to indicate the general 
size of the statistical error bars. 
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Figure 5.8: Behaviour of cq vs. A, based on 
PACS-CS data. The chiral expansion is taken to 
order O (ot^) and a double-dipole regulator is used. 
A few points are selected to indicate the general 
size of the statistical error bars. 
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Figure 5.11: Behaviour of C2 vs. A, based on 
PACS-CS data. The chiral expansion is taken to 
order O (m^) and a double-dipole regulator is used. 
A few points are selected to indicate the general 
size of the statistical error bars. 
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Figure 5.9: Behaviour of cq vs. A, based on CP- 
PACS data. The chiral expansion is taken to order 
(m^) and a double-dipole regulator is used. A 
few points are selected to indicate the general size 
of the statistical error bars. 



Figure 5.12: Behaviour of C2 vs. A, based on 
CP-PACS data. The chiral expansion is taken to 
order (m^) and a double-dipole regulator is used. 
A few points are selected to indicate the general 
size of the statistical error bars. 
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Figure 5.13: Behaviour of cq vs. A, based on 
JLQCD data. The chiral expansion is taken to or- 
der (m^) and a triple-dipole regulator is used. A 
few points are selected to indicate the general size 
of the statistical error bars. 



Figure 5.16: Behaviour of C2 vs. A, based on 
JLQCD data. The chiral expansion is taken to or- 
der (m^) and a triple-dipole regulator is used. A 
few points are selected to indicate the general size 
of the statistical error bars. 





1 1 


1 1 1 






m = = 0.49 GeV= 






















1 1 





0.0 0.4 



0.8 1.2 l.( 
A (GeV) 



2.0 2.4 




0.32 GeV^ 
^ = 0.49 GeV^ 



0.0 0.4 



0.8 1.2 1.6 2.0 2.4 
A (GeV) 



Figure 5.14: Behaviour of cq vs. A, based on 
PACS-CS data. The chiral expansion is taken to 
order (m^) and a triple-dipole regulator is used. 
A few points are selected to indicate the general 
size of the statistical error bars. 




A (GeV) 



Figure 5.17: Behaviour of C2 vs. A, based on 
PACS-CS data. The chiral expansion is taken to 
order (m^) and a triple-dipole regulator is used. 
A few points are selected to indicate the general 
size of the statistical error bars. 
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Figure 5.15: Behaviour of cq vs. A, based on 
CP-PACS data. The chiral expansion is taken to 
order (m^) and a triple-dipole regulator is used. 
A few points are selected to indicate the general 
size of the statistical error bars. 



Figure 5.18: Behaviour of C2 vs. A, based on 
CP-PACS data. The chiral expansion is taken to 
order o(m^) and a triple-dipole regulator is used. 
A few points are selected to indicate the general 
size of the statistical error bars. 
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extraction of an estimate of the variance for A**'^^'^. The function x^^jr defined below 
allows easy identification of the intersection points in the renormalization flow plots, 
and a range associated with this central regularization scale. This function simply 
measures the degree to which the renormalization flow curves match. 

The first step is to plot x^^jr against a variety of regularization scales. The value 
of c is given by the weighted mean formula, evaluated separately for each renormal- 
ized coefficient c (with error 6c) and regularization scale A: 

,-(^)^£L,c(.;A)/(8.(.;A))\ (5 2, 

E5.,l/(8c(i;A))2 

The following x^iof treats relevant degrees of freedom as the extracted chiral coeffi- 
cients with differing values of m\^^y^: 

2 _ 1 A (c(/;A)-c-(A))^ 
W-;^L (g^(..^))2 ' (5.3) 



(=1 



that is, / corresponds to fits with differing values of m^^^g,^- 

The can be calculated as a function of the regularization scale A for each of 
the renormalization plots of Figures 15.1 I through lS .181 This will indicate the spread of 
the extrapolated values at each value of A. In the case of the PACS-CS data, the min- 
imum of the x^^^j- curve will be at the intersection point of the two curves. In the case 
of the JLQCD and CP-PACS data, with more than two curves, there is an interaction 
region on each plot, over a narrow window of A. The minima of x^^^y^ will indi- 
cate the value of A that obtains the best agreement among the renormalization flow 
curves. This central value of A will be taken to be the optimal regularization scale. 
The upper and lower bounds of A obey the condition x^^of < X^o/ min + V i^of) . For 
each of the low-energy coefficients cq and C2, the Xdof curves for a dipole regulator 
are shown in Figures 15.191 through I5.24[ the Xdof curves for the double-dipole case 
are shown in Figures [5. 251 through l5. 301 and the Xdof curves for the triple-dipole are 
shown in Figures [5.3 II through I5.36[ These plots indicate that there exists a statisti- 
cally significant optimal regularization scale at this chiral order, for these data sets. 
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Figure 5.19: Behaviour of x^„/ for cq vs. A, 
based on JLQCD data. The chiral expansion is 
taken to order 0(?m^) and a dipole regulator is 
used. 



Figure 5.22: Behaviour of x^^/ for C2 vs. A, 
based on JLQCD data. The chiral expansion is 
taken to order O(m^) and a dipole regulator is 
used. 
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Figure 5.20: Behaviour of for cq vs. A, 
based on PACS-CS data. The chiral expansion 
is taken to order o(m^) and a dipole regulator is 
used. 



Figure 5.23: Behaviour of x^o/ for C2 vs. A, 
based on PACS-CS data. The chiral expansion 
is taken to order o(m^) and a dipole regulator is 
used. 
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Figure 5.21: Behaviour of for cq vs. A, 
based on CP-PACS data. The chiral expansion 
is taken to order (m^) and a dipole regulator is 
used. 



Figure 5.24: Behaviour of Xjof f°r ^2 vs. A, 
based on CP-PACS data. The chiral expansion 
is taken to order O (m^) and a dipole regulator is 
used. 
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Figure 5.25: Behaviour of x^/^/ for cq vs. A, 
based on JLQCD data. The chiral expansion is 
taken to order (m^) and a double-dipole regula- 
tor is used. 



Figure 5.28: Behaviour of for C2 vs. A, 
based on JLQCD data. The chiral expansion is 
taken to order (m^) and a double-dipole regula- 
tor is used. 
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Figure 5.26: Behaviour of for cq vs. A, 
based on PACS-CS data. The chiral expansion is 
taken to order (m^) and a double-dipole regula- 
tor is used. 



Figure 5.29: Behaviour of x^o/ for C2 vs. A, 
based on PACS-CS data. The chiral expansion is 
taken to order (m^) and a double-dipole regula- 
tor is used. 
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Figure 5.27: Behaviour of for cq vs. A, 
based on CP-PACS data. The chiral expansion is 
taken to order (m^) and a double-dipole regula- 
tor is used. 



Figure 5.30: Behaviour of for C2 vs. A, 
based on CP-PACS data. The chiral expansion is 
taken to order (;«^) and a double-dipole regula- 
tor is used. 
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Figure 5.31: Behaviour of for cq vs. A, 
based on JLQCD data. The chiral expansion is 
taken to order (m^) and a triple-dipole regulator 
is used. 



Figure 5.34: Behaviour of x^^^ for C2 vs. A, 
based on JLQCD data. The chiral expansion is 
taken to order (m^) and a triple-dipole regulator 
is used. 
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Figure 5.32: Behaviour of for cq vs. A, 
based on PACS-CS data. The chiral expansion is 
taken to order (m^) and a triple-dipole regulator 
is used. 



Figure 5.35: Behaviour of x^o/ for C2 vs. A, 
based on PACS-CS data. The chiral expansion is 
taken to order (m^) and a triple-dipole regulator 
is used. 
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Figure 5.33: Behaviour of for co vs. A, 
based on CP-PACS data. The chiral expansion is 
taken to order (m^) and a triple-dipole regulator 
is used. 



Figure 5.36: Behaviour of xfj^j for C2 vs. A, 
based on CP-PACS data. The chiral expansion is 
taken to order (m^) and a triple-dipole regulator 
is used. 
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Furthermore, for each data set and regulator functional form, there is an agreement 
between the cq and C2 analyses as to the value of this optimal scale. This provides 
evidence of the existence of an intrinsic scale embedded in that lattice data. 

5.1.3 Effects at Higher Chiral Order 

Consider the determination of cq and C2 as a function of the regularization scale A, 
for a higher chiral order (^(m^logm^t). As an example, the results for PACS-CS 
and CP-PACS data are shown in Figures 15.371 through 15.401 In this case, no clear 
intersection points in the renormalization curves can be found, and so one is unable 
to specify an optimal regularization scale. This certainly should be the case when 
working with data entirely within the PCR, because all renormalization procedures 
would be equivalent (to a prescribed level of accuracy) and so there could be no 
optimal scale. It has been demonstrated, however, that the data sets used in this study 
extend beyond the PCR. This is further verified by considering the evident scale- 
dependence of Co and C2 in Figures [5.37l through l5.40[ The fact that cq and C2 change 
over the range of A values indicates that the data are not inside the PCR where the 
renormalization must be scale-independent. Furthermore, since no preferred scale is 
revealed, any choice of A appears equivalent at this order. While it is encouraging 
that the scheme-dependence has been weakened by working to higher order, it must 
be recognized that there is a systematic error associated with the choice of A. In the 
case of the CP-PACS results shown in Figures [5.381 and [5.401 it can be seen that the 
statistical errors are substantially smaller than the systematic error associated with a 
characteristic range, Aiower < A < oo, where Aiower is the lowest reasonable value of 
A, taken to be 0.6, 0.4 and 0.3 GeV for the dipole, double-dipole and triple-dipole 
regulator, respectively, as discussed in Section [4.2. 11 

Since it is difficult to identify the optimal regularization scale at this chiral or- 
der, the results for chiral order (m^) will be chosen to demonstrate the process of 
handling the existence of an optimal regularization scale in lattice QCD data. Values 
of A^*^^^*^ for different data sets and regulators, using chiral order o(m^), are given 
in Table 15.11 This table simply summarizes the central values from Figures 15.191 
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Figure 5.37: Behaviour of cq vs. A, based on 
PACS-CS data. The chiral expansion is taken to 
order (m^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 
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Figure 5.39: Behaviour of C2 vs. A, based on 
PACS-CS data. The chiral expansion is taken to 
order (m^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 
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Figure 5.38: Behaviour of cq vs. A, based on 
CP-PACS data. The chiral expansion is taken to 
order (m^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 



Figure 5.40: Behaviour of C2 vs. A, based on 
CP-PACS data. The chiral expansion is taken to 
order (m^) and a dipole regulator is used. A few 
points are selected to indicate the general size of 
the statistical error bars. 
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Table 5.1: Central values of A in GeV, taken from the X^^j analysis for cq and C2, based on JLQCD, 
PACS-CS and CP-PACS data. 







regulator form 




optimal scale 


dipoL 


e 


double 


triple 




A scale 
^Vo JLQCD 

A scale 

JLQCD 


1.44^ 
1.40^ 


-0.18 
-0.18 

-0.02 
-0.03 


i.08i»:|| 


0.96^ 
0.94^ 


-0.09 
-0.09 

-0.01 
-0.02 


A scale 

^V(),PACS-CS 
A scale 

^S,PACS-CS 


1.21^ 
1.21^ 


-0.66 
-0.82 

-0.18 
-0.18 


93+0-^1 
0.93l«:|i 


0.83^ 
0.83^ 


-0.35 
-0.50 

-0.10 
-0.10 


A scale 

^Vo,CP-PACS 
A scale 

^V?,CP-PACS 


1.20^ 
1.19^ 


-0.10 
-0.10 

-0.02 
-0.01 


0.97l«:«j 





through [5.36[ Such excellent agreement between the cq and C2 analyses is remark- 
able, and indicative of the existence of an intrinsic scale in the data. There is also 
consistency among independent data sets. It is important to realize that the value 
of A^'^^^'^ is always the order of ~ 1 GeV, not 10 GeV, nor 100 GeV; nor is it in- 
finity. However, in calculating the systematic uncertainty in the observables cq, C2, 
and the nucleon mass at the physical point due to the optimal regularization scale at 
order o{m'l\ogmTi), two methods are provided. Firstly, the upper and lower bounds 
from the x^^^^ analysis at order (m^) are used to constrain A, and taken to be an 
accurate estimate of the systematic uncertainty in the contributions of higher-order 
terms. Secondly, variation of the observables across the aforementioned character- 
istic range of scale values, Aiower < A < oo are used. The results from both of these 
methods are displayed in Table 15.21 

The final results for the calculation of the renormalized coefficients cq, C2 and 
the nucleon mass Mf^ extrapolated to the physical point (mTt^phys = 140 MeV) are 
summarized in Table 15.31 In this table, the nucleon mass is calculated at the optimal 
scale A^*^^^*^, which is the average of A^.^J^^'^ and A^.^^'*^ for each data set. The extrap- 
olations are performed at lattice sizes relevant to each data set: L^J^^^^ = 1.9 fm, 
^e^rap *^^ = 2.9 fm and L^^^^^^^ = 2.8 fm. The estimate of the statistical error is 
quoted in the first pair of parentheses, and the systematic error, obtained from the 
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Table 5.2: Results at O (m^logmTt) for the systematic error due to the optimal regularization scale, 
calculated using two methods, for the values of cq (GeV), C2 (GeV^') and the nucleon mass Mn 
(GeV) extrapolated to the physical point (m^^phys = 140 MeV). The first number in each column is 
the systematic error due to the optimal regularization scale using the upper and lower bound from the 
X^i^^^ analysis at order o(m^). The second number is the systematic error due to the intrinsic scale 
across the whole range of A values from the lowest reasonable value (A = Aio^er) obtained from the 
pseudodata analysis, to the asymptotic value (A = oo). 



regulator form 



sys. err. 


dipole 


double 


triple 


JLQCD 

O Cq 


0.001 


0.009 


0.001 


0.013 


0.001 


0.016 


2A^PACS-CS 

O Cq 


0.005 


0.006 


0.005 


0.010 


0.006 


0.012 


sAXP-PACS 

O Cq 


0.002 


0.024 


0.002 


0.037 


0.002 


0.045 


gA^ILQCD 


0.02 


0.31 


0.03 


0.38 


0.01 


0.48 


5A^PACS-CS 


0.18 


0.25 


0.16 


0.33 


0.14 


0.43 


gA^CP-PACS 


0.02 


0.40 


0.02 


0.58 


0.02 


0.73 


gA^JLQCD 


0.0004 


0.0051 


0.0003 


0.0073 


0.0003 


0.0090 


s;Aj^PACS-CS 


0.0022 


0.0030 


0.0025 


0.0046 


0.0025 


0.0058 


gA^CP-PACS 


0.0012 


0.0175 


0.0013 


0.0270 


0.0014 


0.0326 



number of values used, is quoted in the second pair of parentheses. Two dif- 
ferent weighted means are calculated. One incorporates the systematic error in the 
optimal regularization scale using the upper and lower bound from the X^^/ ^^^Y" 
sis at order o(m^). The other incorporates the systematic error due to the optimal 
regularization scale across the whole range of A values, from the lowest reasonable 
value (A = Aiower) obtained from the pseudodata analysis, to the asymptotic value 
(A = oo). The weighted means also include an estimate of the systematic error in 
the choice of regularization scale. All errors are added in quadrature. The lightest 
four data points from each of JLQCD, PACS-CS and CP-PACS lattice QCD data are 
used, and the nucleon mass is calculated at the scale determined by the data. 
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Table 5.3: Results at 0(mJ logmjt) for the values of cq (GeV), C2 (GeV~^) and the nucleon mass 
Mn (GeV) extrapolated to the physical point (wn.phys = 140 MeV). WM is the weighted mean of each 
row. The nucleon mass is calculated at the optimal scale A^'^"'^, which is the average of A^^''''^ and 
y^scaie Jqj. gj^pj^ ^^^^ -pjjg extrapolations are performed at lattice sizes relevant to each data set: 

4x^° = 1-9 fm. ^e^sttap"*^^ = 2-9 fm and i-extop^*^^ = 2-8 fm. The estimate of the statistical error 
is quoted in the first pair of parentheses, and the systematic error, obtained from the number of 
values used, is quoted in the second pair of parentheses. Two separate weighted means are calculated 
for each row. WM(1) incorporates the systematic error in the intrinsic scale using the upper and lower 
bound from the x^^^^ analysis at order o(m^). The WM(2) incorporates the systematic error due to 
the intrinsic scale across the whole range of A values, from the lowest reasonable value (A = Aiower) 
obtained from the pseudodata analysis, to the asymptotic value (A = <»). The weighted means also 
include an estimate of the systematic error in the choice of the regulator functional form. All errors 
are added in quadrature. Note that any order o (a) errors have not been incorporated into the total 
error analysis. 



regulator form 



parameter 


dipole 


double 


triple 


WM(1) 


WM(2) 


JLQCD 


0.873(18)(16) 


0.875(17)(16) 


0.891(17)(16) 


0.880(29) 


0.879(32) 


PACS-CS 


0.900(51)(15) 


0.899(51)(14) 


0.898(51)(14) 


0.899(53) 


0.899(55) 


^CP-PACS 


0.924(3) (8) 


0.914(3)(7) 


0.918(3)(7) 


0.918(13) 


0.920(37) 


JLQCD 
^2 


3.09(9)(11) 


3.18(9)(12) 


3.20(9)(11) 


3.16(18) 


3.14(43) 


PACS-CS 
^2 


3.06(32)(15) 


3.15(31)(14) 


3.17(31)(14) 


3.13(39) 


3.12(49) 


^CP-PACS 
^2 


2.54(5) (4) 


2.70(5) (2) 


2.71(5)(3) 


2.66(18) 


2.61(60) 


^JLQCD 

■'*^A',phys 


1.02(2)(9) 


1.02(2)(9) 


1.02(2)(9) 


1.02(9) 


1.02(9) 


»^PACS-CS 


0.967(45) (43) 


0.966(45) (43) 


0.966(45) (43) 


0.966(62) 


0.966(62) 


;,^CP-PACS 

^'^Af,phys 


0.982(2) (40) 


0.975(2) (43) 


0.978 (2) (42) 


0.979(43) 


0.979(50) 



5.2 Summary and Specific Issues for the Nucleon Mass 

Since the chiral expansion is only convergent within the PCR, a renormalization 
scheme such as finite-range regularization should be used for current lattice QCD 
results, which typically extend beyond the PCR. It was found that renormalization 
scheme-dependence occurs when lattice QCD data extending outside the PCR are 
used in the extrapolation. This has provided a new quantitative test for determining 
whether lattice QCD data he within the PCR. 

The optimal regularization scale A^'^^® was selected as the scale at which the 
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renormalized coefficients are independent of the upper bound of the fit domain, 
'^nmax- This also means that the renormalized coefficients must not be identified 
with their asymptotic values at large A. It is also apparent that extremely low val- 
ues of A cause a breakdown of the finite-range renormalization. The cut-off scale 
associated with an ultraviolet regulator must be large enough for the loop integral 
contributions to be finite, so that the chiral physics is not suppressed. 

The mean value of the optimal regularization scale for both the cq and C2 anal- 
yses across each data set is A^^^^^ ^1.3 GeV for the dipole form, A^'^^^ ^1.0 GeV 
for the double-dipole form and A.^^^-^^ ~ 0.9 GeV for the triple-dipole form. Each 
functional form naturally leads to a different value of optimal regularization scale 
due to its different shape of attenuation, as shown in Figure [3TT1 The value of A^^p'® 
is of particular interest in this investigation. In Chapter|71 the magnetic moment and 
the electric charge radius of the nucleon are analyzed with the same procedure, and 
using a dipole regulator. If an optimal regularization scale can be obtained for these 
electromagnetic properties of the nucleon, a comparison can be made with the opti- 
mal regularization scale from the analysis of the nucleon mass, to determine whether 
there exists an intrinsic scale in the nucleon. If the optimal regulators in each case 
are consistent with each other, this suggests the existence of a well-defined intrin- 
sic energy scale in the nucleon-pion interaction. Nevertheless, a robust method for 
accomplishing a chiral extrapolation with a reliable and systematic estimate in the 
uncertainty has been provided. 

In the next chapter, the procedure developed for obtaining an optimal regulariza- 
tion scale and performing a reliable chiral extrapolation is tested, by analyzing the 
quenched p meson mass: an observable for which there does not exist an experimen- 
tal value. This serves to demonstrate the ability of the extrapolation scheme to make 
predictions without prior bias. 



Chapter 6 

Results for the Mass of the Quenched 
p Meson 



"A rigorous theory must begin by specifying the attributes that make a given exper- 
imental device into a measuring instrument." (Omnes, R. 2002. Quantum Philoso- 
phy: Understanding and Interpreting Contemporary Science p. 209) [|Omn02ll 

The quenched p meson mass offers a unique test case for the identification of an 
intrinsic scale, and subsequent extrapolation scheme. It serves to demonstrate the 
ability of the procedure to make predictions with reduced phenomenological bias, 
and also to highlight the difference between quenched and full quantum chromo- 
dynamics (QCD) in making extrapolations of an observable. By using the method 
developed in Chapter |4l an extrapolation is performed using quenched lattice QCD 
data that extend outside the power-counting regime (PCR). 

In chiral effective field theory (xEFT), the diagrammatic formulation can be used 
to identify the major contributions to the p meson mass in quenched QCD (QQCD) 
[|CR98[IAAL"'"06| . The leading-order diagrams are the double and single r|' hairpin 



diagrams as shown in Figures 16.11 and I6.4[ respectively. The constant coefficients 
of these loop integrals are endowed with an uncertainty to encompass the possible 
effects of smaller contributions to order (m^) . Interactions with the flavour-singlet 
Tj' are the most important contributions to the p meson mass in QQCD. This is an 
artefact of the quenched approximation, where the Tj' also behaves as a pseudo- 



83 



Chapter 6. Results for the Mass of the Quenched p Meson 



Hall 84 



/ 



\ 



A. 



1 i 



Figure 6.1: Double hairpin Tj' diagram. 



Figure 6.4: Single hairpin r|' diagram. 




d " " d 

Figure 6.2: Double hairpin quark flow diagram. 




d " " d 

Figure 6.5: Single hairpin quark flow diagram. 




Figure 6.3: Alternative double hairpin quark 
flow diagram. 



Goldstone boson, having a mass that is degenerate with the pion. The dressing of 
the p meson by the r\' field is illustrated in Figures [6^ through l6.5[ Since the hairpin 
vertex must be a flavour-singlet, the mesons that can contribute are the r\' meson, and 
the (0 meson. The contributions from the CO meson are insignificant due to OZI sup- 
pression and the small p-CO mass splitting. However, in QQCD, the T]' loop behaves 
much as a pion loop, yet with a slightly modified propagator. 

In full QCD however, this would not be the case. The r\' masses are large com- 
pared to the pion, and the propagators of the r\' meson are suppressed due to their 
large denominators. If the r\' propagator in full QCD is expanded out, the terms 
can be summed as a geometric series and expressed in closed form, as a function of 
some massive coupling constant Mq between the disconnected quark loops and pion 
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momentum k, as argued in AUton [lAAL+OSll 



Ml 



■X 



Ml 



+ 



Ml 



(6.1) 



Ml 



1 + 



Ml 



k^ + 



m 



31/ J 



-1 



(6.2) 



k^ + ml+Ml k^ + m^, 



(6.3) 



However in QQCD, the first two terms of Equation (16.11) form the whole r\' prop- 
agator, since they alone correspond to the absence of disconnected loops, as shown 



in Figure [6^ 




Figure 6.6: Diagrammatic representation of r|' propagator terms. 
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6.1 Renormalization of the Quenched p Meson Mass 
6.1.1 Chiral Expansion of the Quenched p Meson Mass 

The p meson mass extrapolation formula in QQCD can be expressed in a form that 
contains an analytic polynomial in plus the chiral loop integrals (Tfi): 

Q = a^ + a^ml + 4 m\ + E^,^, (m^ , A) + E^, {m\, A) + (m^) . (6.4) 

The coefficients axe the 'residual series' coefficients, which correspond to direct 
quark-mass insertions in the underlying Lagrangian of chiral perturbation theory 
(%PT). However, the non-analytic behaviour of the expansion arises from the chi- 
ral loop integrals. Upon renormalization of the divergent loop integrals, these will 
correspond with low-energy constants of the quenched effective field theory. The 
extraction of these parameters from lattice QCD results will follow the same course 
as provided in Chapter IH 

By convention, the non-analytic terms from the double and single hairpin inte- 
grals are Xif^n and Xsm^, respectively. The coefficients %i and %3 of each integral are 
scheme-independent constants that can be estimated from phenomenology. That is, 
they can be expressed purely in terms of known constants from experiment, such as 
the pion decay constant fi^ = 92A MeV, and a variety of parameters obtained from 
the underlying effective Lagrangian, as described in Section I6.1.3[ The low-order 
expansion of the loop contributions takes the following form: 

= ^o'"^' + Xim^ + bf^'fnl + I't'^'ml + b'i'^'ml + (m^) , (6.5) 
4 = ^o' + 4^1 + llf^l + bim"^ + {ml) . (6.6) 

The coefficient X3 is obtained by adding the contributions from both integrals, %3 = 
%3^ +X3 • As before, each integral has a solution in the form of a polynomial 
expansion analytic in plus non-analytic terms, of which the leading-order term is 
of greatest interest. In order to achieve an extrapolation based on an optimal finite- 
range regularization (FRR) scale, once again the scale-dependence of the low-energy 
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expansion must be removed through renormalization. The renormalization program 
of FRR combines the scheme-dependent bj coefficients from the chiral loops with 
the scheme-dependent ai coefficients from the residual series at each chiral order /. 
The result is a scheme-independent coefficient c,: 

co = a^ + ^3'^' + ^S', (6.7) 
C2 = fl^ + Z75''i' + ^5', (6.8) 
C4 = a^ + bi"^' + b'i , Qtc. (6.9) 

That is, the underlying af- coefficients undergo a renormalization from the chiral loop 
integrals. The renormalized coefficients are an important part of the extrapolation 
technique. In this chapter, a stable and robust determination of the parameters cq, C2 
and C4 forms the core of the method for determining an optimal scale A^'^"^^ of the 
mass of the p meson. 



6.1.2 Chiral Loop Integrals 

The loop integrals can again be expressed conveniently by taking the non-relativistic 
limit and performing the pole integration for ^o- Renormalization is achieved by 
subtracting the relevant terms in the Taylor expansion of the loop integrals, and ab- 
sorbing them into the corresponding low-energy coefficients, c;: 

V^"^-^)- 37t2 1^^ i^k^^mlf 

-bf-bfml-bfml (6.10) 

- 27i2 y ^ ^2 + ^2 ^0 ^2 

-biml. (6.11) 

The tilde (~) denotes that the integrals are written out in renormalized form to chiral 
order o(m^). The coefficients and %^ are related to the coefficients of the 



Chapter 6. Results for the Mass of the Quenched p Meson 



Hall 88 



leading-order non-analytic terms by the following: 

Xi=^o%W' (6.12) 
Xs^xJ'^^' + xJ'^AoXtiy+Xti/. (6.13) 

In choosing the form of regulator, one must be cautious to avoid any extra 
scheme-dependent, non-analytic terms that might occur at working chiral order. For 
example, consider pseudodata created from the lattice QCD simulations from the 
Kentucky Group, using a dipole regulator created using the scale: Ac = 0.8 GeV. 
The renormalization flow curves are shown in Figures 16.71 through 16.91 The dipole 
regulator induces non-analytic terms proportional to and in the loop inte- 
gral expansion formulae. By writing out the regulator-dependence explicitly in the 
coefficients bj, the following equations are obtained: 

rri'Ti' ^ri'ri' rrj'ri' ■^r['r\' 

Ee;J;P(fc ; A) = A^|]'^' +Ximn + -j^ml + -j^ml + + + o {m\) , 

(6.14) 

Y.^^'\k-K)=K^^ ^Kh\ml^X3ml + ^m\^j^ml^0{mi). (6.15) 

Clearly, the renormalization flow is compromised by the extra non-analytic terms 
appearing at such a close chiral order to the fit parameters. Though it is possible to 
provide additional fit parameters and to contain the contribution from these 
terms, there are often not enough available lattice simulation results to constrain all 
coefficients. Instead, a more effective approach is to choose a regulator functional 
form such that the extra non-analytic terms do not appear in the chiral expansion. By 
selecting an multiple-dipole regulator corresponding to a choice of n > 2 in Equation 
(13.1051) in Chapter [3l the suppression of additional non-analytic terms below the 
working chiral order O(m^) is assured. If one also decides to remove extra m| 
terms, a triple-dipole is sufficient to remove additional non-analytic terms below 
chiral order o(m^). The renormalization flow curves for pseudodata created with 
a double-dipole are shown in Figures 16.101 through I6.12[ The renormalization flow 
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Figure 6.7: Behaviour of cq vs. A based on infinite-volume pseudodata created with a dipole 
regulator at Ac = 0.8 GeV (based on Kentucky Group data). Each curve uses pseudodata with a 
different upper value of pion mass 
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Figure 6.8: 
regulator at Ac 
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Behaviour of C2 vs. A based on infinite-volume pseudodata created with a dipole 
= 0.8 GeV (based on Kentucky Group data). 
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regulator at Ac = 0.8 GeV (based on Kentucky Group data). 
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curves for pseudodata created with a triple-dipole regulator are shown in Figures 
l6.13l through l6.15[ In both cases, the pseudodata are created using the scale: Ac = 0.8 
GeV. 

With the loop integrals specified, Equation (|6.4I) can be rewritten in terms of the 
renormalized coefficients c,: 

mp,Q = co + C2m^ + C4m^ + E^,^, (m^; A) +I^,{m^;A) + (m^) (6.16) 
= CQ + XimTt + C2ml + X^ml + C4m^ + (m^) . (6. 17) 

Equation (16.161) is the extrapolation formula for q at infinite lattice volume. The 
fit coefficients are cq, C2 and C4; and mp^g is obtained by taking the square root of 
either Equation (|6.16l) or (|6.171) . 



6.1.3 Scheme-Independent Coefficients 

The convention used for defining the values of %i, %3 and the various coupling con- 
stants that occur in each, follows Booth nBCF97L For the possible different values 
that coupling constants can take, the definitions from Chow & Rey [ICR98L Armour 
et.al. llAAL"'"06l and Sharpe [|Sha971l are used. The types of vertices available are 



displayed in Figure I6.16[ where the couplings g2 and occur explicitly in the two 
diagrams considered here. Booth suggests naturalness for ^2 ~ 1 , and that g4 ~ 1 /Nc 
[IBCF97II . These quenched coupling constants can be connected to the experimental 
value of gatpTi as per Lublinsky [|Lub97ll by the relation: 



g2 = Ij^gapnfn, (6.18) 

where goipn = 

14 ±2 GeV- 1 and the pion decay constant is again taken to be /jt ~ 
92.4 MeV. Thus g2 is chosen to be 0.65 ±0.09 GeV and g^ is chosen to be ^2/3. The 
coupling between the separate legs of the double hairpin diagram are approximated 
by the massive constant Mq oc m^,. The next-order correction to Mq in momentum k 
defines the coupling to be — Mq +Aok^. These constants can be connected to the full 
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dipole regulator at Ac = 0.8 GeV (based on Kentucky Group data). 
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Figure 6.15: Behaviour of 04 vs. A based on infinite-volume pseudodata created with a triple- 
dipole regulator at Ac = 0.8 GeV (based on Kentucky Group data). 
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§1 



Figure 6.16: Coupling types following convention introduced by Booth IIBCF97I . 



QCD r\' meson mass m^/ by considering the geometric series of terms as illustrated 
earlier, in Figure [6^ As a result, Mq is taken to be 0.6 ± 0.2 GeV^ and Aq is taken 
to be ± 0.2. The central values of each quantity are used in the final analysis. 

The coefficients Xr\'r\' and Xr[' can be specified in terms of the relevant coupling 
constants: 

X,. = -2™p||, (6.19) 

o 

where the couplings are defined relative to trip, representing the p meson mass in the 
chiral limit, which is taken to be 770 MeV. 

The finite-volume version of Equation (16.161) can thus be expressed: 

mp^Q = Co + C2m^ + C4m^ + (E^,^,(m^; A) + 5^^ (m^; A )) 

+ {f^,{ml;A) + diy''{ml;A')) + 0{m^^) . (6.20) 
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6.2 Extrapolating the Quenched p Meson Mass 
6.2.1 Renormalization Flow Analysis 

The data displayed in Figure [QT] are split into two parts. All the data points to the 
left of the solid vertical line are unused in the extrapolation and kept in reserve. This 
is so that the extrapolation can be checked against these known data points. The data 
points to the right of the solid vertical line are used for extrapolation. The full set 
of data is also listed in Appendix O Table IC.41 Note that in QQCD, the simulation 
results are correlated. The correlations have been taken into account in all fits and 
extrapolations. 

In order to produce an extrapolation to each test value of m^, an FRR scale A 
must be selected. As an example, one can choose a triple-dipole regulator at A = 1 .0 
GeV. By using Equation (16.201) . finite- and infinite- volume extrapolations are shown 
in Figure l6T8l The values of selected for the finite- volume extrapolations exactly 
correspond with the missing low-energy data points set aside earlier. The physical 
point m\ = 0.0196 GeV^ is included as well. 

Now the regularization scale-dependence of low-energy coefficients cq, C2 and C4 
is investigated for various upper limits of range of pion masses. The renormalization 
of these low-energy coefficients is considered across a range of A values. Each 
renormalization flow curve corresponds to a different value of maximum pion mass, 
f^n,max- Thus the behaviour of the renormalization of the low-energy coefficients 
can be examined as the lattice data set is extended further outside the PCR. Figures 
16. 191 through 16.211 show the renormalization flow curves for each of cq, and C4. 
Each data point plotted has an associated error bar, but for the sake of clarity only a 
few points are selected to indicate the general size of the statistical error bars. Using 
the procedure described in Chapter |4l the optimal regularization scale is identified 
by the value of the regularization scale that minimizes the discrepancies among the 
renormalization flow curves. This indicates the scale at which the renormalization of 
each Ci is least sensitive to truncation of the data. Physically, this value of A can be 
associated with an intrinsic scale related to the size of the source of the pion cloud. 
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Figure 6.17: Quenched lattice QCD data for the p meson mass provided by the Kentucky Group. 
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By examining Figures 16.191 through 16.211 increasing j^^^ leads to greater 
scale-dependence in the renormalization, since the data sample lies further from 
the PGR. Since the effective field theory is calculated to a finite chiral order, com- 
plete scale-independence across all possible A values will not occur in practice. An 
asymptotic value is usually observed in the renormalization flow as A becomes large, 
indicating that the higher-order terms of the chiral expansion are effectively zero. 
However, these asymptotic values of the coefficients are poor estimates of their cor- 
rect values, as previously demonstrated in the pseudodata analysis in Chapter |4l 
Instead, the best estimates of the low-energy coefficients he in the identification of 
the intersection point of the renormalization flow of these coefficients. It is also of 
note that, for small values of A, the FRR scheme breaks down, as observed for the 
nucleon mass in Section l4.2.1[ The regularization scale must be at least large enough 
to include the chiral physics being studied. 



6.2.2 Intrinsic Scale and Systematic Uncertainties 

The optimal regularization scale A^*^^^*^ can be obtained from the renormalization 
flow curves using a chi-square-style analysis. In addition, the analysis will allow 
the extraction of a variance for A^'^^''^. Knowing how the data are correlated, the 
systematic errors from the coupling constants and A^'^^^'^ will be combined to obtain 
an error estimate for each extrapolation point. Of particular interest are the values 
of mp Q at those values of that are explored in the lattice simulations, but are 
excluded in the chiral extrapolation. The function Xjof constructed in the same 
way as Equations (15.21) and (15.31) . 

The Xdof plots using a triple-dipole regulator are shown in Figures [6.22l through 
16.241 The optimal regularization scale A'*'^^''^ is taken to be the central value A'^'^"^''^' 
of each plot. The upper and lower bounds of A obey the condition x^^j: < x^„f + 
1 / (dof). The results for the optimal regularization scales obtained from analyzing 
each low-energy coefficient, and their associated upper and lower bounds, are shown 
in Table 16.11 It is remarkable that each low-energy coefficient leads to the same 
optimal value of A, i.e. Acentrai = 0.67 GeV. By averaging the results among cq, C2, 
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Figure 6.22: x^^^j- for cq versus A, corresponding to the renormalization flow curves displayed in 
Figure l6?T9l based on Kentucky Group data. A triple-dipole regulator is used. 
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displayed in Figure ISTSTI based on Kentucky Group data. A triple-dipole regulator is used. 
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Table 6.1: Values of the central, upper and lower regularization scales, in GeV, obtained from the 
X^of analysis of co, C2 and 04, displayed in Figures l6.22l through l6.24l 



scale (GeV) 


CO (FigiL22D 


(Fi2J6.23l) 


C4 (Figl624D 


A-central 


0.67 


0.67 


0.67 


A-upper 


0.78 


0.75 


0.75 


A-lower 


0.58 


0.59 


0.60 



and C4, the optimal regularization scale Ascaie for the quenched p meson mass can be 
calculated for this data set: Ascaie = 0-67lo Qg GeV. 

The result of the final extrapolation, using the estimate of the optimal regular- 
ization scale Ascaie = 0-67lQ Qg GeV, and using the initial data set to predict the 
low-energy data points, is shown in Figure I6.25[ The extrapolation to the physical 
point obtained for this quenched data set is: w^p''g('«^phys) ~ ^•'^25lQ Q4g GeV, an 
uncertainty of less than 6%. 

Each extrapolation point displays two error bars. The inner error bar corresponds 
to the systematic uncertainty in the parameters only, and the outer error bar corre- 
sponds to the systematic and statistical uncertainties of each point added in quadra- 
ture. Also, the infinite-volume extrapolation curve is displayed in order to illustrate 
the effect of finite- volume corrections to the loop integrals. 

In Figure [6l26l the extrapolation predictions are compared against the actual sim- 
ulation results, which were not included in the fit. Both the extrapolations and the 
simulation results display the same non-analytic curvature near the physical point. 
Figure 16.271 shows the data plotted with error bars correlated relative to the lightest 
data point in the original set, = 0.143 GeV^. To highlight the importance of 
this application of an extended %EFT, a simple linear fit is included in Figure [6!27l 
By ignoring low-energy chiral physics, the linear fit is statistically incorrect at the 
physical point. All of the missing original data points are consistent with the ex- 
trapolations' systematic uncertainties. After statistical correlations are subtracted, 
the extrapolated points correspond to an error bar almost half the size of that of the 
lattice data points. In order to match this precision at low energies, the time required 
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Figure 6.26: Comparison of chiral extrapolation predictions (blue diamond) with Kentucky Group 



data (red cross). Extrapolation is performed at Ascaie — 0-67^o og GeV, and using the optimal number 
of data points, corresponding to jj^^x = 0.35 GeV^. The inner error bar on the extrapolation points 
represents purely the systematic error from parameters. The outer error bar represents the systematic 
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'"^,max = 0.35 GeV 



in lattice simulations would increase by approximately four times. 

In order to check if scheme-independence is recovered using data within the 
PCR, the low-energy data that were initially excluded from analysis can now be 
treated in the same way. That is, renormalization flow curves can be constructed as 
a function of A for sequentially increasing ni^^^^^. The results are shown in Figures 
I6.28l through l6.30[ Clearly, the renormalization flow curves for each plot correspond- 
ing to Co, C2 and C4 are flatter than those of the initial analysis, indicating a reduction 
in the regularization scale-dependence due to the use of data closer to the PCR. One 
is not able to extract an optimal regularization scale from these plots, as shown in 
the behaviour of X^^p displayed in Figures [6.311 through I6.33[ However, each x^^jr 
curve provides a lower bound for the regularization scale, where FRR breaks down 
UHLYIOII . as discussed in Section [6!2.1[ These lower bounds are: = 0.39 GeV, 

Aiower = 0-52 GeV and A^^,, = 0.59 GeV 

The statistical error bars of the low-energy coefficients corresponding to a small 
number of data points in Figures [6.281 through l6.30l are large, and a statistical differ- 
ence among the curves does not appear until ^^^^ ^ 0. 1 1 GeV^. Thus the identi- 
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Figure 6.28: Behaviour of co vs. A for the initially excluded low-energy data. A triple-dipole 
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Figure 6.29: Behaviour of C2 vs. A for the initially excluded low-energy data. A triple-dipole 
regulator is used. A few points are selected to indicate the general size of the statistical error bars. 
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Figure 6.30: Behaviour of 04 vs. A for the initially excluded low-energy data. A triple-dipole 
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fication of an optimal regularization scale will be aided by incorporating data corre- 
sponding to even larger values of j^^^^. By considering all of the available data, the 
behaviour of Xdof^ displayed in Figures [Oil through [6361 resolve precise optimal 
regularization scales: A^,^^^, = 0.72 GeV, A^'^^^^^, = 0.71 GeV and A^'^^^^^^^ = 0.71 
GeV. The systematic errors obtained from each x^^jr curve seem arbitrarily con- 
strained as a consequence of including more data points, which extend well outside 
the chiral regime, and possibly outside the applicable region of FRR techniques. 
This issue is addressed in the ensuing section. 



6.2.3 Optimal Pion Mass Region and Systematic Uncertainties 

In this section, a robust method for determining an optimal range of pion masses 
is presented. This range corresponds to an optimal number of simulation results to 
be used for fitting. First, consider the extrapolation of the quenched p meson mass, 
which can now be completed. The statistical uncertainties in the values of cq, C2, 
C4 are dependent on m^ ^^^^. As a consequence, the uncertainty in the extrapolated 
p meson mass must also be dependent on m^ jj^g^- Since the estimate of the 
statistical uncertainty in an extrapolated point will tend to decrease as more data are 
included in the fit, one might naively choose to use the largest m^ ^^sx. value possible 
in the data set. However, at some large value of m^ ^^^, FRR %EFT will not provide 
a valid model for obtaining a suitable fit. At this upper bound of applicability for 
FRR %EFT, the uncertainty in an extrapolated point is dominated by the systematic 
error in the underlying parameters. This is due to a greater scheme-dependence in 
extrapolations using data extending outside the PGR, meaning that the extrapolations 
are more sensitive to changes in the parameters of the loop integrals. Thus there is 
a balance point j^^^ = m^^^^x^ where the statistical and systematic uncertainties 
(added in quadrature) in an extrapolation are minimized. 

In order to obtain this value m\ j^^^, consider the behaviour of the extrapolation 
of the p meson mass to the physical point mp''g(m^ pj^yj, as a function of rn^^j^^x- 
Treating the parameters g2, g4, Mq, Aq and A^'^^''^ as independent, their systematic 
uncertainties from these sources are added in quadrature. In addition, the systematic 
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Figure 6.31: X^o/' '-o versus A, correspond- 
ing to the renormalization flow curves displayed in 
Figure 16.281 A lower bound for the regularization 
scale is found: A{°^^^ = 0.39 GeV. 
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Figure 6.32: x3o/' "^2 versus A, correspond- 
ing to the renormalization flow curves displayed in 
Figure 16.291 A lower bound for the regularization 
scale is found: A;I,„, = 0.52 GeV. 



0.35 




Figure 6.33: Xc/„/' for C4 versus A, correspond- 
ing to the renormalization flow curves displayed in 
Figure 16.301 A lower bound for the regularization 
scale is found: a;;^*^^^ = 0.59 GeV. 
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Figure 6.34: xL/' '-o versus A, correspond- 
ing to all available data, including the low-energy 
set. 
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Figure 6.35: X^i^f, for C2 versus A, correspond- 
ing to all available data, including the low-energy 
set. 
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Figure 6.36: x3o/' ^'^^ '-4 versus A, correspond- 
ing to all available data, including the low-energy 
set. 
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uncertainty due to the choice of the regulator functional form is roughly estimated 
by comparing the results using the double-dipole and the step function. These func- 
tional forms are the two most different forms of the various regulators considered, 
since the dipole was excluded due to the extra non-analytic contributions it intro- 
duces. The results for the initial and complete data sets are shown in Figures 16.371 
and l6.38[ respectively. Figure [637] indicates an optimal value m^ „iRx = 0.35 GeV^, 
which will be used in the final extrapolations, in order to check the results of this 
method with the low-energy data. By using only the data contained in the optimal 
pion mass region, constrained by max' estimate of the optimal regularization 
scale may be calculated with a more generous corresponding systematic uncertainty. 
The value A''^'"''' = 0.64 GeV is the average of Af^^"^^, All""^^ and A^.^'^''^ using this 
method. The Xdof analysis does not provide an upper or lower bound at this value 
of /w^ niax- These two estimates of the optimal regularization scale are consistent 
with each other. Both shall be used and compared in the final analysis. Figure [638] 
indicates an optimal value j^^^ = 0.20 GeV^ for the complete data set. A higher 
density of data in the low-energy region serves to decrease the statistical error esti- 
mate of extrapolations to the low-energy region. The corresponding value of A^'^^^'^ 
is unconstrained in this case, since the data lie close to the PGR. 

The values of cq, C2 and C4 for both the original data set and the complete data 
set are shown in Table l6.2l with statistical error estimate quoted first, and systematic 
uncertainty due to the parameters g2, g4, Mq, Aq, A^'^^^^ and the regulator functional 
form quoted second. In the case of the original data set, the value of C4 is not well 
determined, due to the small number of data points used. In the case of the complete 
data set, the results are dominated by statistical uncertainty, and this also results in 
an almost unconstrained value of C4. The coefficients of the complete set are less 
well-determined due to the fact that m^ „iRx = 0.20 GeV^, leaving only low-energy 
results with large statistical uncertainties for fitting. 

The result using the estimate of the optimal regularization scale A^'^^''^ = 0.64 
GeV, with the systematic uncertainty calculated by varying A across all suitable 
values, and using the initial data set, is shown in Figure 16.391 The extrapolation to 
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Figure 6.37: Behaviour of the extrapolation of the quenched p meson mass to the physical point 
'^P^q{'^k phys) '^n max using the initial data set, which excludes the lowest mass data points. In each 
case, Co is obtained using the scale Acentrai (for ^ triple-dipole regulator) as obtained from the X^^j- 
analysis. The error bars include the statistical and systematic uncertainties in cq added in quadrature. 
The optimal value ~ 0.35 GeV^. 
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Figure 6.38: Behaviour of the extrapolation of the quenched p meson mass to the physical point 
"^p^gl"^^ phys) "^nmax using the Complete data set, which includes the lowest mass data points. 
In each case, cq is obtained using the scale Acentrai (for a triple-dipole regulator) as obtained from 
the xff^^j- analysis. The error bars include the statistical and systematic uncertainties in cq added in 

quadrature. The optimal value max — 0-20 GeV^. 
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Table 6.2: The values of Co, C2 and C4 as obtained from both the original data set and the complete 
set, which includes the low-energy data. In each case, the coefficients are evaluated using the scale 
Acentrai (foi" a triple-dipole regulator) as obtained from the x^^^^ analysis. The value of used is 

that which yields the smallest error bar in adding statistical and systematic uncertainties in quadrature. 
For the initial data set, max — 0-35 GeV^. For the complete data set, m^max — 0-20 GeV^. The 
statistical uncertainty is quoted in the first pair of parentheses, and the systematic uncertainty is 
quoted in the second pair of parentheses. For the original data set, C4 is not well determined, with 
only a small number of data. For the complete data set, large statistical uncertainties result in an 
almost unconstrained value of 04. The coefficients of the complete set are less well-determined due 
to the fact that — 0.20 GeV^, leaving only low-energy results with large statistical uncertainties 

for fitting. 





co(GeV2) 


C2 


C4(GeV-2) 


original set 


1.31(5)(17) 


7.9(4) (26) 


-16.2(8)(382) 


complete set 


1.35(4)(241) 


6.8(5)(31) 


-3.3(17)(361) 



the physical point obtained for this quenched data set is: m^Q{m^ ^^^J = 0.922^Q QgQ 
GeV, an uncertainty of approximately 7%. Figure [6401 shows the data plotted with 
error bars correlated relative to the lightest data point in the original set, = 0. 143 
GeV^, using Agcaie = 0-64 GeV, and varying A across its full range of values. This 
naturally increases the estimate of the systematic uncertainty of the extrapolations, 
but also serves to demonstrate how closely the results from lattice QCD and xEFT 
match. 



6.3 Summary and Specific Issues for the Quenched p 
Meson 

A technique for isolating an optimal regularization scale was investigated in QQCD 
through an examination of the quenched p meson mass. The result is a successful 
extrapolation based on an extended %EFT procedure. By using quenched lattice 
QCD results that extended beyond the PCR, an optimal regularization scale was 
obtained from the renormalization flow of the low-energy coefficients cq, C2 and 
C4. The optimal scale is found to be A^'^^^'^ = 0.67lQ Qg GeV. An optimal value of 
the maximum pion mass used for fitting was also calculated, and was found to be 
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Figure 6.39: Comparison of chiral extrapolation predictions (blue diamond) with Kentucky Group 
data (red cross). Extrapolation is performed at Ajcaie = 0.64 GeV, varied across the whole range of 
A values, and using the optimal number of data points, corresponding to in^ j^^^ = 0-35 GeV^. The 
inner error bar on the extrapolation points represents purely the systematic error from parameters. 
The outer error bar represents the systematic and statistical error estimates added in quadrature. 
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Figure 6.40: Comparison of chiral extrapolation predictions (blue diamond) with Kentucky Group 
data (red cross), with errors correlated relative to the point at = 0.143 GeV^. Extrapolation is 
performed at Ascaie = 0.64 GeV, varied across the whole range of A values, and using the optimal 
number of data points, corresponding to m^y^^^ — 0.35 GeV'^. The error bar on the extrapolation 
points represents the systematic error only. A simple linear fit, on the optimal pion mass region, is 
included for comparison. 
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jjjax = 0.35 GeV^. By using only the data contained in the optimal pion mass 
region, constrained by jjj^, a value A^'^^^ = 0.64 GeV is estimated for the optimal 
regularization scale, with a wider systematic uncertainty corresponding to the entire 
range of values of A. These two estimates of the optimal regularization scale are 
consistent with each other. 

The mass of the p meson was calculated in the low-energy region. At the 
physical point, the result of the extrapolation, using A^'^^^^ = 0.67^o o8 GeV, is: 
m^^^(m2 pj^yj = 0.925+Jg5^ GeV. The result of the extrapolation, using A^*^^^ = 0.64 
GeV, with the systematic uncertainty calculated by varying A across all suitable val- 
ues, is: m-ym2 = 0.922^ 

ao60 GeV. The extrapolation also correctly predicts 
the low-energy curvature that was observed when the low-energy lattice simulation 
results were revealed. 

Since there exists no experimental value for the mass of a particle in the quenched 
approximation, this analysis demonstrates the ability of the technique to make pre- 
dictions without phenomenologically motivated bias. The results clearly indicate a 
successful procedure for using lattice QCD data outside the PGR to extrapolate an 
observable to the chiral regime. 
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Chapter 7 

Electromagnetic Properties of the 
Nucleon 

"[W]e can establish the key to our conclusion: the datum and the result are log- 
ically equivalent." (Omnes, R. 2002. Quantum Philosophy: Understanding and 
Interpreting Contemporary Science p. 209) [|Omn02ll 

In this chapter, the focus is turned to the magnetic moment and the electric charge 
radius of the nucleon. The magnetic moment is often studied for the physical sig- 
nificance of its anomalous component, obtained from the Pauh form factor F2 (de- 
fined in Equation (13.761) ). Since electrically charged pions with non-zero angular 
momentum dress the nucleon, they contribute non-trivially to its magnetic moment, 
altering the value from its semi-classical Dirac value. Likewise, the electric charge 
radius, or more precisely, the gradient of the Sachs electric form factor Ge in the 
soft-photon limit, provides a phenomenological test of quantum chromodynamics 
(QCD) theory. The leading-order low-energy contributions from virtual processes 
provide non-analytic behaviour in the chiral expansion. Chiral extrapolations for an 
infinite-volume box agree with experiment at the physical point, as will be evident 
later in this chapter. It is of interest in this investigation to determine if an optimal 
regularization scale may also be extracted from lattice QCD results for these two 
observables. If so, it would provide compelling evidence for the existence of an 
intrinsic scale for the source of the pion cloud of the nucleon. 



Ill 
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In lattice QCD, the isovector combination of the nucleon is often calculated, 
as described in Section I2.1.3I Feynman diagrams including any photons coupling 
to sea-quark loops cancel in the case of the isovector, and the distinction between 
VQCD and full QCD vanishes. This is fortunate, since the calculation of the dis- 
connected loops is computationally expensive. As a result, preliminary lattice QCD 
isovector results for two-flavor (a) -improved Wilson quark action from the QCDSF 
Collaboration are analyzed. 

The magnetic moment and the electric charge radius can each be written as chi- 
ral expansions, ordered in m^, due to the Gell-Mann— Oakes— Renner Relation from 
Equation (|3.43l) in Chapter [3l Each expansion comprises a polynomial residual se- 
ries, and loop integrals that contribute to non-analytic chiral behaviour. The dia- 
grams that correspond to the leading-order loop integrals are shown in Figures 17.11 
through O 

7.1 Renormalization of the Magnetic Moment 

7.1.1 Chiral Expansion of the Magnetic Moment 

Recalling the definition of the magnetic moment of the isovector nucleon in Equation 
(|3.80l) . the chiral expansion is as follows: 

= a^ + a^ml + T^{ml;A) + T^iml;A) + oK) , (7.1) 

for loop integrals denoted (T ) to differentiate them from the self energies. In this in- 
stance, only two free parameters are chosen, since the non-analytic contributions are 
included only to chiral order o(m^logmn). For a process with zero mass-splitting, 
such as that shown in the diagram in Figure mi the leading-order non-analytic term 
is proportional to m^; a lower chiral order than the leading-order term in the nucleon 
mass expansion. As a result, greater chiral curvature is expected, and the automatic 
renormalization process introduced in Chapter |4] will be constructed only to order 
(m^), that is, for the chiral coefficient cq. The fully renormalized chiral expansion 
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Figure 7.1: The pion/kaon loop contribution (with photon attachment) to the magnetic moment 
and the electric charge radius of an octet baryon B, allowing a transition to a baryon B'. All charge 
conserving transitions are implicit. 
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Figure 7.2: The pion/kaon loop contribution (with photon attachment) to the magnetic moment and 
the electric charge radius of an octet baryon B, allowing a transition to a nearby and strongly-coupled 
decuplet baryon T. 
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Figure 7.3: The tadpole contribution at 0{mg) (with photon attachment) to the electric charge 
radius of an octet baryon B. 

may be written to leading non- analytic order O (m^j) as: 

f/r = co + x''Nmn + 0{ml), (7.2) 

where /j is an implicit mass scale. Note also that the diagram in Figure 173] does not 
contribute to the magnetic moment of the nucleon since, in this case, the photon cou- 
ples to spinless pseudo-Goldstone bosons that have no orbital angular momentum. 
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7.1.2 Chiral Loop Integrals 

Each loop integral has a solution in the form of a polynomial expansion analytic in 
m\ plus non-analytic terms, of which the leading-order term is of particular interest: 

T;(m2;A) =Z^'^'^ + X^;m^ + M2''^'^m2 + o(m3), (7.3) 
T/(m2 ;A) = b^^""^^ + b^^'''^ ml + X^^ml\ogm^/ 0{ml) , (7.4) 

where is a mass scale associated with the logarithm. 

The corresponding loop integrals can be expressed in a convenient form by taking 
the non-relativistic heavy-baryon limit, and performing the pole integration for kQ. 
The integral corresponding to the diagram in Figure 17.11 takes the form [IWLTYOTl 
IWLTY09all : 

t;(™|;A) = =^ fa Wxj)VfrA) _^,... ^^^^ 

where q is the direction of the external momentum introduced by an incoming pho- 
ton. The argument for this substitution of the perpendicular part k^ is expounded 
in Appendix IB. 1.11 The function u{k;A) is the regulator, with associated momen- 
tum cut-off scale A. In this case, a dipole regulator will be used (corresponding 
to a choice of n = 1 in Equation (13.1051) in Chapter [3]). Since the working-order 
O (m^logmjt) of the calculation is less than in the case of the nucleon mass analysis, 
there is a reduced possibility of extra scale-dependent non-analytic terms frustrating 
the chiral fit. Thus, ensuring that these scale-dependent non-analytic terms are re- 
moved from the chiral expansion is not so vital, and a dipole form is an acceptable 
choice of regulator. The integral corresponding to the diagram in Figure 17.21 takes 
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the form: 

_ -xl r 3 J:^(2a)(fc)+A)^2(fc;A) _ a,a 

" 37i2 y'^'' 2(o3(^)[a)(fc)+A]2 ^0 ' ^'-^^ 

where (i){k) = A/^2"-i-mf ^^'^ ^ the mass-splitting. The chiral coefficients and 
%^ are constants in terms of the chiral Lagrangian of Equation (13.681) in Chapter [3j 



""871/2 9 - 5Ca • (/-ii) 



On the finite-volume lattice, each momentum component is quantized in units of 
In/L, that is, ki = ntln/L for integers n,-. Finite-volume corrections 5^'^^^ are writ- 
ten as the difference between the finite sum and the corresponding integral. It is 
known that the finite-volume corrections saturate to a fixed result for large values 
of regularization scale nHLYlOL As before, this is achieved in practice by evaluat- 
ing the finite-volume corrections with fixed regularization scale: A' = 2.0 GeV. The 
finite-volume version of Equation (17.11) can thus be expressed as: 

^„ = co + a^ml + {r^{ml;A) + Sf'^lm^; A')) 

+ {r^{ml;A)+df''^{ml;A')) + 0{ml) . (7.12) 



7.2 Evidence for an Intrinsic Scale in the Magnetic 
Moment 

The analysis of the magnetic moment of the nucleon provides an excellent check for 
the identification of an intrinsic scale in the nucleon-pion interaction. Using chiral 
effective field theory (xEFT), it has been demonstrated in Chapter[5]that lattice QCD 
results for the nucleon mass have an energy scale embedded within them. This prop- 
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erty is a consequence of the small size of the power-counting regime (PGR), where 
the expansion formulae of chiral perturbation theory (%PT) hold formally. Since 
a selection of lattice QCD results reasonable for fitting an extrapolation invariably 
extend outside the restrictive PGR [ILTY05|| . the validity of a formal scheme for ex- 
trapolation, and for identifying the leading-order terms in the chiral expansion, is 
compromised. Fortunately, a finite-range regularization (FRR) scheme, in conjunc- 
tion with %EFT as described in ChapterlH provides a robust method for achieving an 
extrapolation to physical quark masses, and identifying an intrinsic scale embedded 
within lattice QCD results. 

Recall that the method proceeds by analyzing the behaviour of the renormal- 
ization of one or more low-energy coefficients of the chiral expansion as a func- 
tion of the FRR scale. Ideally, that is, with lattice QCD results constrained entirely 
within the PCR, the renormalized coefficients should be independent of regulariza- 
tion scale. However, in practice, a scale-dependence is observed; particularly for 
data sets including data points corresponding to large quark masses. By truncating 
the lattice QCD results at different values of j^^^, an optimal FRR scale can be 
identified. This optimal scale is the value of A at which the low-energy coefficient 
under analysis is least sensitive to the truncation of the lattice data. If the optimal 
scale is consistent among the analyses of magnetic moment and the nucleon mass in 
Chapter |5l it provides evidence for an intrinsic scale in the nucleon. 

The preliminary QCDSF results for the magnetic moment at a variety of 
values are displayed in Figure 17.41 The experimental value is also marked. The 
set of data is listed in Appendix O Table IC.5I The lattice sizes of each data point 
vary from 1.43 to 3.04 fm using Nf = 2 and O (a) -improved Wilson quark action. 
A simple linear fit is included in this plot, which does not take into account the 
chiral loop integrals, nor the finite-volume corrections to the data. Therefore, it is 
not surprising that the linear fit fails to reach the experimental value of the magnetic 
moment at the physical pion mass. Since the lattice QCD results extend outside 
the PCR, the result of an extrapolation that includes the chiral loop integrals will 
be scale-dependent. However, the scale-dependence may be ameliorated using the 
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Figure 7.4: Preliminary lattice QCD data for jujf°'' from QCDSF, with the physical value from 
experiment as marked. 



procedure, which obtains both an optimal regularization scale and an estimate of its 
systematic uncertainty, constrained by the lattice results. 



7.2.1 Renormalization Flow Analysis 

In order to obtain the optimal regularization scale, the low-energy coefficient cq from 
Equation (17.121) will be calculated across a range of values of regularization scale A. 
Thus the renormalization flow can be constructed. Multiple renormalization flow 
curves may be obtained by constraining the fit window by a maximum, m\^^y^, and 
sequentially adding data points to extend further outside the PCR. The renormaliza- 
tion flow curves for a dipole regulator are plotted on the same set of axes in Figure 
17.51 Note that each data point plotted has an associated error bar, but for the sake 
of clarity only a few points are selected to indicate the general size of the statistical 
error bars. As more data are included in the fit, a greater degree of scale-dependence 
is observed. There is a reasonably well-defined value of A at which the renormal- 
ization of Co is least sensitive to the truncation of the data: A^'^''^'^ ^1.1 GeV. This 
indicates the optimal regularization scale embedded within the lattice QCD results 
themselves. 
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Figure 7.5: The renormalization flow of cq for obtained using a dipole regulator, based on 
lattice QCD data from QCDSF. 
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Figure 7.6: X^^/ for the renormalization flow of cq for obtained using a dipole regulator, 
based on lattice QCD data from QCDSF. 



7.2.2 Analysis of Systematic Uncertainties 

The optimal regularization scale for a dipole form can be more precisely extracted 
from Figure [73] using the chi-square-style analysis. Such an analysis will also pro- 
vide a measure of the systematic uncertainty in the optimal regularization scale. By 
plotting against the regularization scale A, where dof equals the number of 
curves n minus one for the fit parameter cq, a measure of the spread of the renor- 
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malization flow curves can be calculated, and the intersection point obtained. The 
'^dof P^*^*- corresponding to Figure 1731 is shown in Figure l7!6l The optimal regular- 
ization scale A^'^^^'^ is taken to be the central value A'^'^"^'^^' of the plot, and the upper 
and lower bounds obey the condition x^^^^ < x^^j: ^-^ + 1/ (dof). Thus the optimal 
regularization scale for a dipole regulator is: A^^^^^ = 1.13]']q 2o GeV. This value is 
consistent with the optimal regularization scale obtained for the nucleon mass using 
a dipole form, based on lattice QCD results in Chapter [5l Recall that the mean value 
for the optimal regularization scale from the nucleon mass analysis is: A^^p'® ^1.3 
GeV. This provides evidence that the optimal regularization scale is associated with 
an intrinsic scale characterizing the size of the nucleon, as probed by the pion. 

7.2.3 Chiral Extrapolation Results 

Using the optimal regularization scale, extrapolations or interpolations can be made 
to any quark mass. Consider the behaviour of the magnetic moment as a function of 
the quark mass as shown in Figure 1777] (in physical units). Here, the finite- volume 
expansion of Equation (17.121) is constrained by the lattice results from several differ- 
ent volumes. Extrapolation curves are then plotted for infinite volume and a variety 
of finite volumes at which current lattice QCD results are produced. For each curve, 
only the values for which m-nL > 3 are plotted, provisionally, to avoid undesired 
effects of the £-regime. The infinite-volume extrapolation to the physical point is 
within 2% of the experimentally derived value: = 4.6798/jyv. The finite-volume 
extrapolations are useful for estimating the result of a lattice QCD calculation at cer- 
tain box sizes. This can provide a benchmark for estimating the outcome of a lattice 
QCD simulation at larger and untested box sizes. Note that even a relatively stan- 
dard 3 fm lattice box length will differ significantly from the experimental value at 
the physical point. Since the data points in Figure 1777] are at differing finite volumes, 
the infinite- volume corrected data are also displayed in Figure ITIS] 
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Figure 7.7: Extrapolations of at different finite volumes and infinite volume, using a dipole 
regulator, based on lattice QCD data from QCDSF, lattice sizes: 1.43 — 3.04 fm. The provisional 
constraint m-nL > 3 is used. The physical value from experiment is marked. An estimate in the 
uncertainty in the extrapolation due to A^'^^''^ has been calculated from Figure ItTSI and is indicated at 
the physical value of m^. The curve corresponding to a lattice size of 10 fm is almost indistinguishable 
from the infinite-volume curve. 
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Figure 7.8: Extrapolations of at different finite volumes and infinite volume, using a dipole 
regulator, based on lattice QCD data from QCDSF, lattice sizes: 1.43 — 3.04 fm. The provisional 
constraint m-nL > 3 is used. The infinite-volume corrected data points are shown. The physical value 
from experiment is marked. An estimate in the uncertainty in the extrapolation due to A'*'^"''' has been 
calculated from Figure ItTSI and is indicated at the physical value of m^. The curve corresponding to 
a lattice size of 10 fm is almost indistinguishable from the infinite-volume curve. 
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7.3 Finite-Volume Considerations for the Electric Charge 
Radius 

Reliable extrapolations take into account finite-volume effects, as well as leading- 
order chiral loop corrections. In many cases, calculating the finite-volume correc- 
tions to loop integrals poses no essential problems. Examples of XEFT analyses 
accounting for finite- volume effects can be found in References llBea04b[ IHLYIOL 
However, the treatment of the electric charge radius is more challenging. Once 
form factors have been extracted from the lattice simulation, they are typically con- 
verted directly into charge radii. The essential difficulty lies in the definition of the 
charge 'radius' at finite volume. In order to define the radius, a derivative in the 
momentum transfer = — (at = 0) must be applied to the electric form 
factor. This approach breaks down on the lattice, where only discrete momentum 
values are allowed. 

In this chapter, a method is outlined for handling finite-volume corrections to a 
given lattice simulation result. It will be discovered that the finite-volume correc- 
tions to the loop integrals must be applied before the conversion from form factor 
to charge radius. By applying the finite- volume corrections directly to the electric 
form factor, and ensuring that the procedure preserves the electric charge normal- 
ization, an extrapolation in may be used to construct an infinite-volume charge 
radius. The infinite-volume charge radius can be defined as normal. A finite- volume 
charge radius may also be defined, as long as an allowed value of is used in the 
conversion from infinite to finite volume. 

The first challenge involves the definition of the electric charge radius in terms 
of this derivative in Equation (13.821) . Since only certain, discrete values of momenta 
are allowed on the lattice, the derivative may only be constructed from these allowed 
momenta when calculating finite-volume corrections. This crucial observation be- 
comes apparent when a comparison is made between the loop integrals evaluated at 
allowed, and unallowed, values of momentum transfer, respectively. The comparison 
is shown in Figures |7!9] and ITTTOl for momentum q = {2nn/L), on the lattice. Here, 
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Figure 7.9: Finite-volume correction for the loop integral contributing to Ge, with q = qmm- The 
choice of q/2 = <7rnin/2 is not an allowed value on the lattice. The momentum translated and untrans- 
lated behaviour of the finite-volume correction are inconsistent with each other 
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Figure 7.10: Finite-volume correction for the loop integral contributing to Ge, with q = 2qmm- 
The choice of q/2~ q^i„ is an allowed value on the lattice. Therefore, the momentum translated and 
untranslated behaviour of the finite-volume correction is identical. 
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n is an integer. Note that, if there is a momentum-translation in the loop integrals, 
k ^ k — q/2, the choice of ^ = qmin = {211 /L) (for box length L), means that q/2 
is no longer an allowed value on the lattice, and these finite-volume corrections will 
be inconsistent with the untranslated result. Under such a momentum-translation, 
external momenta of q/2 flow through the loop integral, and one should choose at 
least a value of ^ = 2^min to define a consistent discrete derivative for use in the 
definition of the charge radius in Equation (13.821) . However, choosing a momentum 
transfer of ^ = q^nm for ^ moderate lattice size of 3 fm leads to a relatively large 
value: ^ 700 MeV^. In defining the charge radius, the necessary extrapolation 
to 2^ = will be made more reliable by choosing a value of to be as small as 
possible. This situation differs from the infinite- volume calculation of loop integrals, 
where true momentum-translation invariance is restored, and a continuous derivative 
may be used as normal. 



7.3.1 Chiral Loop Integrals 

Though loop integrals in the continuum limit are invariant under momentum trans- 
lations k k + cq, c eTj (for internal loop momentum k), a finite- volume loop sum 
must not include any values of q less than ^min = {211/ L). Therefore, to obtain a 
suitable charge radius one chooses a definition of the loop integrals such that no fac- 
tors of q/2 appear. In fact, as long as no fractions of q appear in the integrand, the 
finite- volume version will converge correctly to the infinite- volume version as the 
box length is taken to infinity, for q = {2iin/L), n eZ: 

^i{Q^) = S U.^ ^^'-M)u{k-,K)u{k-q,K) 
371 y a)^a)^_.(a)^ + a)^_^^) 

^^^^ 3717 (a^, + A)(a>,_. + A)(a)^, + a>,_.)' 
^.f.(eVf/d^^^f^, (7.15) 
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where (0^ = \J~k?- — m^, and A is the mass-splitting. The chiral coefficients Xn ^"^^ 
%f and %f are summarized by Wang ||WLTY09all : 



167i2/^2 9 - 5cf'", (7.17) 
xf-^-^^-xf« (7.18) 

The integrals which contribute to the electric charge radius, denoted (T^), are ex- 
actly analogous to the integrals (T^) defined in Equations (17.131) through (17.151) . that 
correspond to the electric form factor Ge- To obtain the integrals that contribute to 
the charge radius, one simply takes the derivative with respect to momentum transfer 
2^ at vanishingly small values of Q^. This is allowed in the infinite- volume limit: 

T^= lim -6 ^Xo' . (7.19) 

Note that the ensuing procedure for calculating the finite-volume corrected electric 
charge radius uses only the infinite-volume versions of the chiral loop integrals. Fit- 
ting methods need only be applied at infinite volume. Thus, the external momentum 
derivative in Equation (17.191) need not be discretized, but may remain a continuous 
derivative. 

To achieve a chiral extrapolation, it is convenient to subtract the coefficients 
from the respective loop integrals that contribute to the electric charge radius: 



f^ = r^-ZpJ'^, (7.20) 

ff = rf (7.21) 

7;fd = 7;fd-^^- (7.22) 



This removes the regularization scale-dependence from the lowest-order fit parame- 
ter of the chiral expansion. This technique provides an advantage in easily extracting 
the low-energy coefficient cq from the chiral expansion, described in Section IT. 3. 2[ 
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As emphasized already, Figures fL9\ and ItTTOI show that the finite-volume correc- 
tions to the loop integrals cannot be applied directly to the charge radius itself. The 
momentum discretization ruins the circular symmetry in q except at the values coin- 
ciding with lattice momentum values [Inn/L), n eZ. The finite- volume corrections 
should be applied to the electric form factor Ge{Q^) instead. A momentum conven- 
tion in the loop integral is chosen such that q may be chosen to be ^min = {Inn/L). 
The procedure for achieving the correct finite- volume corrections is outlined below. 

First, the lattice finite-volume charge radius (r^)! must be converted into a finite- 
volume form factor G|(2^), using q = q^i^ = {2n/L). This is achieved through use 
of an extrapolation in Q^. As an example, a dipole Ansatz yields the following 
formula: ^ 

where Q^^j^ = ^^^^ — {Ef^ —M^Y. In many cases, this simply reverses the steps 
used to convert lattice results to charge radii. In this investigation, the electric form 
factor was fortunately obtained directly from the preliminary lattice QCD data from 
QCDSF. The next step is to transform the finite-volume form factor G^{Q^^^ to 
an infinite- volume form factor G'^{Q^i^, so that the infinite-volume charge radius 
can be calculated. This is achieved by subtracting the electric charge symmetry- 
preserving finite-volume correction, defined by: 

AL(eL,0) = 5l [t ^(eL) - ^""(O)] . (7.24) 

The second term of Equation (17.241) ensures that both infinite- and finite-volume 
form factors are correctly normalized, that is, G£°°(0) = 1 . Thus, the infinite- volume 
electric form factor can be calculated using the equation: 

G^eD = Gi(eL) - AL(eL,o). (7.25) 
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7.3.2 Chiral Expansion of the Electric Charge Radius 

The infinite- volume charge radius (r^)^ can be recovered from the form factor by 



tained, a chiral extrapolation can be performed if needed. The chiral loop integrals 
corresponding to the charge radius are those defined by Equation (17.191 ). Using the 
dipole Ansatz: 



This infinite-volume radius, calculated at multiple values of m^, can be used for 
fitting and obtaining coefficients from the chiral expansion: 



{r^E = {cT + 4^1) + T^i^l-^^) + Tl{ml;A) + fi^{ml;A) + {,4), (7.27) 



where the expansion has been renormalized in anticipation of the analysis of the 
renormalization flow of the coefficient cq. This expansion contains an analytic poly- 
nomial in plus the leading-order chiral loop integrals, from which non-analytic 
behaviour arises. 

By evaluating the loop integrals, the fully renormalized chiral expansion can be 
written in terms of a polynomial in and non-analytic terms. To leading non- 
analytic order o(logmn): 



Since the chiral expansion of Equation (I7.28h contains a logarithm, the value of cq 
can only be extracted relative to some mass scale ^, which is chosen to be 1 GeV. 

Finally, the finite-volume charge radius can be evaluated by adding the finite- 
volume correction to the form factor at any box length L, and corresponding mo- 
mentum transfer on the lattice, Qi, '. 



using the extrapolation in Q 



. Once the infinite-volume charge radius has been ob 




(7.26) 



(7.28) 



(7.29) 
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The finite- volume charge radii are obtained from the chosen extrapolation formula at 
box size L. An electric charge radius may be calculated at any desired value of box 
length, based on lattice QCD simulation results. Thus, the finite- volume behaviour 
of the charge radius may be analyzed. 

7.4 Evidence for an Intrinsic Scale in the Electric Charge 
Radius 

The preliminary QCDSF results for the electric charge radius of the nucleon are 
displayed, with the experimental value marked, in Figure 17.1 1[ The set of data is 
also listed in Appendix [Cj Table IC.61 The lattice sizes of each data point vary from 
1.92 to 3.25 fm using Nf = 2 and (a)-improved Wilson quark action. A simple 
linear fit is included in this plot, which does not take into account the non-analytic 
behaviour of the chiral loop integrals, nor the finite-volume corrections to the data. 
Just as for the case of the magnetic moment, the linear fit does not reach the ex- 
perimental value of the electric charge radius at the physical pion mass. Since the 
lattice QCD results extend outside the PCR, the result of an extrapolation will be 
scale-dependent. However, this scale-dependence can be handled by obtaining an 
optimal regularization scale using the aforementioned procedure. 

7.4.1 Renormalization Flow Analysis 

In order to obtain an optimal regularization scale, the low-energy coefficient c^q^ 
from Equation (17.271) will be calculated across a range of regularization scale values. 
Multiple renormalization flow curves may be obtained by constraining the fit win- 
dow by a maximum, m^ ^^^^, and sequentially adding data points to extend further 
outside the PCR. The renormalization flow curves for a dipole regulator are plotted 
on the same set of axes in Figure I7.12[ Note that each data point plotted has an 
associated error bar, but for the sake of clarity only a few points are selected to indi- 
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Figure 7.11: Preliminary lattice QCD data for (r^)!"" from QCDSF, with physical value from 
experiment as marked. 
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cate the general size of the statistical error bars. Note that, unlike the analysis of the 
nucleon mass and the magnetic moment, there is no distinct intersection point in the 
renormalization flow curves. In addition, the regularization scale-dependence of the 
coefficient Cq^'* is very slight, as long as the regularization scale is not too small, as 
discussed in Section 14.2.11 This lack of scale-dependence is a natural consequence 
of the logarithm in the chiral expansion of Equation (17.281) . which is slowly-varying 
with respect to the regularization scale. 

7.4.2 Analysis of Systematic Uncertainties 

An optimal regularization scale for a dipole form can nevertheless be extracted from 
Figure l77T2l using the chi-square-style analysis. The analysis also provides a measure 
of the systematic uncertainty in the optimal scale. By plotting x^^j- against the regu- 
larization scale A, where dof equals the number of curves n minus one, a measure of 
the spread of the renormalization flow curves can be calculated, and the intersection 
point obtained. The function is constructed in the same way as Equations (15.21) 
and (15.31) . The x^^^^ plot corresponding to Figure l77T2l is shown in Figure ITTSl Thus 
the optimal dipole regularization scale for a dipole regulator is: A'*'^^''^ = 1.67]']q33 
GeV. This value, though larger than optimal dipole regularization scale values ob- 
tained from the previous analyses of the nucleon mass and the magnetic moment, 
is nevertheless consistent, with one-standard-deviation agreement. Thus, strong ev- 
idence is found that the optimal regularization scale indicates the existence of an 
intrinsic scale, which characterizes the nucleon-pion interaction. 

7.4.3 Chiral Extrapolation Results 

Using the optimal regularization scale, a reliable chiral extrapolation can be per- 
formed, with the systematic uncertainty in the optimal regularization scale taken 
into account. Consider the behaviour of the electric charge radius as a function of 
the quark mass as shown in Figure 177141 (in physical units). Extrapolation curves 
are plotted for infinite-volume, and a variety of finite-volumes at which current lat- 
tice QCD results are produced. For each curve, only the values for which mj^L > 3 
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Figure 7.13: Xc/o/ for the renormalization flow of cq for obtained using a dipole regulator, 

based on preliminary lattice QCD data from QCDSF. cq is calculated relative to the energy scale /J = 1 
GeV. 

are plotted, provisionally, to avoid undesired effects of the £-regime. The infinite- 
volume extrapolation to the physical point differs from the experimentally derived 
value: (r^)^"^ = 0.88 fm^, by merely 0.5%. The finite-volume extrapolations are 
also useful for estimating the result of a lattice QCD calculation at certain box sizes. 
This can also provide a benchmark for estimating the outcome of a lattice QCD sim- 
ulation at larger and untested box sizes. Note that the result of an extrapolation to the 
physical point, using an optimistic 4 fm lattice box length, will differ significantly 
from the experimental value. Since the data points in Figure [TTH are at differing fi- 
nite volumes, the infinite- volume corrected data points are displayed in Figure ITTTSl 
To highlight the insensitivity of the extrapolation to the regularization scale A^'^^'®, 
an estimate of the systematic uncertainty in the extrapolation to the physical point 
solely due to A^'^^''^ is displayed in Figure 17.161 The size of the error bar at the 
physical point is comparable to that due to statistical uncertainty, as shown in Figure 
17.171 This indicates that, in the case of the electric charge radius, the identification 
of an intrinsic scale is borderline, due to the dominance of the logarithm in the chi- 
ral expansion, and its slowly varying property in the large m^j regime. Therefore, 
chiral extrapolations of the electric charge radius are more robust, in the sense that 
the scale-dependence in the result is suppressed, and the identification of an intrinsic 
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scale is not so vital as in the case of the nucleon mass or magnetic moment. 

7.5 Summary and Specific Issues for the Electromag- 
netic Properties of the Nucleon 

It was discovered that finite-volume corrections for charge radii are ill-defined on the 
lattice. The use of continuous derivatives in constructing the electric charge radius 
leads to inconsistent results for the finite- volume corrections. It was discovered that 
the finite-volume corrections must be applied to the electric form factors rather than 
to the charge radii directly. Therefore, a procedure was developed to apply finite- 
volume corrections to the electric form factor, strictly involving momenta available 
on the lattice. The resultant finite-volume corrected form factor may then be con- 
verted into a charge radius using an extrapolation in momentum transfer Q^. 

The technique for obtaining an optimal regularization scale from lattice QCD 
data has been investigated in the context of the magnetic moment and the electric 
charge radius of the isovector nucleon. By using recent, preliminary lattice QCD 
results from QCDSF, an optimal regularization scale for a dipole regulator was ob- 
tained. This was achieved, in each case, by analyzing the renormalization flow of the 
low-energy coefficient cq of the relevant chiral expansion with respect to the scale 
A, whilst extending the data step-wise beyond the PCR. A regularization scale was 
discovered, for both the magnetic moment and the electric charge radius, for which 
the renormalization of each cq is least sensitive to the truncation of the lattice QCD 
data. The values of the optimal regularization scale were consistent with each other, 
as well as with the results from the nucleon mass analysis. Thus an intrinsic scale 
has been uncovered, which characterizes the size of the nucleon, as probed by the 
pion. 

Using the value of the intrinsic scale, the extrapolation of the magnetic moment 
and the electric charge radius to the physical pion mass and the infinite-volume limit 
is consistent with experiment. The finite-volume extrapolations provide a bench- 
mark for estimating the outcome of a lattice QCD simulation at realistic or currently 
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Figure 7.14: Extrapolations of at different finite volumes and infinite volume, using a 

dipole regulator, based on preliminary lattice QCD data from QCDSF, lattice sizes: 1 .92 — 3.25 fm. 
The provisional constraint m-nL > 3 is used. The physical value from experiment is marked. The 
curve corresponding to a lattice size of 10 fm is almost indistinguishable from the infinite-volume 
curve. 
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Figure 7.15: Extrapolations of {r'^Y^°^ at different finite volumes and infinite volume, using a 
dipole regulator, based on preliminary lattice QCD data from QCDSF, lattice sizes: 1 .92 — 3.25 fm. 
The provisional constraint m-nL > 3 is used. The infinite-volume corrected data points are shown. 
The physical value from experiment is marked. The curve corresponding to a lattice size of 10 fm is 
almost indistinguishable from the infinite-volume curve. 
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Figure 7.16: Extrapolations of at different finite volumes and infinite volume, using a 

dipole regulator, based on preliminary lattice QCD data from QCDSF, lattice sizes: 1 .92 — 3.25 fm. 
The provisional constraint m„L > 3 is used. The infinite-volume corrected data points are shown. The 
physical value from experiment is marked. An estimate in the uncertainty in the extrapolation, due to 
^scaie^ has been calculated from Figure l7?T3l and is indicated at the physical value of m^. The curve 
corresponding to a lattice size of 10 fm is almost indistinguishable from the infinite-volume curve. 
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Figure 7.17: Extrapolations of at different finite volumes and infinite volume, using a 

dipole regulator, based on preliminary lattice QCD data from QCDSF, lattice sizes: 1 .92 — 3.25 fm. 
The provisional constraint mj^L > 3 is used. The physical value from experiment is marked. An 
estimate of the statistical uncertainty in the extrapolation is marked at the physical value of m^. The 
curve corresponding to a lattice size of 10 fm is almost indistinguishable from the infinite-volume 
curve. 
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optimistic lattice box sizes. 

The results clearly demonstrate a successful procedure for using lattice QCD 
data to extrapolate an observable to the low-energy region of QCD. 



Chapter 8 



Conclusion 

"Recall that in our theoretical construction those probabilities appeared simply as 
a logical, or linguistic, tool. It is only at this stage that they finally acquire the em- 
pirical significance they were lacking, and that chance enters the theoretical frame- 
work." (Omnes, R. 2002. Quantum Philosophy: Understanding and Interpreting 
Contemporary Science p. 209) IIOmn02ll 

8.1 Evaluation and Summary Analysis 

Chiral effective field theory (%EFT) offers unique insights into the low-energy be- 
haviour of hadrons. By using %EFT in conjunction with lattice quantum chromody- 
namics (lattice QCD) results, a deeper understanding of the underlying chiral inter- 
actions may be derived. In particular, the mathematical behaviour of the chiral ex- 
pansion of an observable, within a power-counting scheme (PGR), was investigated. 
This led to the development of a method for identifying the PGR, where the renor- 
malization of the low-energy coefficients of the chiral expansion are independent 
of the regularization scale. Novel methods for identifying a preferred renormaliza- 
tion scheme allowed the extrapolation of an observable to the chiral regime, and to 
infinite- volume lattice box sizes, without introducing a regularization scale in an ad 
hoc fashion. 

In this thesis, a procedure was established whereby an optimal regularization 
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scale could be obtained from lattice QCD data. By constructing some ideal pseudo- 
data using a known functional form, and based on actual lattice simulation results, 
the behaviour of the low-energy coefficients, with respect to the regularization scale, 
indicated an optimal value for the scale. By considering pseudodata sets that ex- 
tended increasingly beyond the PCR, there was a value of regularization scale at 
which the renormalization was least sensitive to this extension. This optimal scale is 
the value at which the correct values of the low-energy coefficients are recovered. 

Actual lattice simulation results for the nucleon mass, magnetic moment and 
electric charge radius were also analyzed using the same procedure. In each case, 
the analysis led to a consistent value of optimal regularization scale. In cases where 
multiple low-energy coefficients were analyzed, the optimal scale realized from each 
matched exactly: a non-trivial result. 

The analysis of lattice simulation results for the mass of the quenched p meson 
was used to test the robustness of the method. A reliable technique for determining 
an optimal regularization scale, and performing infinite- volume and chiral extrapo- 
lations, was established. 

Comparing the optimal scales obtained from the nucleon mass, magnetic mo- 
ment and electric charge radius analyses, a consistent optimal regularization scale 
was found. This indicates the existence of an intrinsic energy scale that character- 
izes the nucleon-pion interaction: the size of the nucleon as probed by the pion. 

In the analysis of the nucleon mass, as described in Chapter [51 it was demon- 
strated that a preferred regularization scheme exists only for data sets extending 
outside the PCR. However, it is not always possible to identify this scale. The scale- 
dependence of an observable can be weakened by working to a higher chiral order. 
The aforementioned procedure was used to calculate the nucleon mass at the phys- 
ical point, the low-energy coefficients cq and C2, and their associated statistical and 
systematic errors. Several different functional forms of regulator were considered, 
and lattice QCD data from JLQCD, PACS-CS and CP-PACS were used in the analy- 
ses. By working to chiral order O (m^), an optimal cut-off scale A^*^^'*^ for each set of 
lattice QCD data was obtained, and an estimate of the systematic error in the choice 



Chapter 8. Conclusion 



Hall 137 



of renonnalization scheme was calculated, using a chi-square-style analysis. A mean 
value for the optimal regularization scale of A^^^p^ ^1.3 GeV was obtained for the 
dipole, A^^JjJ^ ^1.0 GeV for the double-dipole and A^^^^ ^ 0.9 GeV for the triple- 
dipole. An analysis of the lowest suitable value for a regularization scale allowed 
the identification of a breakdown region of finite-range regularization (FRR). The 
existence of a breakdown region indicates that the ultraviolet regularization scale is 
low enough to remove or suppress the low-energy chiral behaviour being analyzed. 

The robustness of the procedure for determining an optimal regularization scale 
and performing chiral extrapolations was tested in Chapter [6l In order to estab- 
lish the predictive power of the procedure, the quenched p meson mass was consid- 
ered. Because an experimental value of this observable does not exist, its calculation 
served to demonstrate the ability of the procedure to make predictions without prior 
bias. Using lattice simulation results from the Kentucky Group, the procedure was 
tested, and the interesting low-energy simulation results were predicted correctly. 
By restricting the procedure to use only higher energy simulation data points, the 
low-energy coefficients cq, C2 and C4 were considered and an optimal regularization 
scale was identified: Ap'^^j}^ = O.GV^Q Qg GeV. An optimal value of the maximum pion 
mass used for fitting was also calculated, and was found to be ^^^^ = 0.35 GeV^. 
By using only the data contained in the optimal pion mass region, constrained by 
'^n.max' ^ value Ap'^l'jJfp = 0.64 GeV is estimated for the optimal regularization scale, 
with a wider systematic uncertainty corresponding to the entire range of suitable val- 
ues of A. These two estimates of the optimal regularization scale are consistent with 
each other. 

Upon revealing the omitted low-energy data, the extrapolations were compared 
to the simulation results at each value of pion mass. The correct chiral curvature 
was reproduced by the extrapolations, indicating the non-analytic chiral behaviour 
of the loop integrals. The results of extrapolations using xEFT, and the results of 
lattice QCD simulations were demonstrated to be consistent. The extrapolation to the 
physical point obtained for this quenched data set, using A^p^lf^ = 0.61^qqI GeV, is: 
"^p^y^^iphys) = ^-^^^-oxm GeV, an uncertainty of less than 6%. The result of the 
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extrapolation, using A!ff j^ = 0.64 GeV, with the systematic uncertainty calculated 
by varying A across all suitable values, is: /Wp^gl^^^ phys) — ^-^^^-o oto GeV, an 
uncertainty of only 7%. 

In the case of the electromagnetic properties of the nucleon, preliminary results 
from QCDSF were used. The magnetic moment of the isovector nucleon was an- 
alyzed for a dipole regulator. A well-defined optimal regularization scale was ob- 
tained: A^'^^lp = 1. 13^0 20 GeV, for chiral order o(m^logm7i), and a successful ex- 
trapolation to the physical pion mass and infinite- volume was achieved, and com- 
pared to the experimental value. The infinite-volume extrapolation to the physical 
point was within 2% of the experimentally derived value. 

When considering charge radii, there are subtleties in performing finite-volume 
corrections. In defining the charge radius, the finite-volume corrections must be 
applied before an extrapolation to 2^ = is taken. Thus the finite-volume correc- 
tions must be applied to the form factors directly. Using this method, the electric 
charge radius of the isovector nucleon was analyzed for a dipole regulator. Assum- 
ing the regularization scale is not within the breakdown region of FRR, the scale- 
dependence of the low-energy coefficient (up to some scale ^ of the chiral log- 
arithm) is weak. The leading-order non-analytic behaviour of the logarithm in the 
chiral expansion is slowly varying with respect to the regularization scale. Neverthe- 
less, an optimal regularization scale was obtained: A'^^^f^ = 1.671 GeV, working 
to chiral order o(m^logm,i). A successful extrapolation to the physical pion mass 
and infinite-volume was achieved, and compared to the experimental value. The 
infinite- volume extrapolation was merely 0.5% different from the experimentally 
derived value. 

Figure [STTI coUates the values of the intrinsic scale for a dipole regulator obtained 
from each of the three sets of lattice results from the nucleon mass analysis, the mag- 
netic moment analysis and the electric charge radius analysis. In summary, a method 
for determining the existence of a well defined intrinsic scale has been discovered. It 
has also been illustrated how its value can be determined from lattice QCD results. 
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Figure 8.1: Collated values for the intrinsic scale obtained from each data set for the nucleon mass, 
magnetic moment and electric charge radius, by analyzing a variety of low-energy coefficients. Each 
point, with its associated systematic error bar, is labelled by the low-energy coefficient analyzed. The 
results from the analyses of the nucleon mass are further denoted by the collaboration whose lattice 

/ -\ 

results are used. Cq and Cq denote the intrinsic scale obtained from the analysis of the low-energy 
coefficient cq corresponding to the magnetic moment, and the electric charge radius expansions, 
respectively. A dipole regulator is used. 



The research presented in this thesis encourages several avenues for further inves- 
tigation. In the heavy-baryon formulation of chiral perturbation theory (xPT), pre- 
sented first for the renormalization of the mass of the nucleon in Chapter HI the 
finite-volume corrections to the tadpole contribution are not evaluated. This is due 
to a technical subtlety associated with the large m-n behaviour of the finite-volume 
correction, due to the coefficient occurring in Equations (|4.16l) through (14.181) . 
The tadpole finite- volume corrections diverge as m^. Since it is known that the 
finite- volume corrections must converge [|Bea04bll . and lattice QCD simulations do 
not exhibit any divergence associated with large on a finite volume, higher-order 
terms, for example, those occurring at order (m^), must act to reduce the estimated 
value of the finite-volume correction. 

In the analysis of the renormalization flow of the nucleon mass, it was discovered 
that the scale dependence was weakened by working to a sufficiently high chiral 
order. It was also found, however, that the residual scale-dependence persisted as a 
significant component of the systematic uncertainty. For efficient propagation of this 
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uncertainty, an interesting future direction would be to consider Bayesian methods 
of marginalization over the scale-dependence [ISP09II . 

More generally, this research provides a strong basis for the investigation of 
baryon resonances by analysing lattice QCD simulations. Resonances of the nu- 
cleon, such as the Roper Resonance, are not well understood in terms of effective 
field theory. The structure and behaviour of the resonances lend themselves to a 
fruitful future area of research. Indeed, it is not possible to link the finite- volume re- 
sults of lattice QCD to experiment without understanding their relation to the multi- 
particle states that dress the resonances. FRR xEFT is particularly well-suited to 
exploring this important area of research. 

8.3 Codetta 

"[T]he collective efforts of numerous physicists have revealed some of nature's best- 
kept secrets. And once revealed, these explanatory gems have opened vistas on 
a world we thought we knew, but whose splendor we had not even come close to 
imagining." (Greene, B. 1999. The Elegant Universe p. 386) [|Gre99ll 

The dynamics of quantum chromodynamics provide a rich framework for the 
investigation of the properties of hadrons. In particular, low-energy effective field 
theory allows one to glean insights into the physical behaviour of subatomic parti- 
cles and the structure of matter. By incorporating the fundamental symmetries of 
quantum chromodynamics into the action, chiral perturbation theory provides a ro- 
bust method for the calculation of hadronic observables within the power-counting 
regime. In this thesis, finite-range regularized chiral effective field theory was used 
to develop a procedure for performing calculations beyond the power-counting regime, 
and handling any subsequent finite-range regularization scale-dependence. Using 
chiral effective field theory in conjunction with the non-perturbative approach of 
lattice quantum chromodynamics, chiral extrapolations; finite-volume effects; the 
consequences of dynamical chiral symmetry breaking on subatomic behaviour; the 
importance of strangeness; vacuum polarizations; and many other phenomena yield 
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fruitful understanding into the inner workings of the universe. 
Concluding Statement 

Chiral effective field theory allows the identification of an intrinsic energy scale 
in the nucleon-pion interaction from lattice simulation results. An optimal finite- 
range regularization scale, obtained from analyzing the renormalization flow of the 
low-energy coefficients of the chiral expansion, allows successful extrapolations to 
be made to the chiral regime and to the infinite-volume limit. There is strong evi- 
dence to suggest that the optimal scale characterizes the intrinsic energy scale of the 
interaction between the pion and the nucleon. 

The datum, the results and the rigorous theory integrate to form a strong ar- 
gument. Chiral effective field theory extended beyond the power counting-regime 
allows the identification of an intrinsic energy scale, and leads to a robust method 
for chiral and infinite-volume extrapolations. This is the original contribution of this 
thesis. 

"We have thus achieved the point where the theory may finally be compared with 
experience, and the road leading from formalism to concrete reality is at last com- 
plete." (Omnes, R. 2002. Quantum Philosophy: Understanding and Interpreting 
Contemporary Science p. 209) IIOmn02ll 
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Conventions 

A.l Dirac and Pauli Spin Matrices 

The Pauli matrices are usually chosen as such: 



= I ^ M (A.1) 

1 0, 



-i 

1 



t2=| I (A.2) 



1 
-1 



= I I (A.3) 



There are several conventions for the definition of the Dirac matrices (such as Weyl/Chiral 
or the Majorana Representation). Here is the Dirac Representation: 



7° = I " ^1 (A.4) 

-I, 



, .0 . 

i=\ (A.5) 
-& 



Y5 = 1^ ^1 = /T^yVy"- (A.6) 

I 
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All representations of these matrices satisfy the requirement of Clifford Algebra due 
to the conditions imposed in the derivation of the Dirac Equation ||Pes95L 

{Y',f} = 2g'\ (A.7) 
{y,y5} = 0. (A.8) 

A.2 SU(3) Gell-Mann Matrices 

The generators of the Lie Group SU(3) satisfy the commutator relations: 

=/r^^X^ (A.9) 

This, combined with the relevant Jacobi Indentities for the generators, defines the 
structure constants [|Pes95ll : 



j-ade J'b'^d _|_ j-bde j-cad _|_ J''^de j-abd j^q-j 



A.3 Spinor Fields 

The equal-time canonical anti-commutation relations for Dirac spinor fields are: 

{\\f{x) , \\f{y) }^^=y^ = hd^ (x-y), (A. 11) 

Mx)My)U=yo = 0. (A.12) 

The fields take the form [|Pes95ll : 



and the canonical anti-commutation relations expressed in terms of the Pauli- Jordan 
function: 

{\\f{x),^{y)} = {0^ + m)iA{x-y;m). (A.14) 
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The Grassmann algebra is defined by the anticommutation rule between Grass- 
mann variables V|/ and a commutation rule with non-Grassmann numbers c: 



(A.15) 



For Berezin integration over fermion spinor fields \\f and \\f, the follow rules are 
adopted: 



J d\\fi\\fj = J d\\fi\\fj = cdij , 

f 



(A.16) 
(A.17) 



where the non-Grassmann constant c is chosen, by convention, to be equal to 1 and 
the function / is defined on the Grassmann algebra. As a consequence of Equation 
(IA.17I) . the Berezin integral over unity vanishes: 



y^d\|/= Jd\\f = 0. 



(A. 18) 



A.4 Meson and Baryon Field Definitions 

The SU(3) mixed- symmetric meson octet fields n{x) — 'K,"{x)'k" can be encoded in a 
traceless 3x3 matrix of the form: 



nix) 



V2 



K 



\ 



\ 



K- 



(A.19) 



In SU(2) the pions form the triplet representation {n , Ti:*^, 7r+) which can be written 
by summing over the Pauli spin matrices in Appendix (IA.1|) : 



7r(;c)=x"7r"(jc) 



71° V2TI+ 
V2 -71° 



(A.20) 
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Using the convention for Clebsch-Gordan coefficients from Wang et al. HWLTY09all . 
the the mixed-symmetric baryon octet matrix has the form: 



B{x) 



E+ 

LyO , X A 



(A.21) 



The maximally symmetric decuplet tensor (suppressing Lorentz indices) has ele- 
ments defined by: 



riii=A^ 



7^1 12 = yl^^' ^122 



1 



A", Till = A 



1 



*,0 



1 



(A.22) 
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Integration Techniques 



B.l Magnetic Quantities 



B.l.l Angular Components of Magnetic Moment Loop Integrals 

In anticipation of applying finite-volume corrections to chiral loop integrals by com- 
paring them to their respective summations on the lattice, the time-component of the 
d^^ integral is evaluated using Cauchy's Integral Formula, and a d^k integral remains 
for analysis, as in Chapters |4] through [6l 

When calculating the magnetic moment in the heavy-baryon limit, without ex- 
plicitly specifying a regularization scheme, the one-loop integral (corresponding to 
Figure ItTI) takes the following form: 



It is useful to be able to simplify the angular part of the integral, formed by the 



a numerical coefficient. In order for this to be valid in calculating finite-volume 
corrections, the simplification must hold in both integral and sum forms of the loop 




(B.l) 



cross product of external momentum direction q with the loop momentum k, into 
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diagram. Evaluating the angular part of Equation (IB.ll) yields: 

<r; = d<p dk dx "^^ ^j, (B.2) 

^-^iTdkl'd^^ (B.3) 
% Jo J -I {k^ + mlf 

= dk J. (B.4) 

Now, this one-dimensional integral can be transformed into a three-dimensional 
integral simply by adding in a naive solid angle component, using the identity: 

371 An J ""Jo ""(k^ + ml) 



r/: = -^- d^ dk----^ (B.5) 



d'k— ^ . (B.6) 



37l2y {k^ + ml] 



Comparing Equations (IB.ll) and (IB. 61) shows that the objective has been achieved for 
the integral case. For finite volume sums, the result may not hold in general, and so 
must be checked independently. Define the following sum for box length L: 

Because k^ = k^ is symmetrical in directions ^xj,z> it follows that: 



1? = k^+k^ + k^ = 3kf {i = x,y,z), (B.8) 
and since (qxk) = k\, (B.9) 
it follows kl = 2kf. (B.IO) 
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Thus: 



which is the finite- volume equivalent of Equation (IB.6I) . 



B.1.2 Combinatorial Simplification 

The calculation of the three-dimensional finite sum can be made more efficient com- 
putationally, by transforming it to a one-dimensional sum in terms of the new vari- 
able = k^{2ii/L)^. It does, however, require calculation of the number of config- 
urations of the squares of kx, ky and k^ to obtain each value of n^, denoted C^^\n-^). 
Thus, for an integrand I (k) : 




(B.13) 



where n^^ = /r^^(27c/L)2. 

B.1.3 Sachs Magnetic Form Factors at Finite 

Consider calculations involving the leading-order pion loop contributions to the 
magnetic form factor Gm{Q^) at finite Q^, (allowing non-zero mass splitting A). 
The following integral can be made more easily calculable using spherical polar co- 
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ordinates (using (o(k) = ^/k"^ + m^y. 



k^[(£)(k)+(£)(k + q)+A] 



(ii{k) [co(fc) + A] (i){k + q) [a)(fc + ^) + A] [&{k) + &{k + q)] 
' dk 



dcp/de 

Jo Jo Jo 



k^ k^ sin [(xi{k) +(xi{k + q)+A] 



(0(fc) [(f){k) + A] Oi{k + q) [(Ji{k + ^) + A] [(o(fc) + (o(^ + ^)] J 



(B.14) 



(B.15) 



The integral can be further altered to remove the infinite integral under the change 
of variables k^k/{l—k). For arbitrary function /: 



rdkm = [ 

Jo Jo 



\. f{k/{\-k)) 

(1-^)2 



(B.16) 



10 JO 

Thus, defining (for convenience) p{k) =k/{l—k). Equation (IB.15I) becomes: 



Tl^ Jo 



1 rl 

dx dk 
1 Jo 

a ^2 



k]p^ [(ii{p)+(ii{p + q)+A] 



(d{p) [qo{p) + A] &{p + q) [(0(p + ^) + A] [<ji{p) + (0(;? + q)]{\-kf 



(B.17) 



B.2 Electric Charge Radius Integral Expansions 

For the infinite- volume electric charge radius, the chiral loop integrals must be calcu- 
lated for use with the chiral expansion of Equation (17.271) . Each loop integrand is ex- 
panded out for small Q^, and the derivative in the limit of vanishing is extracted. 
Using the notation of Chapter |71 and a dipole regulator, the one-loop contribution 



Appendix B. Integration Techniques 



Hall 151 



takes the following form: 



6%f /^3 



57C J ^ ^dq^ 



{k + q/2) ■ {k-q/2)uiip(k + q/2;A) Udip{k-q/2;A) 
Gi(k + q/2)0)(k - q/2) [(0(X + q/2) + a)(X - q/2)] 



q2=0 

(B.18) 



6% 



571 



Jd^k (a^ { -0)2 (k) (P + ) ( 1 + 2mlA^ + 5P (2m^ + A^ ) ) 



+ ^"^(21^^ + 16m^ + 5A'^ + 2P(16m^ + 5A2)) cos^ 0} j (^16co^(^)(P + A^)^^ ^ . 

(B.19) 



If a mass-splitting is included: 



571 J 



3^2 



(/; + .7/2) • - .7/2) »dip(^ + ^7/2 ; A) »aip(^ - ^7/2 ; A) 
{(0{k + q/2) +A){0){k- q/2) + A) [0){k + ^/2) + to(^ - q/2)] 

(B.20) 



6X| 

571 



J dh(A^ {-(>}{%) (J? + A^) [13P + 2ml{ml + A{2(o{k) + A))A^ 



+ F(23m^ + 24(£){k)A+ 1 1 + SA^) 

+ P(10m^ + A(8co(^) +3A)A2+m^(20co(^)A+ lOA^ + VA^))] 

+ F[21fc^ + 16m^ + 16m^A(2to(^) + A) +5m^A'^ + A(4to(^) + A)A'^ 

+ k\53ml + 36&{k)A + IVA^ + lOA^) 

+ P(48m^ + a2(8co(^)A + 2A^ + 5A^) + 2m^(32co(fc)A + leA^ + 5A^))] cos^ 6}) 



X {l6(>y'{k){&{k)+A)\'P + A 



.2x6 



Similarly, the tadpole contribution takes the following form: 



(B.21) 



'tad 



7t J 



d^k 



dq^ 



(0{k + q/2)+(jii{k-q/2) 



q2=0 



7C 



16(05(yt) 



(B.22) 
(B.23) 
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B.3 Finite Volume Corrections to Tadpole Amplitudes 

Finite-volume corrections should vanish as becomes large, as observed in lattice 
quantum chromodynamics (lattice QCD) simpulations. This has also been observed, 
in turn, for each of the finite- volume corrections involved in the extrapolation of the 
nucleon mass to fourth-order. However, the tadpole finite- volume correction, ^J^f', 
is different in that it is multiplied by a factor of m^, as evident in Equation (|4.7I) . 
The product C2m\^^^ is not convergent for large m^. Figures IBTI and |R2] show the 
behaviour of the tadpole finite- volume correction for a 2.9 fm box and a 4.0 fm box, 
respectively. 

The finite-volume estimate of C2, denoted c\, is not in general the same value 
as the infinite- volume C2- Thus the finite- volume correction of the tadpole cannot 
be written as simply the difference between the finite volume sum and the infinite 
volume integral, but must distinguish between and C2'. 



cimldfj = ciml (^^EL - ^tadj . (B.24) 

Since C2 is by definition the coefficient of the term in the nucleon mass expan- 
sion, the renormalization of the residual coefficient a2 by the contributions from the 



integrals Eyy, and T^tad, defined in Equations (14.31) through (14.81) . can be written as 
follows: 



C2mi = {a2 + b^+b^ + C2b'2)m^ , (B.25) 

=»C2= ' ' (B.26) 
l-b'2 



An analogous relation exists for the finite volume : 



, iV,N , iV,A 

./ a2 + b^ +b^ 
4= ^\y, ^ • (B.27) 



By simultaneously solving for a2 and c^, the ratio c^/c2 can be calculated in prin- 
ciple, and the tadpole finite volume corrections are tractable. It should be noted 
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Figure B.l: Behaviour of the finite- volume corrections ^^^f' vs. A on a 2.9 fm box 
using a dipole regulator. Results for two different values of are shown. 
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Figure B.2: Behaviour of finite-volume corrections d^^f' vs. A on a 4.0 fm box 
using a dipole regulator. Results for two different values of are shown. 
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however, that this does not resolve the problem of divergent behaviour for large mvt. 
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Table C. 1 : JLQCD llO+OSl lattice QCD simulation results for the nucleon mass Mn at various pion 
mass squared values m^. The lattice spacing is 0. 118 fm and the spatial lattice length is 1 .90 fm. 



m2(GeV2) 


M^(GeV) 




0.567 


1.615(6) 


7.25 


0.386 


1.456(6) 


5.98 


0.273 


1.350(6) 


5.03 


0.191 


1.255(6) 


4.20 


0.135 


1.164(8) 


3.54 


0.084 


1.111(10) 


2.78 



Table C.2: PACS-CS lA+091 lattice QCD simulation results for the nucleon mass Mn at various 
pion mass squared values in^. The lattice spacing is 0.0907 fm and the spatial lattice length is 2.90 
fm. 



m2(GeV2) 


M^(GeV) 




0.492 


1.583(5) 


10.32 


0.325 


1.411(12) 


8.38 


0.169 


1.215(12) 


6.05 


0.087 


1.093(19) 


4.35 


0.024 


0.932(78) 


2.29 



Table C.3: CP-PACS fAK+02l lattice QCD simulation results for the nucleon mass M^, the lattice 
spacing a and the spatial lattice length L at various pion mass squared values m^. 



m2(GeV2) 


Miv(GeV) 


a(fm) 


L(fm) 


ninL 


0.940 


1.809(15) 


0.102 


2.45 


12.03 


0.913 


1.798(4) 


0.130 


3.12 


15.11 


0.704 


1.652(9) 


0.099 


2.38 


10.10 


0.689 


1.643(5) 


0.123 


2.95 


12.42 


0.539 


1.519(9) 


0.095 


2.28 


8.49 


0.502 


1.497(6) 


0.118 


2.83 


10.17 


0.353 


1.348(12) 


0.092 


2.21 


6.65 


0.272 


1.275(7) 


0.111 


2.66 


7.04 



Appendix C. Lattice Simulation Results 



Hall 157 



Table C.4: Quenched lattice QCD data for the p meson mass /Wp at various pion mass squared 
values m^. The lattice size is 20^ x 32, with a lattice spacing of 0.153 fm. 



m2(GeV2) 


mp(GeV) 




3.150 


2.001(1) 


27.53 


2.187 


1.700(2) 


22.94 


1.742 


1.548(2) 


20.47 


1.329 


1.399(2) 


17.88 


1.212 


1.354(2) 


17.08 


1.062 


1.294(2) 


15.98 


0.867 


1.214(3) 


UAA 


0.743 


1.162(4) 


\331 


0.676 


1.133(4) 


12.75 


0.610 


1.103(5) 


12.12 


0.515 


1.060(5) 


11.13 


0.422 


1.016(6) 


10.07 


0.347 


0.985(7) 


9.13 


0.288 


0.960(8) 


8.32 


0.241 


0.938(8) 


7.62 


0.204 


0.926(9) 


7.00 


0.172 


0.914(11) 


6.43 


0.143 


0.908(14) 


5.87 


0.114 


0.899(15) 


5.24 


0.094 


0.899(16) 


4.75 


0.080 


0.896(18) 


4.38 


0.068 


0.898(20) 


4.04 


0.059 


0.902(22) 


3.77 


0.053 


0.903(26) 


3.58 


0.047 


0.907(28) 


3.37 


0.041 


0.913(32) 


3.15 
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Table C.5: Preliminary lattice QCD simulation results from QCDSF for the isovector nucleon 
magnetic moment /yJ^, the lattice spacing a and the spatial lattice length L at various pion mass squared 
values OT^. 



m2(GeV2) 




fl(fm) 


L(fm) 




0.863 


2.394(69) 


0.089 


1.43 


6.73 


0.709 


2.483(45) 


0.073 


1.76 


7.50 


0.688 


2.548(159) 


0.091 


1.45 


6.11 


0.591 


2.621(49) 


0.084 


2.01 


7.85 


0.392 


2.863(86) 


0.070 


1.67 


5.30 


0.357 


2.781(51) 


0.084 


2.03 


6.13 


0.290 


2.840(121) 


0.070 


1.67 


4.57 


0.198 


3.082(120) 


0.081 


1.96 


4.42 


0.159 


3.006(118) 


0.077 


1.84 


3.72 


0.077 


3.711(158) 


0.076 


3.04 


4.26 



Table C.6: Preliminary lattice QCD simulation results from QCDSF for the isovector nucleon 
electric charge radius {r^)E, the lattice spacing a and the spatial lattice length L at various pion mass 
squared values m^. 



m2(GeV2) 


(r2)£(fm2) 


fl(fm) 


L(fm) 


rriTiL 


0.591 


0.303(8) 


0.084 


2.01 


7.85 


0.357 


0.349(6) 


0.084 


2.03 


6.13 


0.349 


0.340(5) 


0.080 


1.92 


5.75 


0.198 


0.384(11) 


0.081 


1.96 


4.42 


0.188 


0.392(12) 


0.068 


2.19 


4.81 


0.074 


0.494(25) 


0.076 


3.04 


4.18 


0.053 


0.586(24) 


0.068 


3.25 


3.79 
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